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CXXIII. Some Properties of Lines on Divergent-beam X -ray Photographs 


By A. G. Peace and G. E. Prrncie 
Department of Inorganic and Physical Chemistry, University of Leeds* 


[Received August 25, 1952] 


SUMMARY 


The quantitative aspect of x-ray extinction is reconsidered in its 
particular relation to divergent-beam photography. Special consideration 
is given to the optical resolution obtainable with the usual experimental 
arrangements, and some experimental data on unresolved lines are 
discussed. This leads to a determination of the mosaic spread of crystals 
previously conditioned by liquid air treatment to minimize primary 
extinction effects, but in general quantitative conclusions are hampered 
by the usual doubt as to how much of the observed extinction is primary 
and how much secondary. 


List oF SYMBOLS 
a=spacing of reflecting planes. 


A=/(A’'A”"). 
A’=qxnumber of reflecting planes intersected by a line parallel 
to the reflected beam. 
A”=qx< number of reflecting planes intersected by incident ray. 
c=ratio in which absorption is diminished by Borrmann effect. 
C,=Euler’s constant. 
C,=a constant, approximately 0-25. 
G(w)=the partition of Q as a function of u. 
G’(u)=the partition of Q’ as a function of w. 
G,(u)=G(u) modified by Borrmann effect. 
I,=Bessel function of imaginary argument, order zero. 
I,=Bessel function of imaginary argument, order 1. 
J —Bessel function of order zero. 
J,—Bessel function of order lI. 
K=ratio of white radiation intensity to characteristic. 
m=variable number of multiple reflections. 
qg=amplitude of beam reflected by one plane when the incident 
beam has unit amplitude. 
Q=volume reflecting power for given set of planes. 
Q’=Q diminished by primary extinction. 
#(x)=intensity of reflected beam at position x in direction of 
incidence. 


* Communicated by Professor E. G. Cox. 
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o¥(z)=function defined by Darwin, expressing the integrated 
reflection by transmission through a mosaic (§ 3). 
TJ («)=intensity of emergent beam at position x. 
t=depth in crystal. 
t’=depth in crystal measured obliquely in direction of reflection. 
t”"=depth in crystal measured obliquely in direction of incidence. 
T, T’, T"=final values of t, t’, t’ respectively at exit surface. 
«x=fractional coordinate of point of exit. 
y=(1—2). 
z—argument of functions ec? and J). 
B=factor expressing source size. 


5 27 wa cos Af sin ay 
fee ear i o+6/° 
u+G@’ up 
PEEL PSU eU ET 
A=wavelength of characteristic radiation. 
n=Mmosaic spread in given plane. 
6= Bragg angle. 
§=angular half width of divergent-beam line for monochromatic 
radiation. 
¢é=inclination of reflecting planes to crystal surface. 
p=ordinary absorption coefficient of characteristic radiation. 
o=contrast of divergent-beam line against neighbouring back- 
ground. 
o,,—=contrast in absence of white radiation. 
€=diffraction width. 
Y=absorption ratio for white radiation/absorption ratio for 
characteristic. 
w(w)=weighting function expressing source shape. 


§1. INTRODUCTION 
Meruops of producing transmission photographs of crystalline specimens 
with wide divergent x-ray beams have been described by various authors 
(e.g. Borrmann 1941, Lonsdale 1943, May 1949). Amongst many 
physically important applications is the study of extinction effects and 
the textural investigation of single crystals (e.g. Lonsdale 1947 a and b). 
One of the differences from parallel-beam photography in practice is 
that divergent-beam photographs can often be obtained in a matter of 
a few seconds and this recommends them as a means of sorting a batch 
of crystals into grades of different microscopic texture. In deciding 
whether the information so used can be pursued to quantitative 
conclusions, we have first to consider the adequacy of extinction theory 
and its relevance to the divergent-beam arrangement. On this basis 
(§ 2) we can attempt to correlate extinction coefficients with measurements 
of line width (§3). The further step (§5) of converting results of 
measurement into knowledge of primary extinction factors and of block 
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dimensions is hampered by uncertainty of the theoretical factors for 


finite blocks, and by the difficulty of separating primary extinction 
effects from secondary. 


§2. RESOLUTION AND LINE STRUCTURE 


In Darwin’s paper (1922) the relevant theory is presented for parallel 
beams of x-rays which are unlimited in lateral extent. The arrangement 
for producing divergent-beam photographs, such as the one we have 


Fig. 1 
\ Film / 
\ Z 


filament 
assembly, 


Siemens tube fitted with transmission target as used for divergent-beam 
_ photography. 


used (fig. 1) depends for its success on employing a small X-ray source, 
smaller than the slab of crystal covering it, so that the simple theory can’ 
be supplemented for this case by studying the behaviour of a narrow’ 
beam of x-rays passing through a single crystal, assumed in what follows 


to be a plate with parallel faces. Without recourse to differential 
equations we can see from the theory of secondary extinction (Darwin 


loc. cit.) that, ignoring rays reflected back into the beam, the apparent 


= 


attenuation of a given monochromatic primary ray incident on a mosaic 


-erystal at an angle differing by a small amount w from a reflecting angle is 


exp [— {p--G/(u)yie, 4 ee wel) 
402 
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where y is the ordinary true linear absorption coefficient, ¢’ is the variable 
path length in the crystal at depth ¢ and G’(u) is the partition of Q’ the 
volume reflection coefficient for the planes in question, i.e. 


G’(u)=G'(0) exp (—0')2n)= a exp (—w?/2n7). . . (2.2) 


The argument of G’ will usually be omitted for brevity in what follows. 
Corresponding unprimed quantities G and @ refer to the ideal mosaic 
conditions with negligible primary extinction. The exponential expres- 
sion (2.1) can be regarded as a ‘ probability of survival’ of a packet of 
energy. Pursuing a statistical argument, we equate the probability of 
reflection per unit length of path to G’(w); then the probability of 
arrival at the exit surface after an even number of reflections, m, can 
easily be shown to be, in the ‘ symmetric Laue case ’ 

(ie 

m | 


exp [—(u+@')7] 


where 7” is the final value of t’ at the exit surface. For the moment we 
confine our discussion to this case, which is mathematically the simplest 
because exp (—Z") appears only as a factor throughout. The primary 
beam together with all even reflections give a total of 


exp [—(u+G’)T’] cosh G’T" 


whereas the ‘ reflected ’ beam, consisting of the sum of the odd reflections 
gives exp[—(u+G")T"] sinh G’T” ; total exp(—pZ"). In a certain sense 
therefore, the two beams are complementary, and the power deficiency 
of the transmitted beam will be equal to the reflected power if the 
background proportional to exp (—7") is taken as the level of reference. 

In order to decide whether the constituent terms comprising these two 
beams are, experimentally, added together, we need to know the lateral 
spread of the reflections in passing through the crystal. We find that 
the beam does tend to spread laterally, but much more slowly than the 
first power of the depth in the slab. This conclusion follows from a 
statistical argument similar to that just used. Adapting this argument 
to the geometry of fig. 2 where x and y (=1—2) are fractional co-ordinates 
of the point of exit in the range BC, we find that the intensities of multiple 
reflections issuing at the point x when unit power is incident in the form 
of a sufficiently narrow beam on the first surface can be expressed for 
the general Laue case as 


> 


(a) Mosaic crystal 
Transmitted, TF (x)=exp [—(u+G@’)(xT"+yT")] 
Y hid 
q’ Ri (Sr )iL26" Vey") cosec 20. 


Reflected, RB(x)=exp [—(u+G@')(aT’+yT")] 
GT [2G4'/(ayT"’T")] cosec 26. 


(2.3) 


re 
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(b) Perfect non-absorbing crystal; in this case, incident intensity x 
(width) is equated to unity. 


Transmitted, Ge (=| g <i (=F) J Rv(qya'a yy} 


2a cos 0 


Reflected, Mold [2./(ayA’A” : ie 

2a cos@ ° 4 ni} ‘ 

where 7” and 7” are the two extreme path lengths in the crystal : 
T’=T cosec (¢+6), 


TIT” =T cosec (6—6), 


Fig. 2 
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RAY 
Multiple x-ray reflections in thick crystal mosaic. 


where 7’ is the true thickness of the crystal plate, and similarly A’, A” 
are the numbers of reflecting planes intersected by 7” and 7”, multiplied 
by q, the reflecting power of one plane ; if a is the spacing of the reflecting 
planes 

A'=qT" sin 0/a, 

A"=qT" sin 6/a. 


J, and I, are Bessel functions of zero order of real and imaginary argument ; 
J, and I, are of the first order. In these expressions we assume no 
limit to the maximum possible number of secondary reflections: if this 
were done we should have to curtail the power-series development of 


2.3 and 2.4 at the appropriate stage. 
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By integrating the reflected amplitude for the perfect crystal over the 
range of values of x and making allowance for phase, we recover Waller’s 
formula (1926) for the power of the reflected beam 


sin? A 4/(1+8?)/1-+6?, 
| where A=A/(Ara~) 
27 ua cos 0 sin d—6 
oa a ene 


The corresponding formula for the mosaic crystal is 
4p Sig pls ye g 
exp [—w+e)(=S*) | /(q) sinh Bvt evate) 


for the reflected power where 


B= GAs CAS: 


Oe piel et 
c= (Eu 24/ (TETS) ? 
which reduces to Darwin’s formula in the symmetric case (e=0) ; in the 
slightly asymmetric case we obtain a fair approximation by 


5 | (a) 0—exp Lory" +7") exp | (AS) } 


neglecting «” and higher powers. 

Formulae 2.3 and 2.4 enable us to find where the reflected beam issues 
from the crystal. In accordance with Borrmann’s observations (1950) 
on calcite, the beam tends to issue, in the symmetric case, at the mid- 
point of the range as if guided by the crystal planes ; however, in the 
asymmetric case the corresponding result is not obvious. If G’ is large 
we may write 


I,(z) = 
so that Z(x) contains the factors 
UV (eyT’T") exp [—p(xT"+yT")] exp {—(@)LV(eT)— VP") P} 5 


‘if therefore 4 exerts only a minor influence, the maximum reflection will 
‘issue where 


bal bir Be : 


that is, as if the beam were constrained to follow the crystal planes. 
Evidently in practice the multiple reflections will emerge over a range 
not exceeding BC, gradually becoming, with increasing thickness, 
relatively more and more concentrated towards the point M. The total 
‘spread of the emergent beam, counting the direct ray, is therefore not 
less than 

T” sin 0 T sin 6 


BM= Sn ee 
sin 6 sin ¢ sin d—#@ 
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The range of values of ¢ and @ under typical experimental conditions 
is shown in fig. 3. It appears that at least in the case where a typical 
Laue pair of lines, light and dark, is present, ¢ will not be far from 90° ; 
in such a case the separation of the deficiency line from the centre of e 
; concomitant reflections will be of order 7 tan 6 but tending to be smaller 
if there is less secondary . extinction. Strong absorption can shift the 
reflections in either direction according to the sign of 0, if there is any 
obliquity (6490°). In the case of a specimen 1 mm thick, if 6 is of 


Fig. 3 


we 


BRAGG 
ANGLE 


by v Ly 
o / EL a 
INCLINATION OF REFLECTING PLANES @ 


Geometrical conditions limiting production of deficiency and reflection lines 
with typical experimental layout. 


order 10° say, 7’ tan @ is about 4mm. This gives an idea of the conditions 
to be fulfilled in order that the direct beam shall be either resolved or 
observed only in company with its multiple reflections, for each element 
of the source will give a repetition of the fundamental pattern, and the 
overlapping will determine the question of resolution. A reasonable 
source would be about + mm in diameter, and ;;mm would be practicable, 
but smaller sources would generally tend to require long exposures of 
the film. A similar condition is imposed if one attempts to resolve the 
angular structure of the beam. The most favourable condition for this 
is to place the film far from the crystal. Whether the minimum distance 
for resolution is practicable depends on the actual mosaic spread of the 
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crystal in the plane of incidence. In some cases the mosaic spread is of 
the order of 1 minute of arc, so that at a distance of one metre, the 
angular structure will be represented by a width of order 0-3 mm, so that 
resolution will begin at about this distance ; actually the line width is 
often considerably greater than the mosaic spread, so that conditions are 
then more favourable. However, even at 1m the exposure is already 
100 times greater than at 10 cm, and at this smaller distance, resolution 
of angular structure will not occur with a source } mm diameter, especially 
if allowance is made for the further blurring due to the multiple reflections. 
We may make a broad distinction between ‘ resolved ’ and ‘ unresolved ” 
divergent beam lines, and this distinction must also be observed in 
choosing the appropriate mathcmatics. We conclude that ‘it is only 
adequate in completely resolved lines to treat the deficiency of the 
transmitted beam as expressed by the factor in expression (2.1). If 
actually none of the structure is resolved, the deficiency relative to 
background will be, in the symmetric case, 


exp (—pT") [ 4(1—exp [—2G"T7"]) du. 
The contrast of the line against the background due to monochromatic 


radiation transmitted in neighbouring non-reflecting directions is due to 
an intensity ratio 


exp (—pT") | 41 —exp [—2@’T])o(u) du f exp (—wT")o(u) du, 
where w(u) is a weighting function depending on the shape of the source. 


If w(u) varies only a little in the effective range of the weighted function, 
the centre-line intensity contrast may be written 


(1/8) J 3(1—exp [—2G"7"]) du, 
where 6 represents a mean range of the function w(w), or source mean 
diameter, and for a circular source is (2/7) angular diameter. Actual 


photographic contrast will also depend on the characteristic curve of 
film blackening as a function of exposure. 


§3. Tur Linz WiptH 


It is natural to attempt to measure the mosaic spread directly by 
divergent beam photography, and to expect the results to correspond 
to the magnitude of secondary extinction, since the extinction coefficient 
Jz is equal to 1/2y/m for a Gaussian mosaic. 

The experimental basis of such an attempt is to take photographs at 
different distances from the crystal, and to rely on the rate of increase 
of the width of lines with distance. We assume that the finite source 
diameter is an additive contribution given by the intercept of the graph 
on the width axis; the mean angular half-width (9) will be the slope 
of the graph. Referring to fig. 4 we see that the contrast for a line with 
resolved angular structure depends on the power ratio, given in the 
symmetric case by 


2(1—exp [—2G"7")), 
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closely resembling a step-function if G@’7" is reasonably large as we shall 
assume. A good estimate of half-width will therefore be given by 


= QT" 
==) —Omaate i" 
or sane (), COnninwic ee gim ne lee 2) 0(3:1) 


-9 re) +8 ies 
DEPARTURE FROM MEAN REFLECTING ANGLE 


Angular distribution of secondary extinction in the symmetric Laue case. 


Thus in general the apparent line width depends upon the amount of 
secondary extinction instead of being simply proportional to the mosaic 
spread. The contrast o,, in unresolved lines for monochromatic radiation, 
in the range of linear film response will be 


oy Be (S—) = 5. Se eee) 


In some measurements on a pentaerythritol crystal, we measured 9 from 
the slope of the width graph, after first correcting for the inhomogeneity 
of the CuKa« line. Contrast o was measured with a Cambridge micro- 
densitometer and corrected for white background by calculation as 
referred to below, §4. The table of values (1) refers to these corrected 
contrasts on (002) and (110) deficiency lines photographed at a distance 
of 6 cm in the symmetric Laue position. 

These results are consistent with our assumption that the lines on 
which o,, was measured were completely unresolved. Knowing Q from 
independent measurements (Llewellyn, Cox and Goodwin 1937), and T, 
we can deduce 7 from eqn. (3.2) and a knowledge of the function of, 
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graphed in fig. 5. This function has been expressed alternatively as the 
power series already given, or in the asymptotic form 


Cy C2 
Hey= V2 In V(2n)2\( + 555 OS aaa) (3.3) 


where OC, is Euler’s constant (~ 0-577....) and C, ~ 0-25. 
Table 1 


(110) 


2-2 10—* 
2-62 10-3 
0-:172-+0-010 0-082 +0-006 


0-183 0-084 
1-75 10-4=0-69’ 1-4 10-4=0-48’ 


Fig. 5 
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20 


INTEGRATED . REFLEXION BY 
TRANSMISSION 5 DARWIN’S SOLUTION 


FOR A GAUSSIAN MOSAIC SPREAD IN 


THE SYMMETRIC CASE 


ee ah a Pa ioe 7 


The latter form shows that for heavy secondary extinction the line 
contrast of unresolved lines or width of resolved lines increases very 
slowly with Q. 

Since 7 is defined as. the r.m.s. of random small rotations about a 
common axis, it is to be expected that it will behave as a symmetric 


tensor of the second rank and will not as a rule be the same for different 
lines of the same crystal. 
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§4. Errnct of Continuous RADIATION 


If in addition to the characteristic line there is energy amounting to 
K times as much in the white radiation, the contrast, o, of lines on the 
photograph is reduced. If in the neighbourhood of the line a typical 
background measurement is obtained free from any selective effects of 
diffraction, the line contrast is simply reduced in the ratio 


o=o,,|(1+KY), 


where Y is the weighted mean absorption ratio exp (—pZ'") for white 
. radiation divided by the ordinary absorption ratio for the characteristic 
radiation. The weighting function is the spectral intensity multiplied 
by sensitivity of the film to the radiation. Both may be represented, 
the one fairly well and one rather roughly, by powers of p so that the 
averaging can be performed with the aid of incomplete [-functions. 
This was done in our experiments so as to find the approximate value 
of (1+KY); the estimated value was 1-49 for the circumstances of 
table 1, and was used to obtain o,, from measured values of o. 

To estimate K accurately was not practicable, but an attempt was made 
to apply existing data obtained by Ulrey (1918) and Unnewehr (1923) 
on ordinary x-ray tubes of that time, together with the white-radiation 
formulae of Kramers (1923) and Davis (1918). Photographic sensitivity 
was taken to be within limits proportional to exposure (Charlesby 1940) 
and its variation with wavelength according to a curve obtained by 
smoothing out the effects of absorption edges in Seeman’s data (1950). 

A further consequence of the admixture of white radiation is that 
instead of increasing indefinitely with crystal thickness, the contrast 
passes through a maximum (Lonsdale 1947a). Detailed calculation 
showed that the optimum thickness should be slightly different for 
resolved and unresolved lines, but was not critical except as to order of 
magnitude, and did not greatly exceed the simplest estimate (1/1) 
_ obtained from the linear absorption p, of the characteristic radiation in 
the crystal. 


§5. Primary EXtTINcTION 


Since eqn. (3.2) can always be solved for 7, given Q’, it does not suffice 
to determine these quantities uniquely. Only in the case of weak 
secondary extinction does the attempt to measure 7 directly succeed ; 
in this case the width of line is more closely connected with mosaic 
spread than in the case of strong secondary extinction. 

For the present purpose we can express the effect of primary extinction 
by means of a factor which diminishes the value of Q to some value Q’. 
This value is only known as a function of block size for very special cases ; 
for reflection from or through an infinite plate of non-absorbing crystal 
(fig. 6), or by reflection from thick plates of slightly absorbing crystal. 
Neglect of the slight absorption introduces errors of order 18% (Darwin 
loc. cit.). It is not to be expected that the extinction factor even for 
finite rectangular blocks can be obtained by simple interpolation but 
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such a procedure may give approximate indications. We may perhaps. 
regard as the appropriate interpolation ratio s,/s, in the diagram (fig. 7). 
However, for comparison with experiment we have kept the two extreme 
cases apart, and the results of calculation are plotted as values of o,, 
against block size, with 7 assumed to be constant for each line (fig. 8). 


Fig. 6 


xe) 
PRIMARY 


PRIMARY EXTINCTION FACTORS 
Os 

FOR NON-ABSORBING CRYSTALS 
PBRAGG CASE” 


Choice of interpolation factor for rectangular block (§ 5). 


In a series of tests in which a crystal of pentaerythritol was repeatedly 
immersed in liquid air, line contrast was found to have changed at first 
after each immersion, but before complete disintegration the contrast 
had reached a more or less settled value for each line. It is the final 
values that appear in table 1. In view of their agreement with the 
measurements of 3, we may regard the factor (1+KY) as having been 
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determined with reasonable accuracy, and similarly 7 (last line of table 1), 
on the assumption that primary extinction was negligible at this stage. 
It then transpires that the measurements at the earlier stages on the two 
strong lines of pentaerythritol can be explained by using a common 
abscissa scale in fig. 9. It happens that the spacing of the two sets of 
planes is also the same; but the block dimensions could in principle 
have quite a different meaning for the two lines and for the Bragg and 


Fig. 8 
© 20 
5 10 immersions 
(LINE 
‘CONTRAST 9 Simmersions 
Om O.1S 6 
0:10 x 
REFLECTION IN 
Wimmersions e \ start PRIMARY BLOCKS (Az L) 
O . . 
pee =onea | | (902) 
vi 
eS 
C) 
0:05 |- Start ~ LAVE (0°2) 
— BRAGG @ heey: 
— LAUE (I 1 o) 
Lt 
Oo 1 a, ce} 4 Ly & JO” REFLECTING PLANES 


“Changes produced in the contrast of two lines of pentaerythritol on liquid air 
immersion. — Calc. effect of primary extinction. O Experimental points 


plotted on the basis of fig. 9. 


Laue cases. In the Bragg case the relevant dimension is parallel to the 
-erystal face, in the Laue case perpendicular ; only in the latter case is 
a block dimension common to both lines involved. It seems that the 
‘results are at least consistent with the view that during most of the 
-experiment the blocks were of such a shape that the Laue case was the 
more relevant. This is probable for (110) from the tetragonal symmetry, 


-since 6 in fig. 7 is fairly small. 
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The estimated values of 7 for the two lines are, not unexpectedly, 
somewhat different though the larger value was observed for (002) 
contrary to a remark by Llewellyn, Cox and Goodwin (loc. cit.). Had we 
attributed the contrast changes purely to changes in y, we should have 
had to assume an increase from 0:3’ to 0-6’ for (002) and from 0-24’ to 
0-48’ for (110). If, on the other hand, the changes were purely due to 
splitting of blocks, the initial mean block dimension is determined to 
order of magnitude as about 50-60. Although this is an appreciable 
fraction of the crystal thickness it is probable that the extinction formulae 
remain substantially applicable. An attempt to detect such blocks 
microscopically was unsuccessful. 


Fig. 9 


Xlo* PLANES 
3:0 


ASSUMED FOR 
BOTH REFLECTIONS 
IN FIG. & 


(@) VG LS Se SG ey ae So ma OmN LIQUIDEATR. 
IMMERSIONS 


Block dimensions assumed, to fit curves in fig. 9. 


§6. Tor BorRMANN Errect 


It has been observed by Borrmann (1950) and theoretically confirmed 
by von Laue (1949) that an abnormally low absorption can occur in 
perfect crystals in the case of an x-ray beam at the reflecting angle. 
This reduction is likely to be incomplete in small blocks, but we may 
investigate the effect in mosaics by supposing p reduced to some unknown 
smaller value cu, taking this to apply as a mean in the reflection range. 

A transmitted ray spends a fraction €G/Q of its path in blocks oriented 
so as to reflect, where € denotes the diffraction width.; different values of 
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u are distributed almost uniformly within that range so we may conclude 
that the apparent absorption neglecting multiple reflections is 


Msp (1 =) +e “ AeG! 


lc 

=p+ (1- ) G=y+G, say. apes yc eae OnLy 
The effect therefore resembles a small fractional change in G(uw) inde- 
pendent of w and therefore acts in the same sense as primary extinction. 
The transmission is now seen to be in the symmetric case 

On=1l—exp (—G,7") cosh GT”. 
If GT" is large, the condition for line reversal, the ‘ deficiency ’ becoming 
negative, reduces to 
/(2:) In2 7 
l—c_ €’ 
this demands a very much greater thickness than for a perfect crystal, 
for which the analogous inequality is 
In 2 
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having regard to the probably much smaller ¢ in the second case. 

Further calculations have shown that for organic crystals it would be 
difficult to notice the effect on the contrast of lines, though it would 
become more important for weak lines than for strong ones, and more 
important for more strongly absorbing substances. 


§7. FurTHER EXPERIMENTAL POINTS 
Lonsdale (1947 b) in a reference to Ehrenberg mentions that source 
dimensions may have an influence in blurring divergent beam lines. 
Detailed calculations ($3) have shown that this effect may be crucial 
for their visibility. In the first place we have shown that f, the effective 
source dimension, appears in the denominator of the contrast of un- 
resolved lines. It will also determine the distance beyond which different 
parts of a Laue pair will, through inhomogeneity of the characteristic 
line, increasingly overlap and so suffer partial cancellation. Illustrating 
these points, we found first that a source with diameter 0:23mm gave 
better lines than one with 0-3mm. From several crystals of ammonium 
dihydrogen phosphate however, only one specimen gave even faint lines ; 
others gave none even after having been immersed in liquid air. However, 
with a hole of 0-06 mm diameter accurately drilled in a tantalum plate, 
fairly well defined lines were obtained after 1} hours from a crystal 
0:25 mm thick ; with the larger source this specimen gave none. 
Calculation shows that with a given experimental arrangement the 
contrast passes through a maximum as a function of 4 and is small for 
small values of 7. This conclusion is relevant to the interpretation of 
the results of liquid air immersion. This treatment is known to be 
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effective in diminishing the effect of extinction on integrated reflection 
(e.g. Llewellyn, Cox and Goodwin 1937) and has also been found by 
Lonsdale to be often effective in producing divergent beam lines. Success 
is usually attributed to the effect of sudden cooling in producing innumer- 
able cracks and so, by reducing the blocks to a finer texture, diminishing 
extinction and increasing the reflecting power Q’. The process does not 
always succeed, however. Lonsdale has shown (1947 a) that the texture 
of both mosaic and perfect diamonds is apparently immune to any change 
except by much more drastic treatment. Whilst other crystals may also 
have the thermal and mechanical properties to simulate the resistant 
behaviour of diamonds, this cannot be deduced directly from the contrast 
on a single divergent-beam photograph, without careful consideration of 
the resolution. If the mosaic spread is extremely small the unresolved 
lines will probably be faint or absent, and this cannot be taken as a direct 
proof that primary extinction is operating, without detailed calculation. 
So then, if treatment results in no improvement, it may be because the 
texture is already fine enough to prevent primary extinction. Again, 
we may safely assume that the mosaic spread will be either unchanged by 
immersion, or increased, perhaps by the release of internal strains pre- 
viously present in the crystal. The resulting effect on contrast will 
then depend on which domain of angles includes y, and will be favourable 
or unfavourable to line production according as 7 is below or above the 
angular subtense of the source. 


ACKNOWLEDGMENTS 


We wish to thank Professor E. G. Cox for suggesting the enquiry and 
for constant advice and criticism, Professor K. Lonsdale for helpful 
discussion at a later stage, and Professor R. Whiddington and Dr. G. E. 
Brindley for the use of a microdensitometer. 


REFERENCES 


BorrMann, G., 1941, Phys. Z., 42, 157 ; 1950, Z. Phys., 127, 297. 

CHARLESBY, A., 1940, Proc. Phys. Soc., 52, 657. 

Darwin, C. G., 1922, Phil. Mag., 48, 800. 

Davis, B., 1918, Phys. Rev., 11, 433. 

Kramers, H. A., 1923, Phil. Mag., 46, 836. 

Lave, M. von, 1949, Acta Cryst., 2, 106. 

Lirewe yn, F. J., Cox, E. G., and Goopwiy, T. H., 1937, Chem. Soc., 183, 
883. 

Lonspae, K., 1943, Natwre, Lond., 151, 53; 1947 a, Phil. Trans. A, 817, 
240, 219; 1947 b, Min. Mag., 28, 14. 

May, W., 1949, Nature, Lond., 163, 569. 

Seeman, H. E., 1950, Rev. Sci. Instrum., 21, 314. 

Urey, G. T., 1918, Phys. Rev., 11, 401. 

UNNEWERR, E. C., 1923, Phys. Rev., 22, 529. 

Water, I., 1926, Ann. Phys., 79, 261. 


- 


[ 1248 ] 


CXXIV. Programming a Digital Computer to Learn 


By Antuony G. OrrrincER* 
University Mathematical Laboratory, Cambridget 


[Received June 11; revised August 20, 1952] 


ABSTRACT 


Techniques for programming a digital computer to perform simple 
learning tasks are discussed with reference to two concrete examples. 
These examples have a particular bearing on questions proposed by 
Turing (1950) and by Wilkes (1951). By the application of the 
techniques described, digital computers can be made to serve as models 
in the study of the functions and of the structures of animal nervous 
systems. 


§1. INTRODUCTION 


In recent publications, Turing, Wilkes, Shannon, von Neumann and 
others examine in theoretical fashion several facets of the problem 
of the design of machines having properties usually associated only 
with the nervous systems of living organisms.{ Turing (1950), writing 
on ‘Computing Machinery and Intelligence’, examines the question 
“Can Machines Think?’ and, finding this question too vague to be 
answerable as such, proposes that it be replaced by a related but less 
ambiguous question based on an imitation game played with a machine 
and a human interrogator. In addition, he suggests that the term 
‘machine’ be restricted to mean ‘universal digital computer’, since 
such machines have the important property of being able, when provided 
with a suitable programme, to mimic arbitrary machines in a very general 
class. Wilkes (195la) discusses the possibility of giving a digital 
computer the ability to learn. Stating that under special conditions 
* it is easy to see how a machine could be programmed so that it appeared. 
to learn’”’, he then poses the more speculative question of “ whether 
it would in principle be possible to construct a generalized learning 
programme which would enable an operator, if he had sufficient patience, 
to ‘teach’ the machine any subject he chose”. Shannon (1950) 
outlines the requirements for a machine capable of playing chess, on the 
assumption that the rules of the game and the necessary principles 
of strategy are made known to the machine from the start. von Neumann 
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(1951) expresses the feeling that “it is unlikely that we could construct 
automata of a much higher complexity than the ones we now have, 
without possessing a very advanced and subtle theory of automata and 
information ”; he describes some of the difficulties to be faced when 
designing complex automata, and briefly sketches some potential 
components of a general theory. 

In spite of the wealth of conjectures now available in the literature, 
attempts to produce concrete realizations of these conjectures have been 
sparse. Ashby (1948) and Grey Walter (1950) have reported the 
construction of some electro-mechanical devices specially built to 
perform a limited number of the functions of living organisms. The 
theoretical importance of universal digital computers has been discussed 
by the writers quoted earlier. This paper describes some practical 
techniques by which a universal digital computer can be made to exhibit 
a variety of modes of animal-like behaviour. Two concrete examples 
which have a particular bearing on the questions proposed by Turing 
and by Wilkes are described in detail. In each of these examples a 
machine modifies its behaviour on the basis of experience acquired 
in the course of operation. The techniques employed may be of interest 
to all those who are concerned with the non-numerical applications 
of digital computers, and to those who, like psychologists and neuro- 
physiologists, are interested in the potentialities of existing digital 
computers as models of the structure and functions of animal nervous 
systems. 

A brief description of those properties of universal digital computers 
which are important in the present application will be given as a 
preliminary. In that description particular reference will be made to 
the properties of the Edsac,* the machine used for these experiments, 
but this entails no essential loss of generality. Information, coded 
in the form of holes on standard teleprinter tape, is introduced into the 
Edsac from the outside by a tape reader which, by means of light beams 
and photocells, transcribes sequences of holes into trains of electrical 
pulses. The pulse trains formed in this fashion are held in the machine’s 
store, sometimes called its memory. The store is divided into bins 
or storage locations, each capable of holding seventeen pulses. The pulses 
are usually interpreted as binary digits, the presence and absence of a 
pulse in a given position representing 1 and 0 respectively. For easy 
identification the storage locations are numbered serially. Information 
can be extracted from the store and sent to the arithmetic unit, where 
operations are performed on it; the result may then be returned to an 
appropriate storage location. Finally, information held in the machine 
can be sent to the outside world by using the electrical pulses to actuate 
a teleprinter or a tape punch. All these procedures are carried out 
under the direction of a control unit. The solution of a given problem 


* The Edsac (Electronic Delay Storage Automatic Calculator) is in operation 
at the University Mathematical Laboratory, Cambridge. 
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proceeds step by step as a series of elementary operations, each of which 
is specified by an order. When a mathematician assembles the set 
of orders required to work out the solution of a problem he is said to be 
programming this problem for the machine. Eight of the twenty orders 


available on the Edsac are listed in table 1, together with their respective 
functions. | 


Table 1. Selected Edsac Orders 


Add C(n)* to C(Acc)}. 

Subtract C(n) from O(Acc). 

Transfer C(Acc) to storage location n, and clear the accumulator. 

Transfer C(Acc) to storage location n, and do not clear the 

accumulator. 

n If C(Acc)>0, execute next the order which stands in storage 
location n ; otherwise proceed serially. 

n Collate the number in storage location n with the number in the 
multiplier register and add the result into the accumulator : that 
is, add a ‘ 1’ into the accumulator in digital positions where both 
numbers have a‘ | ’, and add a‘ 0’ in other digital positions. 

Pa Read the next row of holes on the input tape, and place the result 

in storage location n. 

On Print the character now set up on the teleprinter, and set up on the 

teleprinter the character represented by the five most significant 

digits in storage location n. 


sse3es 


Q & Gia 


* “ C(x) ’ means ‘ the contents of storage location x ’. 
+ ‘ C(Acc) ’ means‘ the contents of the accumulator ’. The accumulator is a 
special storage location in the arithmetic unit. 


In table 1, as on the programmer’s worksheet, an order is represented 
by a function letter which specifies the function of the order, and by a 
number n, or address, which identifies the storage location holding the 
information to be operated on. From the programmer’s worksheet 
orders are transcribed on to tape in a suitably coded form, and hence 
into the store, where they are held in the same way as other information. 
Normally the control unit selects orders from the store in the proper 
sequence, and supervises the performance of the corresponding 
operations. The orders are transferred from the store to the control 
unit by a special channel. However, when desirable, it is possible 
to treat an order as if it were a number. An example of the numerical 
interpretation of orders is given in table 2, where C (201) represents 


Table 2. An Example of the Transformation of Orders 


Location Content 
100 A 201 
101 A 200 
102 T 201 
200 1 
201 AND 

4P2 
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the order A 25. When C (100) is interpreted as an order by the control 
unit and is obeyed, C (201), interpreted as a number, is added into the 
accumulator, presumed to have been previously cleared. After C (101) ~ 
has been obeyed, O (Acc)=C (201)+C (200)=the number representing 
the order A 26. C (Acc) is finally transferred to 201. Hence when 
C (201) is selected by the control unit as an order, it will be interpreted 
as A 26. Only the address was modified in this example, but the same 
principles apply to the modification of function letters, since these can 
also be interpreted as numbers. 

In normal operation all the information held in the machine, other 
than orders, is interpreted as numerical information, but this need not 
always be the case. The validity of the result of any of the non- 
arithmetic operations C, I, O, ZT’ and U evidently does not depend on 
a numerical interpretation of the operand. With some restrictions 
this is true also of arithmetic operations such as A and S. Hence a 
non-numerical problem can be treated by numerical means whenever 
it is possible to establish a correspondence between propositions of the 
language in which the problem is formulated, and propositions in numerical 
language, in a manner such that the manipulation of numerical 
expressions yields expressions corresponding to those which would result 
if the whole process had been carried out in the original language. For 
example, in the ‘shopping’ programme described further on, a system 
of ‘shops’ is represented by a matrix. In the ‘ response-learning ’ 
programme, stimuli of varying intensity correspond to numbers of 
similarly varying magnitude. 

In normal operation the control unit selects orders serially ; i.e., after 
using the order C (n), it selects C (n+1), and so on. The insertion 
of an H# order into a programme radically alters this procedure. It is 
the function of an EH order to decide whether operation should proceed 
serially, or whether an order held in another storage location, m, should 
be obeyed instead. In the latter case serial operation is resumed 
in a new sequence. The great importance of the # order arises from 
the fact that the choice between these alternatives is based on the results 
of earlier operations. Hence, while the programmer must essentially 
foresee the need for a choice and provide # orders to meet this need, the 
actual decision is made by the machine itself on the basis of information 
obtained in the course of operation either from its own store or, by means 
of the I order, from the outside world. 

Since a computing machine like the Edsac is equipped with an 
information store, it is trivial to programme it to absorb information 
from the outside world and to reproduce this information verbatim 
when required. For example, to make the Edsac first ‘ memorize’ the 
English alphabet and then ‘recite’ it, by printing it on the teleprinter, 
only a series of J and O orders and an input tape holding the alphabet 
in coded form would be required. More interesting types of learning 
behaviour can be obtained by giving the machine a programme which 
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provides for the transformation of its own orders in the fashion described 
above. By the liberal provision of E orders in the programme the 
machine is enabled to organize new information meaningfully and to 
select alternative modes of behaviour on the basis of this organization. 


Two programmes in which these techniques are applied will now be 
discussed. 


§ 2. PRoGRAMMING THE Epsac To Go SHOPPING 


The properties of the shopping programme are summarized in the 
following illustration. The shopping programme may be imagined to 
define a small child sent on a shopping tour. If this child were asked 
by its mother to buy different articles in various shops, it would not 
know at first where to find these articles, and would hunt for them 
by going from shop to shop in a random fashion until it came to the 
desired one. Having found an article once, it would remember in what 
shop, and would go directly to this shop the next time the same article 
were requested. In addition, as its curiosity would prompt it to note 
the whereabouts of articles for which it had not had a specific request, 
when such a request did come, it would often be able to go directly to 
the right shop. 

Before giving the abstract specification of the shopping programme 
it is helpful to digress for a moment to define the boundary between the 
“child machine’, defined by the shopping programme, and the outside 
world. Obviously the simplest definition would be to consider the whole 
of the Edsac as the machine, and to represent the outside world by a 
suitably coded input tape. Unfortunately, tape reading is the slowest 
operation of the Edsac. It is therefore convenient to make a conceptual 
division of the Edsac into two parts, one to play the role of the 
experimental subject and designated s-part or s-machine, the remainder 
to serve as an extension of the outside world. By means of this 
expedient it becomes possible to conduct experiments at the highest 
possible machine speed, since less frequent use of the input mechanism 
need be made. 

Provision is made to enter a description of the shops into the part 
of the Edsac store belonging to the outside world as an m xn matrix. 
This matrix has elements a;,; such that a;,=1 if shop 7 has article 9, 
and a,,=0 otherwise.* The matrix will be called the stock matrix S. 
The ith row of S may be regarded as a row vector whose components 
define the stock of shop 7, and will be called a shop vector. The ith 
shop vector is held in storage location m-+t, where m is a reference 
address. Similarly, the jth column has components representing 
the shops where article 7 may be found, and it will be called an article 
vector. The stock matrix shown in table 3 has eight rows (shops) and 


* A representation of this type was used by N. G. Parke (1949) to specify the 
membership of a set of committees. 
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seven columns (kinds of articles). 


Table 3. A Typical Stock Matrix and the Corresponding Shopping Tours 


Shopping Tours 


i 


(1) 0000010 
0000080 


(2) 0000010 
0000080 


(3) 0000100 
0000480 


(4) 0000100 
0000480 


(5) 0000001 
0000480 


(6) 0100000 
0550480 


(7) 0010000 
‘ 0555480 


(8) 1000000 
1555480 


(9) 0100000 
1555480 


(10) 0010000 
1555480 


(11) 0001000 
1555480 


(12) 0000100 
1555480 


(13) 0000010 
1555480 


ee a ee eee 


= Pseudo-random ’ numbers good enough for these experiments are generated 
by squaring a certain constant and by selecting a number of digits from the 


2 
13465728 


1234 


62413875 


45 


3541 


A shopping order for article j is 
given to the s-machine by a row vector (order vector) consisting of 
zeros everywhere except in column j, where there is a 1. 
of column 1 (table 3) shows an order vector for article no. 6. 
s-machine proceeds to select a shop number 7, at random,* and to form 
the address m+i, of the corresponding shop vector. 
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Stock Matrix 


Shops 


Articles 


> 
1001000 
1010000 
0110000 
0000100 
0111000 
0001000 
0111100 
0000010 


result. The middle digits of the square then serve as a new constant. 


The top line 


A C order 
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(table 1), Cm-+i,, is constructed with this address and, when obeyed, 
collates the shop vector with the given order vector. After the C order 
has been obeyed the content of the accumulator is different from 0 if, 
and only if, shop 7, has article j in stock, that is if a, ;=1. If this is not 


the case, another shop number i, is selected, and so on with Cte ee 
until either the article is found or all the shops have been scanned. 
In the first shopping tour of the example given in table 3 the shops were 
scanned in the order indicated in column 2; column 3 shows that article 
no. 6 was found in shop no. 8. When an article is unavailable, as is the 
case with no. 7 in the fifth shopping tour, the s-machine tries all shops— 
the order of search being indicated in column 2—then prints zero in | 
column 3 to indicate the absence of the article. Once an article has been 
found, the number of the shop in which it was first discovered is recorded 
in the appropriate column of the s-memory, as shown in the case of the 
first shopping tour by the second line of column 1. Before attempting 
a search for an article the s-machine examines its memory to determine 
whether or not the location of this article is already known to it. If 
so, it proceeds directly to the shop where the article was found on the 
earlier trial. The second shopping tour provides an illustration of 
this mode of behaviour. When an article j is found in shop 7 a second 
collation with the shop vector is made, using an inspection vector, 
a row with zeros everywhere except in the two columns adjacent to 
column j. In this fashion the presence of the articles j+1 and j—1 
in shop ¢ can be detected, and recorded in the memory. In the sixth 
tour, for example, article no. 3 was found to be in shop no. 5 while 
shopping for article no. 2. In the seventh tour the request for article 
no. 3 was met without search, as a result of the discovery made during 
the sixth tour. 

If the s-machine memory were cleared and a new arbitrary m’ xn’ 
(m'< m, n'< n) matrix S’ substituted for the matrix of table 3, then, 
proceeding as above, the s-machine would soon discover the location 
of any article requested. Most observers would hesitate to attribute 
this ability to intelligent thought on the part of the machine, and in 
the absence of agreement about the meaning of the word ‘thought ’ 
it is difficult to decide whether or not a machine can be said to think. 
Turing (1950) rightly regards as inappropriate attempts to find the 
meaning of ‘thought’ by an examination of the common usage of the 
word. Instead, he suggests that the question ‘Can Machines Think ?’ 
be replaced by a related but less ambiguous one in terms of what he 
calls the ‘imitation game’. He postulates a game played by a man A, 
a woman B, and an interrogator C. C knows the man and the woman 
only by the labels X and Y, and the object of the game is for C to make 
the correct identification, that is to match X and Y correctly with A 
and B.. A, B, and C are in separate rooms, and communicate only by 
means of teleprinters. If, when A is replaced by a machine, C is wrong 
in his identifications as often as when A was a man, the man and the 
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machine become indistinguishable to C. The question ‘Can Machines 
Think ?’ may therefore be replaced by the question ‘ Are there imaginable 
digital computers which would do well in the imitation game ?’. 

The imitation game can be played with the shopping s-machine taking 
the part of A, if A and B are both given the same stock matrix, and if 
the questions are restricted to shopping orders of the form ‘ In what shop 
may article j be found ?’ coded as order vectors. The answer expected 
from X or Y would be the appropriate shop number, together with a 
list of the shops that were tried, ie., the data given in column 2 and 
column 3 of table 3. Under these conditions the interrogator C would 
find it difficult to make the correct identification. Potentially, of course, 
the man could answer a far wider range of questions than the machine, 
but it is clearly possible, within limits, to find that a man thinking and 
a computing machine obeying a learning programme are behaving 
identically as far as an observer is concerned. Hence, by Turing’s 
criterion, a denial of a machine’s ability to think which hinges on a 
question of the sort ‘ Yes, this machine certainly gives satisfactory 
answers to questions of such and such a category, but what can it do 
if I ask it about something else ?’ appears to be not so much a denial of 
machine intelligence as an admission that this intelligence exists, 
although in a very limited form. 

The shopping programme has some severe limitations. The size 
of the stock matrix has an upper limit and its position in the store is 
fixed, and the method by which the stock matrix and the s-memory 
are scanned is explicitly prescribed. Unless a technical fault develops 
in the machine in the course of operation, the programme does precisely 
what it was designed to do, and information presented to the machine 
is interpreted either as a shopping order, in which case the machine 
proceeds normally, or as nonsense, in which case the machine soon stops. 
There is, thus, only an extremely small chance that the s-machine could 
understand and sensibly answer any request other than a shopping order, 
and an even smaller chance that such a happy accident would be repeated. 
While the shopping machine appears to learn, it certainly cannot be taught - 
any subject the experimenter chooses. 

Several improvements suggest themselves. Some, as the introduction 
of memory decay or forgetfulness, would give the learning process a 
greater degree of verisimilitude, that is a greater resemblance to the 
mental operations to which we are accustomed. Others would increase 
the learning capacity of the s-machine. For example, a more general 
programme could be designed so as to make it possible to teach the 
s-machine a method of scanning the shop matrix and the s-memory. 
In the programme described above this scanning method was explicitly 
prescribed by a fixed sub-programme. The more general programme 
would allow the s-machine to construct a scanning sub-programme, 
by trial and error, subject to indications of progress by an external 
“teacher ’. The method of trial and error itself would now have to be 
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specified, and hence the set of sub-programmes which could be tried 
by the s-machine, as well as the order of trial, would essentially be defined 
implicitly by the initial programme. What would remain for the 
s-machine to do would be to operate until it obtained explicitly the 
particular scanning method which met with the approval of the ‘ teacher ’. 
Once the scanning method had been learned this machine would proceed 
as before.. If the approved sub-programme were not contained in the 
trial set there would be no way of teaching the s-machine the 
corresponding scanning method. An s-machine capable of being taught 
a scanning method would be capable of competing successfully in a less 
restricted form of imitation game, a form in which the simpler s-machine 
would fail; but it is more powerful than the simpler machine only in 
the sense that its limitations are shifted to a higher level of abstraction. 
Making it possible to teach the shopping machine the scanning method 
which it must use in turn to learn the contents of the shops is only one 
of many imaginable generalizations of the shopping programme. 
Additional generalizations, and generalizations of these generalizations, 
could broaden the range of subjects which can be taught to the s-machine, 
but in every case the method of operation of the machine at the highest 
level of generalization would have to be prescribed, and the range of 
operation consequently limited. 


§ 3. A ReEsponse-LEARNING S-MACHINE 


The shopping programme defines an s-machine capable of performing 
functions which, in living organisms, are considered to be the result 
of intelligent behaviour. The response-learning programme, on the 
other hand, defines an s-machine operating at a level roughly corresponding 
to that of conditioned reflexes. 

The response-learning s-machine will be described in two stages, 
with reference to one specific example. The first stage of the description 
is given in § 4, from the point of view of an experimenter who can control 
the input to the machine and observe the output, but for whom the 
machine is essentially a closed ‘black-box’. This point of view is 
analogous to that of an experimenter who is attempting to deduce the 
structure and the internal mode of operation of an organism from 
controlled observations of its functions. In the second stage (§5), 
the factors determining the behaviour of the s-machine are revealed, 
and analysed from the privileged point of view of the designer of the 


learning programme. 


§ 4. THe FUNCTIONS OF THE RESPONSE-LEARNING S-MACHINE 


This s-machine is provided with a sensory device capable of detecting 
a stimulus of varying intensity represented by a number whose magnitude 
corresponds to the intensity. Having defined this representation, it 
simplifies further descriptions to refer to ‘the number corresponding 
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to a stimulus’ simply as ‘the stimulus’. The same applies to the 
notions of response and of approval introduced below. When the 
s-machine begins to operate, a stimulus s, > 0, applied at a time f, 
initiates at random one of a set of possible responses R,, 1=1, 2, . . . , 5. 
The machine signals the occurrence of the response Rk, by printing the 
number 7, after which the operator may express his approval or 
disapproval. In the example of table 4, s, is displayed in column 1, 


Table 4. A Typical Response-Learning Experiment 


1 pe 3 4 ig. 6 ¥ 8 
t St Ri ay Sy Soy Ss Sa Ss 
02 xX 03 03 03 03 03 
02 xX 04 03 03 03 03 
02 2 00 03 O7 03 03 03 
02 xX. 03 04 03 Ol 03. 
5 02 3 00 03 04 O7 03 03. 
02 1 00 06 03 03 02 02 
02 x 03 03 02 02 02 
02 xX 02 03 03 02 02 
02 3 00 03 03 08 02 02 
10 02 eX. 03 04 04 02 02 
02 x 02 02 03 Ol Ol 
02 x 02 02 02 Ol Ol 
03 2 00 02 05 03 Ol Ol 
03 5 00 02 03 03 Ol O5- 
iL5) 03 3 00 02 03 08 Ol 02 
03 1 Ol 06 02 03 Ol 01 
04 3 00 03 02 03 Ol Ol 
03 3 00 03 02 04 Ol 02 
03 5 00 03 02 05 Ol 06 
20 03 3 00 03 02 05 Ol —02 
Ol xX 02 Ol 05 Ol Ol 
03 3 —O0l Ol Ol 05 Ol Ol 
03 XxX 02 Ol 03 Ol Ol 
03 3 —02 02 Ol 05 Ol Ol 
25 03 3 —03 02 Ol 04 Ol Ol 
03 xX Ol Ol Ol Ol OL 
03 1 02 06 Ol Ol Ol Ol 
03 3 —(Q] 04 Ol 05 Ol Ol 
03 1 02 04 Ol Ol OL Ol 
30 02 l Ol 07 Ol 02 Ol Ol 
02 ] 00 O08 Ol Ol 05 Ol 
02 1 00 ll Ol Ol Ol Ol 
Ol 1 00 10 04 Ol Ol Ol 
P Ol i —03 ll Ol 05 Ol OL 
35 Ol 1 —03 09 Ol Ol Ol Ol 
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i 2 3} 4 5 6 al 8 
t oo 7 Ri, ay Sy So, S3z Sa, S54 
Ol 1 —(4 06 Ol Ol OL Ol 
Ol xX 04 Ol OL Ol Ol 
02 “XG 04. Ol Ol Ol 02. 
03 I —(4 04 Ol Ol OL 00 
40 03 4 02 Ol Ol Ol 04 Ol 
03 4 Ol Ol Ol 02 04. OL 
02 EX: Ol Ol Ol 04. OL 
03 4 02 Ol Ol 02 06 OL 
03 4 Ol 00 Ol Ol 09 OL 
45 02 4 Ol Ol Ol 02 11 OL 
Ol 4 —Ol Ol Ol Ol 10 OL 
Ol 4 —Ol 00 Ol Ol 1133 Ol 
Ol 4 00 —Ol Ol Ol 13 Ol 
Ol 4 —02 Ol Ol —03 14 OL 
50 Ol 4 —02 06 Ol Ol 13 Ol 
Ol 4 —(04 Ol Ol Ol ibe 05. 
Ol 4 —04 00 Ol Ol 09 OL 
Ol 4 —(04 Ol Ol Ol 06 OL 
Ol xX Ol Ol Ol 03 Ol 
55 02 Exe Ol Ol —(4. 03 OL 
03 IG 02 Ol Ol 02 Ol 
03 exe 00 Ol Ol 02 Ol 
03 XxX Ol Ol Ol 02 03 
03 xX Ol 02 Ol 02 Ol 
60 03 2 00 Ol 06 Ol 02 Ol 


the resulting R; at time t, Ry, appears in column 2, and the intensity, 
a, of approval or disapproval is given in column 3. At t=3, s,=2 
initiated the response R,, and at t=17, s,=4 initiated R,. An X in 
column 2 at time ¢ indicates that s, was too weak to initiate any response 
whatsoever at that time. In the interval 1<i< 12, s,—2 is frequently 
too weak to elicit any response, and those responses that are made 
occur at random. It is possible for the experimenter to train this 
s-machine to give one particular response only, by expressing his approval 
(a; > 0) when this response occurs, and his unconcern (a,—0) or his 
disapproval (a,< 0) otherwise. Conversely, a response can be 
discouraged by repeated disapproval. Table 4 shows that the approval 
signals aj=2 given to R, at t=27 and t=29, and a,=—1 given at t=30 
were sufficient to train the machine to respond to every stimulus 
with R,, except at t=28, when, at an early stage of training, A, occurred. 
An earlier attempt, made at {=16, to teach the s-machine to make the 
response F, failed for reasons whose explanation will be given later. 
As the training proceeds, errors become less frequent, and the learned 
response may be initiated by a progressively weaker stimulus. The 
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attempt to teach the s-machine to give the response R, begins at t=27 
with the stimulus 3; at £=30 the stimulus 2 was tried and found to be 
sufficient, and at t=33 R, occurred when the smallest stimulus, 1, was 
used. At the same time the need for approval diminishes. At t=27 
and t=29 a,—2, but a,=0 at t=31 and at t=32, and R, nevertheless 
is still initiated by the small stimulus 1 at t=33. From t=34 onwards 
R, is discouraged, until it disappears for the stimulus 1 at t=37. At 
t—38 even the stimulus 2 will not initiate it, but at i=39 the stimulus 3 
brings it back. One last sharp disapproval (—4) finally inhibits it. From 
t=40 onwards the same procedure is repeated with R,. Note how, 
at t=42, a premature attempt to reduce the stimulus from 3 to 2 produced 
no response at all. 

It has already been mentioned that a response may be learned by 
the s-machine if encouraged by the experimenter, but if the experimenter 
is neutral and expresses unconcern (a,=0) for every response, it is 
nevertheless still possible for some particular response to occur more 
and more frequently. Eventually, occurring to the exclusion of all 
others, this response becomes a habit. The high frequency of R, from 
t=15 onwards is due to this effect. To train the machine to give R, 
it was first necessary actively to discourage R;, which showed promise 
of becoming a habit. At ¢=21 the stimulus was reduced to 1 to test 
how strong a habit R, had become by that time. As the stimulus 1 
produced no response, s,=3 was used again at t=22. The reappearance 
of R, then indicated the necessity for disapproval. It is also possible, 
under similar conditions, for the s-machine to decay into a lethargic 
state, making increasingly infrequent responses. In table 5, columns 1, 
2, and 3 for 1 <¢ < 15 are identical with the corresponding columns 
of table 4. However, at t=17 and thereafter s, (table 5) was reduced 
to 1, and the frequency of responses dropped sharply. An interaction 
of the habit-forming effect with that of decay and of the stimulus 
strength takes place throughout the operations of the s-machine. 
Examples of this interaction will be given, but first the response-learning 
programme itself will be described. 


§ 5. Tum Response-LEARNING PROGRAMME 


To control the occurrence of responses, a threshold state number S is 
S;50 is associated with each response R; Columns 4 through 8 of 
table 4 display the threshold state numbers at time ¢t, S,, for every i. 
The set of threshold state numbers is held in the s-machine store in 
consecutive storage locations. When a stimulus is introduced at time t 
by suitable use of I orders these storage locations are scanned serially 
and cyclically, and their contents are compared with a scanning constant 
which is decreased by unity at every cycle. When the first largest 
Si, say Sy, is found the scanning stops. S,, found, the s-machine 
forms S,,+-s, and tests for the condition Sy+s, > T by means of an 
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H-order. T is a fixed number, 7 in the present experiments, called the 
triggering threshold. A response can occur at time ft only if the relation 
Sj+s8, > T is satisfied. When this is not the case X is printed in column 2, 
and the machine prepares to receive the next stimulus, Si14- When 
Sys; > T the scanning proceeds, to find and count those Sj, satisfying 
the relation S,j~=S,,, kj. When no such S,, exists, response &, occurs, 
and the number j is printed in column 2. When two or more responses 
compete with each other, the winning response is selected at random. 

After a response R,; has occurred, the machine signals to the 
experimenter with a buzz of its tape reader, and asks for approval. 
The experimenter then introduces a number a, of appropriate magnitude 
and sign into the s-machine. Given a,, the machine proceeds to form 
Sii:7 from 8S, for all 3. When &; has occurred at time f, 
Si 41 =SytG+1+Ni,—N, ;1—d (Sy, t); the terms of this expression 
will be described in turn. 

By increasing, decreasing, or leaving S,, constant, the addition of the 
factor a, correspondingly modifies the probability that S; 44, > S; 2 
7A~j, and hence the probability that the scanning process will stop at the 
storage location holding S,,,,. Adding unity to the threshold state 
of the response which has been initiated at time t increases the probability 
that this response will occur again at time ¢’>#. This device accounts 
for the habit-forming effect described earlier. It is this effect which, 
together with the chance selection of R,; at t=17, accounts for the 
s-machine’s delayed learning of the response R,. Attempts to teach 
this response were begun at t=16, but were unsuccessful until t=29, 
after R, had been effectively discouraged. 

N, and N,;-, are both pseudo-random numbers. In each interval 
(¢, t-+1) a pseudo-random number N,, —5 < N,< 5, is added to one 
S,, selected at random, and N,_, is subtracted from the S,, (kj, or k=7) 
to which it had been added in the interval (t—1, t). In this fashion 
random fluctuations are superimposed on the average level of the 
threshold states. Because of these fluctuations the machine can make 
mistakes, that is, it can occasionally make a response other than the 
one it has been taught to make. More important, provided that no S; 
is excessively large, each S; has a reasonable probability of being greater 
than the others at some time. This makes possible the teaching of a 
new response, or, when no response is favoured over the others, produces 
an interesting variety of responses. 

The last factor, d(S,,, t), produces a decay trend of all threshold states 
toward 1. d is different from zero only in the intervals between t=5n 
and t=5n+1, where n=0, 1, 2,... In these intervals d=+1 when 
S4>1, d=—1 when S;,<0, and d=0 when S,=—1. The effect of d>0 is 
self-explanatory. The negative decay is provided when Sy<0 for the 
purely practical purpose of preventing the ‘ death ’ of the s-machine. , As 
illustrated in table 5, the decay introduces some lethargy into the behaviour 
of the machine, by causing all S, to drop, hence requiring ever increasing 
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stimuli to produce a response. Were all the S; to become so strongly 
negative that S,,+-s,< 7 for alli, given the most favourable random effect 
and stimulus, no further response could be elicited from the machine. 
The use of negative decay effectively makes 1 the average minimum 
level of the S;. The effect of decay and of the random variations on the 
threshold states can be observed most clearly in table 5, when ¢>15. 


Table 5. The Effect of Decay and Random Fluctuations on the Threshold 
State Numbers 


1 2 3 4 5, 6 7 8 

t Sy Ri, ay Sy So, Ss Sy S54 
02 3.6 03 03 03 03 03 

02 AG 04 03 03 03 03 

02 ye 00 03 O7 03 03 03 

02 BXG 03 04 03 Ol 03 

5 02 3 00 03 04 O7 03 03 
02 1 00 06 03 03 00 02 

02 Xx 03 03 02 02 02 

02 xX 02 03 03 02 02 

02 iS 00 03 03 08 2 02 

10 02 aX 03 04 04 02 02 
02 xX 02 02 03 Ol Ol 

02 ».¢ 02 02 02 Ol Ol 

03 2 00 02 05 03 Ol Ol 

03 5 00 02 03 03 00 05 

ED 03 3 00 02 03 08 Ol 02 
03 1 00 06 02 03 Ol Ol 

Ol ».¢ 02 02 03 Ol Ol 

Ol aX 2 2 03 Ol Ol 

Ol XxX 02 02 03 Ol 02 

20 Ol 5 00 02 02 03 Ol 06 
Ol 5K Ol Ol 02 Ole Se=ae 

Ol XxX Ol Ol 02 Ol 02 

Ol ».¢ 00 OL 02 Ol 02 

Ol nG Ol Ol 00 Ol 02 

25 Ol xX Ol Ol 02 Ol 02 
Ol xX Ol Ol Ol Ol Ol 

Ol ».4 Ol Ol Ol Ol Ol 

Ol 1 00 06 Ol Ol Ol Ol 

Ol xX 02 Ol 05 Ol Ol 


For all those R,; which did not occur at time ¢ 
83, 141=Si tN Ny atd(Sy 0). 
In this expression a, and the habit-forming term 1 do not appear. The 
essential features of the preceding description are thus summarized by 
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the three relations which govern the operation of the response-learning 
s-machine : 


1. To initiate a response at time ¢: 
Si+s,>T for some j ; 

2. When R, has occurred at time t : 

Sit =Seta+ 1+ Ny—N,11—d(Sin t) ; 

3. When &; has not occurred at time t: 

Sern =S t+ Ny, 1 1—d(8y,, t). 

The interaction of habit forming, decay and stimulus strength has 
already been mentioned at the end of § 4, and the promised examples will 
now be given. The three experiments recorded in table 6 were performed 
with identical initial threshold states (S,,=1, all i) and approval (a,=0, 
all R;), but with a variable stimulus. 

For A, s,=1, for B, s,=3, and for C, s,—4, at all times. In A the decay 
effect predominates and responses are rare. In B the increased stimulus 
stength elicits more frequent responses, while in C R, soon becomes a 
habit. A and B were, in fact, continued to t=81, and table 7 shows the 
response frequencies obtained. Hence for a given distribution of the 
random numbers J,, it is possible to vary the stimulus so as to obtain either 
one of two unstable states, decay or habit formation, or the steady state of 
frequent and varied responses. It is clearly in the latter state that the 
s-machine exhibits the most varied behaviour, and is most receptive to 
training. 

The reader may find it interesting now to follow once more the descrip- 
tion of teaching processes given with reference to the first three columns 
of table 4 in § 4, but this time using the threshold states given in columns 
4 through 8, and fig. 1, where some of the data of table 4 have been plotted. 
Graph A of fig. 1 displays S,, S;, and S,,, and the triggering threshold 7’. 
In graph B, s, and a, are plotted against the same time scale as the Sj. 
a,=—0 is represented by a dot on the t-axis and a,40 by a vertical line of 
appropriate magnitude and direction. The number of the response to 
which each approval was given is written directly beneath the t-axis of 
graph B. It must be remembered that, while he is training the s-machine, 
the experimenter does not know the threshold states, and must rely on his 
recollection of his own past actions and of the responses the machine made 
to them. This is the data given in the first three columns of table 4. 
With this limited insight, training the s-machine is an interesting occupa- 
tion for which, however, no less patience is required than for work with 
animal subjects. It wasfortunately possible to mechanize the experiments 
of table 6 and to carry them out at machine speed, since the stimulus and 
approval were fixed. 

The behaviour pattern of the response-learning s-machine is sufficiently 
complex to provide a difficult task for an observer required to discover the 
mechanism by which the behaviour of the s-machine is determined. By 
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Table 6. The Interaction of Habit Forming, Decay and Stimulus Strength 
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examining the data of the first three columns of table 4 such an observer 
could easily find regularities in the response pattern, and might even 
develop empirical rules for predicting responses with tolerable accuracy. 


Table 7 


A B 
freq freq 


74 57 
3 
1 
2 
0 
1 


Total trials: 81 
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He would find it very difficult by this means to obtain a good approxima- 
tion to the description of the response-learning programme given above. 
Switching the machine off to dissect it would be of limited value only, 
since this action would make the programme vanish from the store. 
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The response-learning s-machine could, as the shopping machine does, 
hold its own in a restricted version of the imitation game, and, like the 
shopping machine, it suggests generalizations. It is interesting to envisage 
the use of the response-learning programme as an elementary building 
block of a wide variety of more sophisticated systems, including, for 
example, the shopping programme. In the response-learning s-machine as 
described the stimulus and the response are external. However, it is 
possible to form a set of s-machines in which the responses of some serve as 
stimuli for others. The threshold 7’ of any unit could be allowed to vary 
under control of some response or conjunction of responses, thus varying 
the sensitivity to stimulation of this unit. Or again, a response or con- 
junction of responses could provide approval or disapproval for a member 
of the system. Such a chain of controls might be open, or it might be 
closed and regeneration be introduced into the system. No system of this 
type was actually put into operation, but not, as might have been expected, 
because of storage capacity limitations. Before this study was begun it 
was believed that information storage facilities would be the major 
limiting factor in the design of learning programmes of increasing 
complexity, but in fact neither the shopping programme nor the response- 
learning programme require the full Edsac store, and it would be possible 
to programme the Edsac to act as a system of two, and perhaps of three, 
response-learning s-machines. However, it was felt that a generalization 
of the response-learning progremme in this way would illustrate no new 
principles and would be subject to all the limitations already described in 
connection with the shopping programme. Once a learning programme 
has been designed, and the s-machine is operating properly, this s-machine 
is of little further interest to the experimenter unless the functions it 
performs are of interest in their own right, or unless it can serve as a 
model of value in the investigation of parallel problems of animal 
behaviour. 


§ 6. Untversat DicrraL Computers as MopELs 


It is important to distinguish carefully between machines intended to be 
models of the physical structure of animal nervous systems and machines 
designed to perform specific functions of animal nervous systems. The 
fact that the s-machines which have been described can play restricted 
versions of the imitation game implies a correspondence between their 
functions and some functions of the human nervous system. But as there 
are usually many means of synthesizing a given pattern of functions it is 
dangerous to reason from functional to structural correspondence. A special 
machine could be built, for example, to do only the job of the shopping 
s-machine. The physical structure and the internal mode of operation of 
this machine would need to bear no relation to the Edsac and, in turn, to 
the human brain, other than that of superficial correspondence of certain 


functions of input and output. 
4Q2 
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The readiness with which the Edsac, and other universal digital 
computers, can play different roles at a moment’s notice is one of their 
important properties. When the role a universal digital computer is 
called upon to play is markedly different from those its designers had in 
mind, a necessarily large proportion of the orders in a programme must be 
devoted to specifying this role. In the response-learning programme an 
aggregate of elementary Edsac orders is required, for example, to imitate 
the firing of a response, since this is an operation which does not correspond 
to any single order. This drawback is balanced by some important 
advantages. Given the Edsac, the only additional equipment required to 
turn it into a learning machine is a length of teleprinter tape, and changes 
in s-machine design and the rectification of errors require at most the 
punching of a new tape. Only a few seconds of input time are required to 
convert the Edsac from an orthodox computing machine into an experi- 
mental learning device, and deconversion is equally rapid. For these 
reasons a flexible universal digital computer can serve with advantage as a 
proving ground for a wide variety of models. 
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CXXV. On the Destruction of Superconductivity by Large Currents 


By C. G. Kurer* 
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ABSTRACT 


A theoretical study is made of the restoration of resistance in super- — 
conducting wires by large currents. The experiments of Shubnikov and 
Alexeevsky and of Scott show a marked deviation from the theory of | 
F. London. The present theory uses the same geometrical model of 
the structure as London, i.e. a string of superconducting conical domains 
along the axis of the wire. But it takes account of the effect on the 
resistance of the scattering of electrons at successive normal-super- — 
conducting interfaces, when their separation is comparable to the mean 
free path of the electrons in the normal phase. The scale of the structure 
is determined by minimizing the ‘ Magnetic Gibbs’ Function’. It is 
shown that the interfacial surface energy may be neglected. The theory 
is in fair agreement with the rather scanty experimental data. The only 
parameter of the theory is the mean free path of electrons, and the value 
assigned is consistent with that obtained from other phenomena. 


§1. IyTRoDUCTION 


In the macroscopic theory of the intermediate state of superconductors 
(London, e.g. 1948) it is assumed that the intermediate state is homo- 
geneous, with |H|—#,, while | B| has some value between 0 and H.,. 
Various bulk properties can then be predicted. However, when 
observations of such properties are made on small specimens (the meaning — 
of ‘small’ in the present context will be discussed in § 4), deviations 
from London’s theory are observed. In the case of a specimen under 
the influence of an external magnetic field, these small deviations are 
fairly well understood and arise from the surface energy, which will be 
written as 4H,°/S7 per unit area, associated with normal-superconducting 
interfaces (see e.g. Landau 1937 and 1943, Andrew 1948, Kuper 1951, ’ 
Lifshitz and Sharvin 1951). But when the intermediate state is set up 
by currents, the deviations are so large that this explanation is inadequate. 
After a brief discussion of the discrepancies, we will develop a model 
which will account for the experimental facts. 

We will use the following notation: H=magnetic field, B—magnetic 
induction, H,=critical field, J—=electric current, I,—4H,p)=critical 
current (see next paragraph), R=resistance of wire per unit length, 
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&,,=normal resistance per unit length, po=radius of wire, p=radius of 
‘core’ in the intermediate state (see fig. 1). Instead of taking explicit 
account of all the supercurrents, we will use the picture of a superconductor 
in which the permeability »=0. We will assume that p= py when J=I,,. 
In London’s theory of the destruction of superconductivity by currents 
in cylindrical. wires (London 1948, p. 120) the model assumed is as shown 
in fig. 1. The field in the core is everywhere critical (|H|=Z,), and 
outside the core, |H|>H,. Thus the outer sheath is normal, and the 
core is in the intermediate state. Silsbee’s (1918) suggestion that 
I,= 3H ,py (introduced to explain the original observation of Kamerlingh 
Onnes (1911) that superconductivity is destroyed by large currents) is 
confirmed, and it is predicted that the resistance changes discontinuously, 
at Silsbee’s critical current, from zero to half the normal value. The 
relation London obtains between resistance and current is (cf. fig. 2) 


Peer tema (lei yieen(i 1.) oe ee 2. (1) 


The structure of the wire in the intermediate state. The coordinate systems 
to be used are indicated. 


(A result equivalent to this one was given in a private communication by 
Langevin to Silsbee, but London’s was the first published derivation.) 
The angle % (see fig. 1) between a domain wall and a cross-section of the 
wire is undetermined, but provided it is small it will not affect the 
resistance curve. 

In marked contrast, experiments (Shubnikov and Alexeevsky 1936, 
Alexeevsky 1938, Scott 1948) while further confirming Silsbee’s hypo- 
thesis, show a larger discontinuity than is predicted by London. Thus 
when J increases through the value [,, R rises discontinuously from zero 
to about 0-8R,. At J=TI,, the right-hand derivative of R with respect to 
I is finite, while according to eqn. (1) it should be infinite. Also for 
decreasing currents, some hysteresis is observed (see fig. 2). 

Three effects capable of influencing the configuration are left out of 
consideration in London’s treatment. These are: (i) the geometrical 
spreading of the lines of current flow, (ii) thé fact that the current is 
flowing in regions small compared with the mean free path of the electrons 
in the normal metal, and (iii) the surface energy associated with the 
interfaces between normal (N-) and superconducting (8S) phases. 
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(ii) will be referred to as the ‘size effect’. Scott tentatively suggests 
that the deviations from London’s theory may be due to effect (iii), but 
it will be shown to be unimportant compared with spreading and the 
size effect. 
§ 2. OUTLINE OF THE THEORY 

We will assume the same model as London (fig. 1), but will attempt to 
take account of the increase of resistance resulting from the size effect (i). 
It then appears that the discontinuity in the resistance is not independent 
of the angle %; the value that this angle takes will be discussed in 
§§4 and 5. The value of % should be determined by the condition 


Fig. 2 


1-0 


o2 


O06 O08 1-0 12 1-4 -6 ; 
Tae 8 2-0 


The resistance-current curves. The broken curve is that given by London, 
the full curve that predicted by the present theory, using the value 
1=2-8x10-% cm. The points are those found experimentally by Scott 
for his specimen 3, of radius 17-°83x10-' cm. The dotted curve is the 
extrapolation of the theoretical curve into the hysteresis region. 

that the ‘Magnetic Gibbs’ Function’ G=/{(U—TS—B .H/4z) dr is a 

minimum (§4). It is difficult to calculate the effect of the spreading (ii) 

of the current stream lines on G when % is large, but an approximate 

calculation valid for small % will be performed. This will lead to %-values 
around 7/6, unfortunately not small, which will be seen to vary only 
slowly with the radius of the core. (The error in @ results from neglect 
of terms in #°, and is likely to be about 10 to 20%.) To a first approxi- 
mation, we shall thus take y=7/6, independent of p. It will be shown 

in § 5, that effect (iii) leads to only a small change in G. 

An exact calculation of the increase of resistance caused by the size 
effect is also impracticable. But by the use of an approximate inter- 
polation formula (3), we are able to derive an expression for the resistance. 
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The shape of the observed resistance curves is well reproduced (see fig. 2), 
and the variation of R(I,)/R,, with the wire radius Po is in fair agreement 
with Scott’s data (fig. 3). The hysteresis shown by some of Scott’s 
specimens will be discussed in § 5. 


§3. THE Size anp SPREADING Errects 
(a) The Size Effect 


It will be assumed that electrons are diffusely scattered when they 
impinge on an N-S interface. This assumption receives some support 
from the experiments of Pippard (1950), which suggest a ‘ range of order ’ 
of electrons in the S-phase, of the order of 10-4 em, which is small compared 
_ with the mean free path in the N-phase. The existence of this order 
suggests that an electron impinging on the interface, from the N-phase, 
will be scattered in a distance comparable with the range of order. The 
analysis in the present paper would not be modified in substance by 
assuming instead that a large fraction p was scattered, and the remainder 
transmitted. This would merely introduce p as another parameter, and 
in view of the paucity of experimental data, would add little to the 
physical content of the theory. 

We will take a system of cylindrical polar coordinates (7, ¢, z), with 
z as the axis of the wire, and the origin at the point of contact of two 
S-domains, or where it is more convenient for certain approximations, 
spherical polar coordinates (rp, 9, 4), (see fig. 1). In conformity with the 
usual notation, we shall use n for the number of conduction electrons 
per unit volume, e and m for the electronic charge and mass, v for the 
velocity of an electron at the Fermi surface, and / for the mean free path 
of an electron in the N-phase. 

Owing to the geometry of the configuration, it appears that an exact 
calculation by Fuchs’ (1938) method of the ‘ size effect ’ contribution to 
the resistance isimpracticable. The following rather crude approximation 
will be used instead (Thomson 1901, Lovell 1936). The conductivity o 
is assumed to be a point function. Let o be written o,A/l, where 
op=nel/mv is the bulk conductivity. Then A at a point P is to be 
evaluated as follows: Construct a sphere of radius / about the point P. 
Draw a cone of small solid angle dw, with P as vertex. Let A be the 
distance from P to the surface of the sphere or the bounding surface of 
the domain, whichever is the less. Then 


tN aol ans ian tuto. (2) 


Even this approximation is intractable except in two special cases. One 
occurs when P is the origin of coordinates. Here A=/siny%. The 
other is where the sphere does not intersect the domain boundaries, i.e. 
far from the axis of the wire. Here A=/. A convenient interpolation 


formula for intermediate points is 


o/oy=A/l=1—(1—s) exp (—rs/]), - 6 + 6s (3) 
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where s=siny. This yields A/l=sin ¢ for r=0, and has risen to near 
unity when the path length 2rs>/. Equation (3) assumes o to be constant 
along a line of constant r. This is not true, for when P is near a domain ~* 
boundary c/o, is appreciably less than when z=0. We can interpret 
eqn. (3) as giving an average of o along a stream line. In view of this, 
we are able to fix the coefficient of —rs/l in the exponent as unity, since by 
eqn. (2) when rs ~ 1, o(z=0) ~ 9 but o(z=1) ~ 405, so that o/o, should be, 
say, 0-8. Taking s ~ 4, eqn. (3) gives o/oy ~~ 0-75. The detailed choice of 
the form of the interpolation formula has little influence on the final result. 
(6) Spreading Effects 

We shall make two assumptions, true only in the limiting case of ¢ 
small. These are that o is a function of rp only and not of @, and that A 
is constant along a stream line. With the latter assumption the electro- 
static potential V satisfies Laplace’s equation. For* if j is the current 
density, the equation of continuity gives 


a a aes £ j=0. 
Defining a stream function Y, such that 
1),=0P 02 ; rj,—=—9OWP/odr, 
OVS i) 0 tee eee 
9a ~|ijor fil. ee” 
where ¢/da means differentiation along a stream line. Also by Ohm’s 
Law, 


we have 


jr=8, 5 j=08,. 
Hence OV /oz=orL, ; 0 /or=—orEL,,. 


As we have assumed o=ne"A/mv is constant along a stream line, we may 
define y=Jo d¥ and then 


ox OV ox OV 
dz or? ro gion 
whence V?V=0, as stated above. 

If we assume that when the current is greater than I, (and the core 
radius correspondingly less than p,), the current in the core still remains 
strictly confined to the radius of the core, then we may apply the boundary 
condition dV/dr=0 to the surface of the core. The other boundary 
conditions may be chosen as V=0 on the normal cross-section through 
the vertex of the cone (fig. 1); Ve on the conical domain surface. 
In order to solve the Laplace equation we will simplify the conditions by 
replacing the cylindrical boundary r=p by the spherical one rp=p. 
Then the solution is such that the lines of current flow are ares of circles, 
whose centre is the origin. Then by the axial symmetry, the magnetic 
field lines are circles about the axis, and the field at the point (rp, 8, d) is 


9 
a 


'p i 2 Pp cos ws 
H= -| 3) fa i di — —aal ) ; 4 — 
TJ 7p j(1p)drp rp sin BJ 7p j(rp)drp 


sin§@ Hy, . (4) 


*T am indebted to Professor D. R. Hartree for this proof. 
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where H, is the field for 6=7/2—y. Thus the only effect of spreading, 
in this approximation, is to introduce the factor cos %/sné@in H. Thus 
the spread of current does not alter the resistance, but it will affect the 
magnetic Gibbs’ function considerably (§ 4). In a better approximation, 
we should expect spreading to produce a slight modification of the 
resistance as well. 


§4. Tor Macneric Grsss’ Funcrion 


Owing to the presence of irreversible effects (Joule heating and heat 
conduction) it is not immediately obvious that the usual thermodynamic 
methods are applicable, and that the equilibrium configuration may be 
got from the condition that the magnetic Gibbs’ function 

G=f (U—TS—BH/4r) dr 

is minimal. While it should be possible to justify it by the methods of 
irreversible thermodynamics, we may do so more simply by appealing 
to experiment. For Scott (1948) has shown that the form of the resis- 
tance curve (fig. 2) is nearly independent of the temperature, though the 
actual value of the critical current varies through a wide range. Since 
the Joule heating and thermal conduction are both proportional to j?, this 
invariance of the curve of fig. 2 with the value of J, shows that conditions 
in the wires and their heat-baths are nearly isothermal. By performing 
the experiments close enough to the normal transition temperature 75, 
we can make the irreversible effects as small as we choose, and hence we 
can apply the usual thermodynamic formulae. The invariance of the 
shape of the curve then confirms the continued validity of this procedure, 
even at temperatures well below 7’). 

The only terms in G which depend on the configurational parameter 
x are (where 4//,?/87 is the surface energy) 

1 1 


ae { Fe ee | AH? dS, 
87 N-regions 87 J N-S interfaces 
of which the surface energy term will be shown to be unimportant. To 
evaluate 7 y= — | FH? dr/87, we require an expression for H. The 
N-regions 
current density is 
WESC Ae) se ee: Ree ee a9) 
2 ecg 
and thus Hi { Oy, Age ae etre Meee ee re) oh (0) 
0 


where H, is the value of H at 0=7/2—+ (see eqn. (4)), Le., 


Anage {' G , SHE (7 fan, is 
eee ae es [ (a) exp (—rs/0)] dr 
ATraye is i. “yi | 
== ib reat exp ( sx) } Z 


where z=r/1; s=siny. Hence 


' ie 
ete ; AS | 1— =" f—exp (—s1)} | 
sind ob SX 
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Now the potential difference ¢ across half a domain is related to the mean 
electric field by the equation H=e/p and therefore 
1—8 
es aa dao, tie E —{1—exp (—sz)} |. int ete 


At large distances from the axis of the wire, we assume that H tested 2 80 
that we find for p the equation 


A =4trey E prints de. ae) ee 


Thus 
—8r Gy=2r { 


or 
— G y= }p*(1—s*) tog, /(). (4770 Hp)? 


ft tag 2 
x | 2 Ee (1— o*) | Bx*dx per half-domain, 


ar Cone, 
Haridpant ieee I, 


740 | Hyer dr 


N-domains 


or from eqn. (8), if y=p/l, 
1 y l—s z - 
— G y=}? p3(1—s?) 108, | (= )e a J, 1— Spree “| Hols 


per half-domain 


1 1 JI— l—s 1+3s] 
=}H2o s1—siystert tog, /(j 7) j-> 1--- =| 


8 sy sy 28"y 
l=s)ael- 1—s)? ’ 
aaa ; 2 . 84 — ( Dye” ow} per unit length 


darn taneriat ByeRt!) 
(since the length of a half-domain is p tan b= ps//(1—s*)). 


Table 1. —4Gy/H2p? 


1 0-021 0-043 | 0-062 0-0808 | 0-0894 | 0-0880 | 0-071 


5 0-061 0-100 | 0-1276 | 0-1410 | 0-1481 | 0-119 


In table 1, the values of —4% y/H,2 p? are tabulated for varying s, 
fory=landy=5. The maximum is seen to occur at s=0-64 for y=1, and 
at s=0-57 fory=5. Thus the approximation of small ys cannot be valid, 
and little reliance should be placed on the exact numerical equilibrium 
s-values. On the other hand since the error in Gy due to the approxi- 
mation in (4) is O(/°), we should not expect the values to be in error by more 
than about 20%, and we can have rather more confidence in the qualitative 
assertion that s,,,,;, varies very slowly with p. In comparing the theory 
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with experiment, we shall thus (as an adequate approximation) take 
Sequii=%, independent of the radius p=yl of the core. We shall see, 
in § 5, that this value of s leads to a perfectly sharp transition for wires of 
vanishingly small radius, as we should expect. 

We have not yet taken account of the surface energy, 4H,?/87 per unit 
area, associated with the N-S interfaces. The contribution to the free 
energy is 
(Alls es 


Go | 
4 8a N-S interfaces 


ds =} AH 2pL sin 


per half-domain, where L=ptan=the half-domain length; or per 
unit length, 


Grup tH «Bye wee ess (10) 


H,’p/8s is comparable in order of magnitude to Gy, but the factor A/p is 
of order 10-3 (see, e.g. Kuper 1951, for a discussion of the value of J). 
Thus the contribution of the surface energy to the Gibbs’ function is not 
significant ; the minimum of Y will not be shifted far in consequence. 

The meaning of a ‘small’ specimen in §1 is now clear—in those 
intermediate state phenomena where surface energy is important ‘ small’ 
means of the order of magnitude of 4, while where mean free path effects 
dominate, as in the present situation, ‘ small’ means of the order of /— 
about 100 times larger than 4. 


§ 5. THz RESISTANCE—CURRENT CURVE 


The current is made up of two parts, that in the core, J,, and that out- 
side the core, J,. The expressions for J, and I, are 


I,=pH(p)/2; I,=709H (po’—p?) « eee ae Gl tL) 
From eqns. (7) and (8), the sum J=I,+J, may be written 


ve 
I=I1,4+1,=70,H p+? [12 — ae) | /167ro,L, 


l—s: 
or with the MOUAULON foe oe dae ted Se See a PS) 


I=no, Ep? +H 2f?/16709H. ee ea Pee (8) 


(Compare London 1948, p. 122, eqns. (3) to (5).) Following London, we 
can solve (13) to give the resistance, R per unit length, 


R=B[I= 14+ Vl—Lep PBR Lt VOLPE) (Ed) 
ity ee 4) 


We can now account for the observed features of the resistance curves 
which are not reproduced by London’s theory. For while the current I 
is increased from zero to J,, the resistance remains zero, but as soon as it 
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passes I,, the resistance rises abruptly, on to the curve of (14), approxi- 
mately of the form of London’s curve, but with a ‘ critical current ’ less 
than the true critical current. This will therefore account both for the 
restored resistance in excess of 1R,,, and for the finiteness of the right-hand 
derivative of the resistance at J=J,. With increasing current, the core 
will contract, and f? will consequently approach zero. This will result in 
a rather more rapid approach of the resistance to R,, than would be the 
case if f? were constant, but it will not affect the general features of the 
curves. To derive an expression for f? in terms of J, we insert in (12) for 
y=pl/l its value from (8), y=H,/470,HI=H,/4n0,RII=(1,/1)(R,,/R)yYo 
where %=py//. This leads to 


2 —s)IR 
r=ye,(i+/[1—% 1-2 a (1—exp—s1.R,/1R))} }) : 


To solve this implicit equation for R, we approximate, inserting R=R, 
on the r.h.s. to give (putting s=}; see below), 


Sel! bela a els — ae Yole/2I i 15 
be +h, 1+ 1 iE 1 ig S } oe eae koe ( ) 


We can see that using thevalues,,,,,,;=+} leads, for very small wires, to the 
result that €+4, or f?+0. Thus, in the limit as p)—0, we find the whole 
resistance restored at J=I,. This is in accord with an extrapolation from 
the existing experimental data. The retention of s,,,,,—=} for such smal 
wires is partially justified by the following considerations: Although 
eqn. (9) would give 8,9,;>1 as p+0 this is just the condition for the 
approximations made in deriving (9) to have no validity at all. On the 
other hand in the region y~1 to y~10, 8,.,;; a8 derived from (9), si) for 
brevity, is not very different from 4, and the errors in Yy are of order ¥°, 
or about 20%, say. Thus throughout the region of p-values where the 
approximations made in deriving ¥y are at all tolerable, s@~4. Further, 
though these approximations break down in the limit p—0 there is no 
apparent reason why the model itself should do so. Finally, the result 
that the N-S transition should be quite sharp for small enough specimens 
is to be expected, by analogy with a similar extrapolation in other inter- 
mediate state phenomena (e.g. Andrew 1948). 

The form of (14) or (15) suggests the occurrence of hysteresis, for when J 
decreases through the value J,, the curve can formally be extrapolated 
to where it has a vertical tangent. In the present model, it is true, this 
extrapolation has no physical meaning, for it would imply p> pp, which is 
absurd. But if the S-domains have slightly rounded edges at the boundary 
of the core, and consequently do not quite touch the surface of the wire 
when JJ, (this is a small modification of the model, to be expected when 
the surface energy is taken into account) then the extrapolation, for some 
distance at least, will be physically significant. (The greatest possible 
range of currents between which hysteresis can occur is roughly from 
I=fI, to I=I..) 
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§ 6. COMPARISON WITH EXPERIMENT 


The only experiments in which the variation of the form of the resistance 
curves with the radius of the wire has been investigated are those of Scott 
(1948). Three indium wires of differing radius were investigated. The 
mean free path of electrons in indium is unfortunately not known with any 
precision, but has to be determined from the observations as a parameter. 
The continuous curves in fig. 3 are the graphs of R(/,)/R,, and f?(J,) as 
functions of y) (from eqn. (12)). The first three rows of table 2 indicate 
Scott’s experimental values of R(I,)/R,, and the values of f?(I,) deduced 
from them, for varying py. The value of / for best fit is =2-4x 10-3 cm. 
If we express the values of p, in terms of y, using this / then the experi- 
mental points lie as indicated in relation to the curves in fig. 3. 


Fig. 3 


eS) 


2 4 6 8 10 20 40 80 


rae 


Graphs of R(J,)/R,, and f?(I,) against yy. The experimental points of Scott are 
shown, taking J=2-4x10-' cm. After the break, a logarithmic scale is 
used. + R(I,)/R,; Of?(Z,). The rising curve is of f?(J,). 


Table 2 


Specimen Substance Radius (pp) R(L,)/Ry aces 


Scott, 1, la In 53.10" cm 0-854 0-498 
Scott, 2 In 14:3 x 10-3 cm 0:795 0-651 
Scott, 3 In 17-8 x 10-3 cm 0-778 0-689 


Alexeevsky I, I Sn 5-6 x 10-3 cm 0-80 0-64 


SS 
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In fig. 2, we have drawn the curve of eqn. (15) for the value 
J—2-8 x 10-8 em, and indicate Scott’s experimental points for his specimen 
3, of radius pyp=17:8X 10-8 cm. (This value of / is chosen as the best for 
this particular specimen, rather than the mean value, 2-4 x 10~% cm, for all 
the specimens.) 

Although Alexeevsky has not studied the variation of the resistance 
curves with p, in detail, but has only given adequate data for one specimen 
radius, he has fortunately worked on tin, where a good estimate of / 
may be found from the experiments of Chambers (1950) on the anomalous 
skin effect. For this specimen (see the last row of table 2), f*=0-64 
corresponding to yy=5, or since pp=5-6x 10-3 cm, to J=1-1x 10-% cm, 
This is to be compared with the value of / found from Chambers’ result that 
o,/l=constant=9-5< 10-1 mho/em? for Sn. For this specimen we 
estimate o,-1=2:4 to 3-6 x 10-8 ohm—cm, whence /=0-35 to 0-6 x 10-3 cm. 
The agreement is satisfactory, in view of the approximations we have made. 


§ 7. DiscussION 


There are several difficulties in the present treatment. Thus, as in 
London’s theory, we are still left with a singularity of the magnetic field on 
the axis of the wire (as a consequence of the infinite current density at the 
point of contact of two S-domains). But when the surface energy is 
taken into account, we should expect the domain boundaries to become 
rounded, so that instead of touching, they have a separation of order 4. 
The phenomenon of hysteresis observed by Scott has already been inter- 
preted in terms of the rounding of the domain boundaries at the core 
surface. 

The effect of spreading is interesting. Even in the absence of size 
effects, spreading will reduce the magnetic field inside the N-domains, 
though not on the domain boundaries. If the field is H, on the boundary, 
it is H, cos w/sin 6, less than H,, inside the N-metal, which is therefore not 
in true thermodynamic equilibrium. A similar difficulty occurs in 
Landau’s (1937) theory of the intermediate state under external fields, and 
has not really been satisfactorily resolved. Either the intermediate 
state is not in true thermodynamic equilibrium, as seems likely, (though the 
departure from equilibrium is small) or the model is inadequate to deal 
with this difficulty. The ratio H/H, inside the S-metal is, in any case, 
never less than cos 7, or about 0-9. In fact the approximations used in 
§ 3 (b) overestimate the effectiveness of spreading, and H/H, is really still 
closer to unity. 

The ‘ tails ’ in the resistance curves of both Alexeevsky and Scott are 
not accounted for by the present theory. The fact that not all of Scott’s 
specimens showed this behaviour suggests that this behaviour represents a 
departure from ideal conditions—perhaps a consequence of strains or 
impurities. The smallness of the electronic mean free path (l<10-3 em) 
in Alexeevsky’s specimens confirms the impression that the specimens 
were strained or contained impurities, since in very carefully prepared tin 
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specimens J can become about 10-2 cm. Experiments on such specimens 
(and also on specimens in which J had been deliberately reduced by addition 
of impurities) would be of great interest. 
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CXXVI. On the Born—Green Theory of Binary Mixtures 


By G. 8. RUSHBROOKE 
King’s College, Newcastle, University of Durham* 


[Received September 10, 1952] 


SUMMARY 


[he Born—Green equations for binary mixtures are related to certain 
approximations to the integrals occurring in Mayer—Montroll expansions 
of radial distribution functions in powers of density and composition. 
The significance of the vanishing of the denominator in the resulting 
formulae is discussed. 


$1. INTRODUCTION 


THE statistical theory of classical fluids, formulated in terms of the radial 
distribution function and based on the superposition approximation of 
Kirkwood and Boggs (1942), was recently discussed by Rushbrooke and 
Scoins (1951)+ who showed, in particular, that the resulting equations 
for the three distribution functions, g 4 ,(7), 9, p(7) and gpp(7), appropriate 
to a binary mixture (V, A-systems, NV, B-systems in a volume V at 
temperature 7', pressure p, etc.) are mutually consistent only if these 
equations are first linearized in the way that Born and Green (1949) 
and Rodriguez (1949) have linearized the corresponding equation for a 
pure, one component, fluid. This linearization involves writing 


gis(r)=exp [—$4(r)/kT] exp [fi(7)]= exp [—F(r)/AT]A+fi(r)),  ~ (1) 
i=A, B; j=A, B; where ¢ is the interatomic potential function and 
f(r)>0 as r>oo. With this normalization g, ,(r)=g,,(r). And the two 
forms of (1) are effectively equivalent since only terms linear in f,,(r) 
are retained in the linear equations concerned. 

These equations, satisfied by /;,(7), are (eqn. (34) of I) 


rf (7) =2 7p 4 I. (9 (t-+r)[fialE-+r)+€:a%ia(t-+7)] dt €; sa, 4(s) ds 


+2npy] [Hen Lfinlt+r)+ecnean(t+0)] dl espsayn(6) ds ~ (2) 


0 
i=A, B; j=A, B; where 
a(r)=exp[—4(r)/kT]—1, py=Na/V, pp=Nz/V 
and the e’s are parameters which will be discussed more fully below. 


In I we did not give explicitly the solutions of these equations, which 
are easily found by the method of Fourier transforms, for their 


* Communicated by the Author. 
t Referred to hereafter as I. 
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interpretation was felt to be obscure. They have, however, been 
published independently by Fournet (1951). In our notation, these 
solutions are 


__ 1 f?  pya(s)+-ppe?(s)—p 4 pn0(s)[a(s)b(s)—e2(s)] _ sin rs 
Paalr) Gay). (I=p,a(s))(I—ppb(s))—pappelsy) °° 7% 
an? 4 age) 


z | * psals)c(s) + px0(s)c(s)—p pre(s)[a(s)b(s)—c(s)] _ sin rs d 
V(2n)J—o (I= psals))(—pp(6)) paps) 
and similarly for fgp(r). Here py=(27)32p,, pp=(27)2?pp_ and 

a(s)=€4 Py a(S), 0(8)=epnbpn(s), C(s)\=eanPan(s) - - - (5) 
with 
1 oS i 
Bi(s)= a/ (am) tbe ‘. a(t —— dr. 

It is the purpose of the present paper both to throw some light on the 
interpretation of these equations and to relate them to the cluster-integral 
theory of Mayer and Montroll (1941), as has been done recently for a 
pure, one component, fluid by Rushbrooke and Scoins (1952).* In 
particular, we are interested in the significance of the vanishing of the 
denominator in (3) or (4). For a one component fluid, when the 
corresponding denominator is 1— peA(s), it is supposed that condensation, 
from gas to liquid at a given temperature, occurs at that value of p for 
which 1—pe6(0)=0; ie. that the significant root of the denominator 
is at s=0. Support is certainly lent to this assumption by the work 
of II, and we shall assume that here too, in the case of a binary mixture, 
we are concerned only with the roots of 


(1—paeasPaals)) 1—PpespPapl$))=PaP seas Pap (Ss) - (6) 
at the point s=0. 


§ 2. CONNECTION WITH THE THEORY OF MAYER AND MONTROLL 


Analogously with the expansion of g(r) in powers of p given by Mayer 
and Montroll (1941) we have, for a binary mixture, the three expansions 


gylr)=oxpl—d (lk I] 1+ E patontatilr | eta) 


where (ij) stands for (AA), (AB) or (BB) and m and v are positive integers. 
The functions g',(r) are related to the irreducible cluster integrals 
familiar in the virial coefficient theory of imperfect gases: and it is 
appropriate to refer to them as open cluster integrals. Graphically, an 
open cluster integral involved in g7,(") has two fixed major points, 
distance r apart, occupied by species A or B according to the upper 
specification (7j), and m-+-n other points, over which spatial integrations 
are performed, of which m are occupied by A-species and n by B-species. 


* Referred to hereafter as IT. 
4R2 
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gi,(r) is the sum of all distinguishable such open cluster integrals, the 
integrand in each case being a product of the appropriate «(r) factors 
between pairs of the m-+-n+2 points (excluding any direct link between 
the two major points). In each case there must be sufficient linkages 
for the corresponding cluster diagram to form an irreducible cluster 
diagram of the type familiar in virial coefficient theory were a link to ke 
included between the two major points. 


Fig. 1 


(Z)e Gf) 


Fig. 2 


(z) (f) 


In each case there are m-+n ©’s of which m are occupied by A’s and n by B’s. 


Now it is quite impracticable, on this basis, to derive a closed 
expression for the right-hand side of (7); but there are two simple 
approximations to the functions g'3,,(r) which enable us so to do. They 
consist in replacing g/’,(r) by either c¥,,(r) or d,(r) (the latter approxima- 
tion being the better, but more elaborate) where the integrals involved 
iN Cy,(r) are simple chain integrals, fig. 1, while those involved in d’4,(r) 
correspond to netted chains, fig. 2, In neither case, of course, is full 
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justice done to g(r) except when m+-n=1; but we have shown in IT 
that for an imperfect gas either approximation gives the second and 
third virial coefficients, B and CO, correctly, while the netted chain 
approximation takes very satisfactory account of the fourth virial 
coefficient, D. The analysis for mixtures is a little more difficult, in 
that linear Ising lattice problems are involved, but equally practicable : 
it is hoped that the following rather condensed account will sufficiently 
indicate the method of derivation of the final formulae. 


(i) The Simple Chain Approximation 


According to Lemma J of I, the chain integral involving p A-A links, 
q B-B links and X A-B links is simply 


il ro) 
Tam Omiertme | sPRas\Pb(e)BEa(6) 


sin 7s 
r 


Coe ean) 


In evaluating 2 p,"p,"c¥,,(r), therefore, we require to find 
m+n>1 


ZO AaB i als) Pug@)ban(8)s0 eee - > «+ (9) 


where the sum is over all appropriate values of m, n, p, q and X. 
Consider, for definiteness, the approximation to g,,(r). Then we are 
concerned with linear arrangements of the type, fig. 3, in which 


Fig. 3 


X=2Y-+1 and, apart from the terminal blocks, there are Y blocks 
of A’s and Y blocks of B’s, each containing one or more systems. Thus 
for given Y, >0, the contribution to (9) is 


1 ( Pa at PB i e¥ +1 I mee (10) 
(I—paBaa)\1—paBas/ \l—pphsp (l—ppP sp) 
In the particular case Y=0, however, when we have terminal blocks 


only, the contribution is 


Bas £6 
Th oo ae Oe 


it being necessary to remove the contribution, Bap corresponding to 
m+tn=0. Summing (10) and (11), Y>1, we find that (9) reduces to 


ee iat EAT <P s SS e Te 
(1—paPaa)l—Pshss)—PaPshas Pan Le 
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Consequently, on this simple chain approximation, we find (combining 


(7), (8), (12) and (1)) 


a Bax(s) ] 
Janl)= atl ula —pxPy(8))l—pphppl$)) —pappPan(s) Ban(s) 
pin oe en a 


which is exactly equivalent to eqn. (4) if in (4) we put x A—=€Bp= "ABU 
In the same way, the present approximation to fxa(r) yields eqn. (3) 
with €,,—€pp=€apn=—!. In other words, if in the eo theory 
we put all the parameters « equal to unity the ensuing equations are 
just those which result from replacing, in the Mayer—Montroll expansion 
formulae, every g¥,(r) by the simple chain approximation c!,(r). 


(ii) The Netted Chain Approximation 
The corresponding formulae based on the netted chain approximation 
are derived by essentially the same argument. We require the generating 
function which replaces the expression (12). Along the length of the 
principal (shortest) chain joining the major points, themselves occupied 
by A and B systems respectively, we have a sequence of A’s and B’s 
falling, as before, into blocks, fig. 3. Let there be X=2Y-+-1 divisions 
between these blocks. Now consider one of the Y central blocks with 

A’s along the main chain: this can comprise as many 


A 


B 
Loa ane Oe 
A 


A A A 


AEA 


units as we like in any order, which entails a generating factor 


——___?s____ (14) 
1—pyBya—Pa'Yasa—PaPBYABA 


the first term corresponding to a single A (we have still to introduce the 
appropriate factor for all the A-B links in the main chain). Here 
Yapa, for example, is defined a 


you stean 
YABA=YaABA(S => ea tupiy(u u) ———— du dt, . (15) 


again using Lemma J of paper I. 
Similarly, since each A—B link in the main chain can be 
A 


B 
‘ ne a eA 


A B A B 


A 


in place of 6,3(s) in (12) we now have BaptPaYaaBtPBYABB 
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Tn fact, quite generally, changing from simple to netted chains, i.c. from 
c2(r) to d4,(r), alters our previous results to the extent that 


Ba 4 becomes BaatPaYaaatPBYABas 
By, becomes BuptPaYBaBt+PBY BBR? . og fENG) 


Bay becomes Bypt+Payaapt+Payape 


Consequently f, 4(r) and f,,(r) are now given by eqns. (3) and (4) 
with €,,4, €pp and ie replaced by €, 4(8), €p p(s) and €4 p(s), where 


€4a($)Ba a(8)=Bya(8)+Payaaa(S)+Ppyapals)s 
€pp($)Ppp(8)=Bpn(8)+PaypanlS)+Paypppls), ¢- ~~ (17) 


€4p(5)B 4 n($)=Ban(8)+Payaan(S)+PByapp(s)- 


Thus we find that the netted chain approximation, in which g(r) is 
replaced by diJ,,(r), is not quite the same as the Born-Green linearized 
integral equation theory, for there the parameters « are taken to be 
independent of s. The work of Rushbrooke and Scoins on pure fluids, 
paper II, shows that the netted chain approximation is in fact better 
than the Born-Green theory, in that a closer approach is made to the 
fourth virial coefficient of an imperfect gas. But the singularities of the 
two theories are the same, i.e. at the same values of p and 7’. In the 
one-component case these conclusions rest on the approximation (see I), 
e=1-+(3C/2B)p which has to be replaced in the present, binary mixture, 
problem by the equations 


€ya=14+(3C ga apat3Cyappp)/2Baa, 
Beams es nen ay Or Kei (18) 
€ap=1+(3Cgappat3Cappps)/2B ap. | 


Here the B’s and C’s are the coefficients in the virial expansion of the 
pressure of a binary gas mixture (see §3 below). Then at s=0 eqn. (16) 


reads 
(1+2B 4 ,pyt3C 44 apa t3C 4 apPaPs) 


(1+2Byppp+3C ,pppaPa+3Cpppps) 
=papp(2Bapt3CaapPat2Cappes), ~- - (19) 
since (27)3/2B, ,(0)=—2B,4, etc. ; and this is also the condition for the 
vanishing of the denominator in the netted chain theory, provided this 
also occurs at s=0, since 
(277)8y 4 4 a(0)=—3C yaa, (277)?y 4 pa(0)=(277)?y,an(0)=—3C yap, - (20) 


etc., the eqns. (20) following from (15) and Mayer’s virial coefficient 
theory, see I. It remains, therefore, to interpret eqn. (19). 
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§3. SIGNIFICANCE OF THE VANISHING OF THE DENOMINATOR 


We are now in a position to interpret the vanishing of the denominator 
in eqns. (3) and (4), assumed to occur at s=0, for we have seen that this 
implies eqn. (19) which is evidently related to the equation of state 


N 1 : 
p=kT [a + pe [VP Bast2N NpBaptNp Bz) 


+793 WV 2Canat 3 s2N 30 ant 8N NC appt Ws"Cnss | (21) 


Equation (21) is simply the virial expansion for a binary gas mixture 
including terms up to those in p?, but no further. From this equation 
we can determine the two chemical potentials, wu, and py (see 
Guggenheim 1952), and find 

Ma=kT Inge + = [IV Bays tNpBazl] 
A 


3kT 
“t: 2V2 [N 2C ana t2N NpCaaptN eC aps, 
. . (22) 
N 2kT 
Vz + 7 IN Bast NpBzsl 


3kT 
+ 3p [NPC sant2N NpCapptNpCzppl, 
where ¢, and ¢, depend only on 7. It is then easily verified that 
eqn. (19) simply expresses the condition 


(ma, HB) (Mas Lp) 
O(N, Vp) — (FOr ~O(p, @) 3 | Derr 3 (23) 
where p=(N,+N,)/V and a=N,/(N,+Np). 

For a given temperature, eqns. (22) and (23) together define a curve 
in the (p,2) plane. For values of the density and composition 
corresponding to points on this curve small isothermal changes, dp and dz, 
of density and composition can be found which leave the chemical 
potentials, 1, and zp, unchanged. The denominator with which we are 
concerned in the Born—Green theory, or the netted chain approximation, 
for a binary mixture vanishes, at a given temperature, for just these 
values of p and 2. 

It is easy to prove that as 2-1 eqn. (23) becomes simply dp/dp=0. 


pp=kT In 
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ABSTRACT 


Among 18 000 photographs taken during an experiment with a large 
cloud chamber at the Jungfraujoch (3 580 m) 13 examples of the decay of 
charged V-particles have been observed. A summary of the measure- 
ments of the 13 particles is given. Two of them are produced in nuclear 
interactions of remarkably low energy. 

The average time that the particles would have taken to cross the 
chamber, if they had not decayed in flight, is ~10-® sec, and the distri- 
bution of the decay events does not differ significantly from that which 
would be expected for particles with a mean life greater than 10-® sec. 


§1. InTRODUCTION 


In the course of a cloud chamber experiment carried out at the 
Jungfraujoch, 13 photographs have been obtained which show the decay 
of charged V-particles. We have found no evidence which conflicts with 
the assumption that charged V-particles are identical with the y- and 
«-mesons found in the photographic emulsions. 

The first part of this paper is devoted to the interpretation and dis- 
cussion of the 13 decay events. In § 5, the mean life of V-particles is 
considered. 

§2. APPARATUS AND RESULTS 


The apparatus used in this work is a large cloud chamber, situated in a 
magnetic field of 5 800 gauss ; the illuminated volume of the chamber is 
54x 54x9cm3. The equipment is installed on the Jungfraujoch (3 580 m) 
where the mean atmospheric pressure is 49 cm mercury. 

The photographs discussed here were obtained in about a year of 
effective operation. Three different arrangements of Geiger counters have 
been used to control the chamber expansions. Figures 1 (a) and (b) show 
the disposition of the counters and absorbers for selection systems A and B 
respectively. Selection system C is a modification of system B; the 
absorber inside the chamber has been removed ; each of the counters in 
the middle tray has been shielded from its neighbours by ¢ in. of lead ; and 
expansions are controlled by five-fold coincidences comprising two counters 
from the top tray, two from the middle tray and one from the bottom tray. 
ee ee i ee eee 

* Communicated by Professor P. M. 8. Blackett, F.R.S. 
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The absorber above the chamber with system C has been either 12cm of 
copper or 30cm of lead. We have observed no outstanding difference 
between the results with the different absorbers, but no detailed analysis 


has yet been made. 
The results obtained with the different selection systems are summarized 


in table 1. 
Fig. 1 


PS 
KX 
KS 
Kh 
RT 


(a) Selection system A. The counter immediately above the cloud chamber 
is 50 emx4 cm. All other counters are 65 cmx4 cm. A six-fold 
coincidence is required to expand the cloud chamber. 

(b) Selection system B. The counter tray above the chamber contains 6 
counters 40 emx1-5 cm. The remaining counters are 65 cm x4 cm. 
Six-fold coincidence is required. 


Table 1. Statistical Analysis of Results from Three Selection Systems 


Running Penetra- | Probable 
Selection | Time Naso ting | Penetrating Neuial  Onere 


Ph s 
Ges, 10tographs haere ih Shaan V-decays|V-decays 
A 14 5814 502 249 10 1 
B 14. 5881 783 1229 29 4 


C 12 6631 1318 1807 40 8 
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For systems A and B a shower was classified as penetrating if at least 
two of the secondary particles penetrated the 3 cm lead plate in the 
chamber without producing an electron cascade. This stringent condition 
is not satisfied by many high energy nuclear interactions above the cloud 
chamber. Table 1 gives, in column 5, the number of photographs which 
show convincing evidence for such interactions yet do not satisfy the 
penetrating shower criterion. For selection system C, with no lead plate 
in the chamber, the allocation of any event to column 4 or column 5 has 
inevitably to be left to the discretion of the -observer ; it is possible, 
therefore, that the number of events given for system C in column 4 has 
been over-estimated. 

We have observed 40 examples of tracks undergoing sudden deflections 
in the gas through angles greater than 3°. An analysis of these kinks is 
given in table 2. 


Table 2. Interpretation of 40 Sudden Deflections of Tracks in the Gas 


V+-decays 13 
m™—-decays 13 
p—-e-decays Mi 
Elastic scattering 3 
Indeterminate 10 


§3. MEASUREMENTS AND INTERPRETATION OF CHARGED V-DECAYS 


The measurements made on the photographs of the 13 events classified 
as charged V-decays are given in table 3. Photographs of events 3, 4, 
6 and 11 are reproduced in Plates 1 and 2. 


Table 3. Measurements Made on Charged V-Particle Tracks 


Primary Secondary 
______—__—__—______| Angle of =U 
i iZ- Deca oniz- 
rent. Sign Momentum Eee Mass (D ee Momentum | ation Mass 
(108 ev/e) | (7/75) (m,) (10° ev/e) | (1/75) (mm) 
(1) + >8 <2-5 _ 1542 | 7t010 | <2:5 a 
(2) — | 2-3+15% |2-5to5| 630 to 1420] 95-7 —— — — 
(3) + 15 to 7 |3to6|550to04200| 80+7 — = as 
(4) + >5 < 2-5 —- 90+5 | 0-47+8% | 4 to 7 200 to 330 
(5) er 4:2 10% | =< 2°) |. = 1490 7545 Ti <2°5 a= 
(6) — | 3-0-+15% | 2-4 | 690 to 1570) 53-2 2:8418% | <2°5 — 
(7) Q >10 25 — eae ~- is : — 
8 $ a <2°5 — 57+2 = 2- ae 
‘9 2 >4 <2:5 = 5+1 210 eae alae 
(10) -- >10 5 — Vib) |) 640%, 25 7s 
(11) 33 >4 2@ 5 — 2842 | 49410% | <2-5 idk 
(12) ae >3 <2-5 — 18+2 | 5:2+10% | <2-5 — 
2- ~ 45 +5 —- — oe 
(13) 2 >8 <2:5 


———$<— SS 
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The errors in the momenta given in table 3 are probable errors resulting 
from gas motion in the cloud chamber, uncertainty in the magnetic 
field strength, spurious curvature due to multiple scattering of the particles 
in the gas and errors in measuring the track curvature. The ionization 
values are the limits of the visual estimates made by three observers. In 
some cases the momentum and ionization values enable the masses of the 
particles to be deduced. The mass limits given in table 3 have been 
obtained by combining first the minimum values of momentum and 
ionization and then the maximum values. The secondary particles in 
events 2, 3 and 13 are travelling almost directly away from or towards the 
camera, and nothing more than their direction of travel can be determined. 
In event 3 a convection current in the cloud chamber distorted all the 
tracks immediately below the lead plate. Measurement shows that 
the distortion increased the apparent momentum of the primary particles, 
and the measured curvature is assumed to give the upper limit to the 
momentum. The lower limit has been found by estimating the effect of 
the convection current on neighbouring tracks and assuming similar 
distortion of the V-track. 

Ionization estimates have been made by comparison with neighbouring 
tracks produced by energetic particles. The work of Ghosh, Jones and 
Wilson (1952) shows that the reference tracks may have an ionization 
density as much as 30% above minimum. Bearing in mind the errors of 
our visual estimates we do not consider that we can distinguish with 
certainty tracks for which J/I)<2-5 from tracks of minimum ionization. 


Interpretation 


In interpreting the 13 events as decays of charged V-particles we have 
followed closely the arguments used by Armenteros e¢ al. (1952). In all 
cases interpretation in terms of a large angle scattering in the gas can be 
excluded by the absence of a recoil track. Events 2, 3 and 6 cannot 
be decays of light mesons because the masses of the primary particles are 
far too large. In the remaining events the momenta of the particles and 
the angles of decay are incompatible with the known dynamics of 7 and 
decays. In most cases P;, the component of the secondary momentum 
transverse to the line of flight of the primary, is much greater than the 
maximum possible in a m or mw decay (for 7, Pymax=29 Mev/c; for p, 
P tmax= 55 Mev/c). In event 13 the angle between the tracks is much greater 
than that possible in the decay of a 7 or w meson with such a high primary 
momentum. Whichever particle is regarded as the primary of event 4 
(Plate LX XVIIT) interpretation in terms of 7 or 4 decey can be excluded 
by one or the other of the foregoing arguments. This event can however 
be interpreted as the decay of an upward moving V,°-particle. No case 
has so far been reported of a V°-particle producing a slow positive light 
meson ; and V°-particles moving upwards are rarely found. In view of 
the high momentum of the lightly ionizing particle, and in the absence of 
any obvious origin fora V°-particle, we consider that it.is morereasonable to 
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interpret the photograph as representing the decay into a light meson of a 
charged particle travelling downwards in the chamber. The primary 
particle is inclined at 35° to the plane of the cloud chamber and we presume 
that it comes from a secondary nuclear interaction in the magnetic pole. 


§ 4. Discussion ofr THE DEcay Events 


The average energy of the penetrating showers causing the expansion of 
our cloud chamber is estimated at 1 to 5x10!ev. The 13 events 
discussed here were not in general associated with penetrating showers 
noticeably different from the average. However, three of the charged 
V-particles were produced in secondary nuclear interactions of the 
particles in the main shower. Two of these (events 2 and 3) took place 
inside the cloud chamber and in these two cases the total energy of the 
secondary interactions appear to have been very low (<5x10%ev). In 
these two events heavily ionizing V-particles were produced. The third 
example of the decay of a heavily ionizing particle occurred in a shower 
which does not appear to be exceptionally energetic (event 6, Plate 
LXXIX(a)). 

These observations may be important for the design of experiments to 
study charged V-particles, but it must be remembered that there is a 
strong instrumental bias against observing the decay of very energetic 
V-particles. The more energetic the particle, the less likely it is to decay 
in the chamber ; if it does decay, the smaller will be the angle between the 
two tracks and therefore the more difficult it will be to recognize and 
interpret the event. 

In two of our events (1 and 5) the geometry of the decay was such that 
the y-rays produced by a secondary 7°-meson might be expected to enter 
the 3 cm lead plate inthe chamber. In neither case is there any indication 
of an electron cascade from this cause. However the probability of such a 
cascade being observed is rather low owing to the wide angular divergence 
of the photons from the decay of alow momentum 7®-meson. If7°-mesons 
were produced in these two events we estimate the probability of our 
observing evidence for them as ~30% for each event. 


§ 5. THE MEAN Live oF V-PARTICLES 

In accordance with the recommendation of Butler and Wilson (1952) 
we present the measurements relevant to mean life determination. 
In table 4 we list the visible path lengths of the particles before decay (0), 
their possible path lengths iv the cloud chamber had they not decayed (L) 
and the corresponding times measured in the rest systems of the particles 
(tand 7’). 

The relation between ¢ and / is 

: l=(P{M) cb, 

where P and M are the momentum and rest mass of the particle, measured 


in energy units. 
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Estimated values of P/M are given in table 4, column 4. For events 2, 
3 and 6 they are made by taking the mean of the ionization estimates. 
For the remaining events they are obtained from the measured momen- 
tum associated with a value for the mass assumed to be 1 200 m,. For all 
but four of the events only lower limits can be given. 


Table 4 
l L ; t fy 

EVENS (cm) (cm) a! (101° sec) (10-2° sec) Zee 
1 6-5 21-0 21537 -<al lay <5-4 0-31 
2 13-0 33-0 0-47 9-2 23-4 0-39 
3 10-0 22-0 0-40 8:3 18-3 0-45 
4 8-0 16-5 0-82 <3°3 <6-8 0-48 
5 13-0 19-0 0-7* 6-2 9-0 0-69 
6 38:3 60-5 0-53 24-1 38-0 0-63 
7 11-1 15-0 >1-6* <2°3 <3-1 0-74 
8 20:2 41-0 >0-7* <10-4 < 21-2 0:49 
9 9-6 30-0 ==0°7% <4-9 <15-3 0-32 
10 33-7 44-5 >1-6* <6-9 <i9- il 0-76 
i 8-0 18:5 0:75 <4-l <9-5 0-43 
12 6-6 26-0 >0°5* <45 ee 7310) 0-25 
13 22:2 43-5 =>1-3* 051, oA) 0-51 


* Assuming mass of V+=1 200 m,. 


The quoted ratios ¢/7’ are, however, independent of the measured 
momenta. They are listed in table 4, column 7. All the distances have 
been measured to the surface of the illuminated volume. The correction 
suggested by Butler and Wilson has not been applied. The mean value 
of t/T' for the group is about 0-5, and no reasonable ajustment of a ‘ fiducial 
surface’ would make the mean value of ¢/7' differ significantly from 0-5. 

The expected value of 1/n 2'(t/T), for particles with a mean life much 
greater than their average time of flight across the chamber, is 0-5; our 
result therefore indicates that there may be present in our data some 


particles at least whose mean life is long compared with 7’, the average 


value of 7’. The value of 7’ depends on the uncertainty in the estimates of 
P/M, and for 9 events we have only been able to obtain estimates of an 
upper limit of their 7’-values. We can, however, get a rough indication of 
the value of the average 7’ by assuming (i) that in every indeterminate 
case the 7'-value is zero ; and then, (ii) that in every indeterminate case 
the 7'-value is given by the upper limit. We thus obtain 

(7<7T'<14) x 10-9 sec. 


It should be emphasized that, in spite of the extreme assumptions, 
the above values of the limits for 7’ are themselves subject to error. 
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Some of the estimates of P/M, although given as definite values in table 4, 
are nevertheless only known very roughly. But the main contribution 
to the average value of 7’ comes from the three heavily ionizing particles, 
for which the value of P/M has been estimated most directly. 

The mean value of ¢/T is also affected by errors in measurement. The 
length J can usually be obtained relatively accurately, but the values of L 
given in table 4 may in a few cases differ from the true value by about a 
centimetre. The effect of this error is, however, small compared to the 
statistical uncertainty in the mean value, due to the paucity of the data. 

We do not attempt here to assess any quantitative confidence limits ; 
but we consider that the most reasonable interpretation of our data is that 
the particles represented by our observed tracks have a mean life 27107 * sec; 
An investigation of the mean life by a more powerful method will be 
described in a later paper. 


§ 6. CONCLUSIONS 

Our apparatus is sensitive to penetrating showers with average energy 
between 1 and 5x 101° ev. The decay events considered as a group come 
from showers not noticeably different from the average ; but we have 
found two examples which show that V*-particles can be produced in very 
low energy secondary interactions. 

In two events where evidence for 7°-mesons among the decay products 
might have been observed it has not been found. This fact has little 
significance until combined with further evidence, as the probability of 
detection of 7°-mesons was only ~30% for each event. 

The distribution of the decay events even in our large chamber does not 
differ significantly from that which would be expected for particles with a 
mean life greater than 10-® sec. This result is consistent with those 
obtained in other cloud chamber experiments (Annis et al. 1952, Barker 
et al. 1952). 
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ABSTRACT 


The scattering of positrons by electrons has been studied in nitrogen 
by a cloud-chamber method. In 6000 photographs 5600 positrons 
were measured, giving a total track length of 792 metres between the 
energies 0-5 Mev and 1-1 Mev. The experimental results, when compared 
with the exchange theory of Bhabha, show satisfactory agreement. 


§ 1. IntRopuUcTION 


Tue theory of electron-electron scattering as formulated by Moller 
(1931, 1932) is generally accepted for comparison with experimental 
results. The theory has been developed taking into account the effects 
of exchange and of retardation in collisions between two free electrons. 
Because of the fact that sources of fast electrons are easily available, 
considerable experimental work has been done on this scattering problem 
by many workers, notably by Champion (1932), Hornbeck and Howell 
(1941) and Shearing and Pardue (1942). The results of all these workers 
showed agreement with the theory of Moller. In the case, however, of 
more recent experiments performed by Groetzinger, Leder, Ribe and 
Berger (1950), the results agree equally well with the theories of 
Moller, Mott (1930) and with Rutherford’s relativistic formula, having 
regard to the statistical fluctuation. 

The investigation of positron—electron scattering is not so extensive, 
due, perhaps, to the difficulty in obtaining suitable positron sources, 
but in some respeets it is more interesting than the electron—electron 
scattering owing to the fact that it is possible to distinguish between the 
incident positron and the secondary electron in a cloud-chamber by 
applying a magnetic field. . 

The effect of exchange phenomena arising out of the interaction of a 
fast positron with an electron when a positron passes through matter 
was first pointed out by Bhabha (1936). Due to this exchange, the number 
of fast secondary electrons ejected in the collision process between the 
positron and the free electron will be more than would be expected if it 
were a simple scattering process in which the positron is considered as 


an independent positively charged particle. 


Rees ee ee ee 
* Communicated by the Author. 
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The formula of Bhabha, after von Ritter, Lieseberg, Maier—Leibnitz, 
Papkow, Schweiser and Bothe (1951), can be written as 


Wy, €) dl de=(e?/m ec?) N dl de(P AP ytP yt Dy—-1)P, . . 
where Pj=2y*/P—2y*— Net (y—1, 


Pus () {y+ 2y+3—2(y—1)e(1—e)}, 


Lire 2(25) {y(y+2)/e—2(—1) + (y—1)°e}. 


Formula (1) gives the probability, as calculated by Bhabha, that a 
positron with total energy ym,c* transfers the ratio « to «+de of its 
energy to an electron, for a track length dl in a gas with N number of 
electrons per unit volume. 

The term P, coriesponds to the ordinary scattering of the positron 
when considered as an independent positively charged particle, according 
to the Dirac hole theory, while the effects of exchange are taken into 
account in the terms P,, and P;. The physical significance of P,, is 
that the positron in interaction with the atomic electron annihilates only 
to create a new pair of electrons of opposite charge. P,; gives the 
interference between the ordinary scattering term, P;, and the pair- 
production term, P;;. 

The aim of the present experiment is to study the scattering of 
positrons by electrons and to compare the results with the exchange 
theory of Bhabha. 


§2. EXPERIMENTAL ARRANGEMENTS 


A fully automatic cloud-chamber was employed in this investigation. 
The positron source "Rh was kept inside a curved container to allow 
only high energy positrons to appear in the chamber, which was filled 
with nitrogen at near-atmospheric pressure, and to this a mixture of 
water and alcohol was added for the formation of good tracks. A pair of 
coils supplied the magnetic field used to determine the energy of the 
tracks. The strength of the field was 739-6 gauss and was uniform 
within 0-4°% across the 24cm of the diameter of the chamber. About 
6000 photographs were taken, out of which, 5 600 tracks, between the 
energies 0:5 Mev and 1-1 Mev, were selected for measurement. The 
number of tracks per expansion was kept to a minimum in order to 
ensure that the error due to spurious effects remained negligible. These 
effects, which may give the impression of secondary electrons originating 
from primary positron tracks, may arise from the crossing of the track 
after reflection from the wall of the chamber or from the appearance of 
a 6-track near the primary positron caused by annihilation radiation 
from the source container, 
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§3. MEASUREMENTS AND CALCULATIONS 


To determine the energies of the interacting particles the tracks were 
first reprojected to their original size. The measurements were then 
taken of the curvatures before collision, and of either the ranges and 
curvatures, or the curvatures and angles, depending on the type of 
energy-transfer, after the collision. The curvatures were measured by 
fitting circles of known radius through the beginning and end of the 
track so that there was a minimum mean square difference between 
the two. The maximum error involved in this method was estimated 
to be 5%. When the energy-transfer was very great or very small the 
energy of the positrons or electrons after impact had to be ascertained 
from their range, the minimum value of which taken into account 
was 15mm. This lower limit was decided having regard to three 
. considerations: (a) that, according to the theory of Bhabha, the 
electrons are free; consequently the binding energy of the electron 
must be negligible compared with the energy-transfer; (b) that the 
determination of the energy will be inaccurate by a higher percentage 
as the range is diminished; and (c) that in fixing the lower limit at a 
comparatively high range it is sometimes possible to make a counter- 
check on the energy-transfer by measuring the change of curvature. 
The estimated error involved in the range measurement is about 9%. 

Accurate determination of the energy of the primary positron and of 
the positron and electron after collision indicated, within experimental 
error, that the energy was conserved between the interacting particles. 

Since the photographs were taken by one vertical camera, some of 
the scattering events may be altogether missed from the measurements. 
It is, however, possible to correct the observed number in the following 
way. Let 0 be the angle of scattering in the laboratory system, A its 
projection on the horizontal plane, and ¢ the angle made by the horizontal 
plane and the plane of scattering. We then obtain . 


cos é=tan A/tan 0. 
By limiting the observation to A>Aq ¢ is also limited to 
d< are cos (tan A,/tan #)=¢o. 
Therefore the probability of not observing a scattering is 
1— 249]. 


For this reason it is necessary to multiply the number of events measured 


at angle 0 by 


7 


~ 3 
Again (see Bhabha 1936, pages 198 and 203), 


2e 1/2 
wnt={Tp asa) 


482 


Ce =r] {2 are cos (tan Ap/tan 6) }. 
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C(0)=7 / 2 arc COS [ {eet i tan | \ ; 


where y=E/m,c? and « is the ratio of the energy transfer. 


Therefore 


§4. ResuLtts anD Discussion 


From eqn. (1) probabilities were calculated for the energy transfer 
to lie in the regions 0-03-0-1, 0-1-0-2.... etc., with y in the intervals 2-0-2-4, 
2-4-2-8, 2-8-3-0. 

Three histograms, figs. 1 (a), (b) and (c), have been plotted with the 
number of collisions as ordinate and ¢ as abscissa. 312 metres of track 
were measured between the interval 0:5 Mev and 0-7 Mev, 270 metres 
between 0-7 Mev and 0-9Mev and 210 metres between 0-9 Mev and 
1-1 mev. The experimental results, given by the thick line, and the 
theoretical results, the dotted line, show satisfactory agreement. 
Calculation made from the formula with the term P, only gave results 
which did not accord with those of the experiment, hence confirming 
that, in addition to the ordinary scattering of positron and electron given 
by the term P,, the exchange effect plays a part in the collision process. 
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Experimental results of the positron-electron scattering are given by the thick 
line, whereas the dotted line shows the values as calculated from Bhabha’s 


theory. 
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The first experiment to investigate the Bhabha exchange was 
performed by Ho ZaH-Wei (1946) with energies of positron varying from 
0-025 mev to 0-8 mev and a total track length of 170 metres. She 
observed discrepancies between the value, as calculated from the theory, 
and the experimental value in the region «>0-5. 

Results published recently by von Ritter, Lieseberg, Maier-Leibnitz, 
Papkow, Schweiser, Bothe (1951) on positron-electron collisions where 
positrons of low energies—from 0-1 Mev to 0-4 Mev—were used do not, 
however, confirm these discrepancies. 

Similarly, in the present experiment, in which the positrons were of 
high energies, where the Bhabha exchange is important, we noted no 
disagreement with the theory. 
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ABSTRACT 


The technique for obtaining aligned Cobalt nuclei has been applied to 
the isotope Cobalt 58. An examination of the angular distribution of 
the 0-805 Mev y-ray shows that this transition cannot be dipole but is 
probably quadrupole. The temperature dependence of this angular 
distribution gives for the nuclear magnetic moment a value of 3-5-0-3 


nuclear magnetons. This result is discussed in relation to the nuclear 
shell-model. 


Tue technique for the alignment of Cobalt nuclei has now been established 
(Daniels e¢ al. 1951, Gorter et al. 1951) and demonstrated for the radio- 
active isotope Cobalt 60. We have used the same technique to align 
the isotope Cobalt 58 in its ground state. The disintegration scheme of 
this nucleus has been investigated by Deutsch and Elliott (1944) and by 
Strauch (1950) (see fig. 1). They showed that decay was by K electron 
capture (~ 85%) and by positron emission (~ 15%) to an excited state 
of Fe 58, and thence by a single 0-805 Mev y-ray to the ground state. 
Strauch concluded from measurement of the internal conversion co- 
efficient «, that this y-ray was M1: however, it seems unlikely that this 
. measurement could exclude #2 which is more probable on the basis of 
existing knowledge of the first excited states of even—-even nuclei. Since 
there is only one y-ray and the positron spectrum has allowed form, no 
further information can be expected from angular correlation measure- 
ments. However, a measurement of the angular distribution of y-rays 
from aligned nuclei can certainly distinguish dipole from other multi- 
pole transitions (Steenberg 1951). Hence this is a case where nuclear 
alignment offers special advantages. 

The experimental technique was the same as that employed for 
Cobalt 60 except that scintillation counters fitted with lin. cube 
NalI(T1) crystals were used to discriminate against annihilation radiation, 
scattered 0-805 mev radiation, and low energy y-rays from Cobalt 57 
which was present as weak contamination. 

The fact which was immediately observed, that the y-rays are emitted 
preferentially at the nuclear equator, shows that the y-ray transitions 
cannot be dipole. From the presence of coincidences between x-rays 


* Communicated by the Authors. 
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and y-rays (Deutsch e¢ al. 1944) and from Strauch’s value for ag we 
conclude that the radiation is quadrupole. It follows from the allowed 
form of the positron spectrum that the angular momentum carried away 
in the B transition, i,+1. This allows the spin assignments to the 
levels to be made in three possible ways :— 

(A) Ig=2, [,;=2, 1g=0, 

(B) 1g=3, I,=2, 1=1, 

(C) Ip=2, [,;=2, t=1, 
where I,, J, are the spins of the nuclear levels before and after the 
K-capture process. 


Fig. 1 


85-5 %K. 
14-5%8* O47 


Measurements of y-ray intensity (W) were made along the K,-axis 
(W(0°)), the K,-axis (W(90°)) and along the bisector of these (W(45°)). 
(The K, and K, axes are the axes of maximum and minimum magnetic 
susceptibility of the Co++ ions in the salt used.) For quadrupole radiation 
the polar diagram at any temperature (Spiers 1948) 7’ is known to be of 
the form 

W(@)=1—A(T) cos? 0—B(T) cos? 6 
and three different expressions for the functions A and B can be obtained 
theoretically for the three different spin assignments given above. 

The magnitude of the anisotropy 

ais W(90°)— W(0°) 
i A(T)+B(T)= due (00°) saa 
was obtained from measurements of the y-ray intensity, and the results 
are shown in fig. 2, where 1/e is plotted against 1/7’. On the same diagram 
the three theoretical curves are plotted, fitted at the point 1/7=52. 
The 1/7' scale of these curves is arbitrary by a factor proportional to the 
magnetic moment of the parent nucleus ; this factor is chosen so that 
the curves fit the experimental results at the point 1/7=52. It is seen 
that the best fit over the whole range is obtained from Curve C. Further 
confirmation of the correctness of this spin-assignment is obtained by 
determining experimentally values of the function 


B(T)= {4W(45°)—2W (0°)—2W(90°)}/W (90°). 
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; These are shown in fig. 3, along with the three theoretical curves, and 
it is seen that the experimental results again favour curve C. Assuming 
the decay scheme C to be correct, we find, from the fit of the theoretical 
curve to the experimental points, that the magnetic moment of Co 58 
is 3-5+0-3 nuclear magnetons assuming that no reorientation of the 
nucleus occurs following the K-capture process. 


Fig. 2 
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ANISOTROPY OF y-RADIATION FROM Co®®Rb,(SO,)3.6H,O 


Fig. 3 


The nuclear shell model (Mayer 1950) suggests that in the ground state 
of both Co 58 and Co 60, the odd proton and odd neutron are in f7/2 and 
p3/2 states respectively. Pryce (1952) has shown that the coupling 
energy is greatest, and nearly equal, for the parallel and antiparallel 
configurations. These give states of 5+ and 2-+- respectively. 
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This model satisfactorily accounts for the established features of 
Cobalt 60, namely the isomeric transition with 4J=3, the ground state 
spin of 5 deduced from the decay scheme, and the magnetic moment 
(Bleaney et al. 1952) of the ground state, 3-5 nuclear magnetons, which is 
close to the Schmidt limit of 3-8. In our previous paper (Bleaney e¢ al. 
1952), owing to an error, the magnetic moment was given as 3:0-+0:5 
instead of 3-5-+-0-5 nuclear magnetons. 

In Cobalt 58 also, the ground state is fed by an isomeric transition with 
AI=3, and our assignment of spin 2 to the ground state is in keeping 
with one of the alternatives offered by the shell model. However, the 
Schmidt limit for the magnetic moment of this state is 6-2 nuclear 
magnetons, which is considerably larger than our experimental value 
of 3-5. It has been suggested (Grace and Halban 1952) that some 
reorientation of the nucleus may take place after the K-capture process 
due to the large field of the odd K electron. A detailed examination (to 
be published) of this effect shows that it is unlikely to account for this 
difference. 


We wish to thank Professor Lord Cherwell for his continued interest 
and encouragement. We are grateful to Professor F. E. Simon for 
many helpful and stimulating discussions. We are indebted to 
Dr. W. I. B. Smith and the Nuffield Cyclotron Group at Birmingham 
for their cooperation in the irradiation of the source, and to Dr. P. F. D. 
Shaw of this laboratory for the chemical separation of the Cobalt 58. 
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ABSTRACT 


CdI, crystals growing in thin plates, have been subjected to mechanical 
deformation, by an indentation process. As a result of this deformation, 
the growth properties of the crystals are altered, and several other striking 
effects take place, which have not previously been observed. It is 
concluded from the experiments that the mechanical deformation is the 
vital factor for these growth effects to occur. Variation in concentration 
by disturbing the mother liquor has no effect. Slip processes, leading to 
the creation of dislocation or increase of dislocations already present in the 
crystal, may account for the spiral nature of the growth sometimes 
observed. 


§ 1. INTRODUCTION 


CdI, CRYSTALS growing from aqueous solution have been studied by Forty 
(1951, 1952) and the spiral features as well as the step heights are in support 
of Frank’s hypothesis of growth by a dislocation mechanism. When 
CdI, grows from solution one may also observe thin crystal plates of an 
extremely uniform thickness growing in extension to about 4mm with 
little or no observable change in thickness. Buckley (1952) has drawn 
attention to these plates and raised the question of how this correlates 
with the spiral growth. 

It was thought that these plates may have nearly perfect (0001) planes 
and if exposed ledges were produced on these planes they might grow faster 
in the [0001] direction. 

Experiments were performed to see whether an external force could 
deform the crystal in such a way that exposed ledges were formed and 
rapid growth started. Although the method of deformation used has not 
been developed to give controlled stresses, the results obtained are very 


informative. 
The experimental technique will be briefly described and the results 


discussed. 
§2, EXPERIMENTAL TECHNIQUE 


The crystals were grown in a drop of solution placed on a microscope 
cover slip. This was placed on the moveable stage of a Cooke-Troughton- 
Simms microscope, so that it could be observed in transmission and 


reflexion. 
* Communicated by Professor S. Tolansky, F.R.S. 
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An indenter was made by drawing out a glass rod to a fine point. This 
was adjusted above the cover slip in the centre of the field of view. It. 
could be raised or lowered by a vertical screw motion. In this way the 

_ transmitted light was not affected by the presence of the indenter when the 
latter was not in contact with the crystal. As the indenter was lowered it 
came into the field of view as a shadow, so that the position of contact of 
the indenter could be controlled, and the complete behaviour of the crystal 
before, during and after indentation could be observed. 

A Leica adaptor was used in place of the eyepiece with a 35 mm Leica 
camera body, so that the process of growth could be observed simultane- 
ously while photographing. 

In general a crystal growing as a flat plate was selected and a photograph 
taken. The indenter was very carefully lowered until slight contact was 
made and then carefully raised out of the field of view. Often no growth 
was observed at the point of contact for a few seconds and this emphasizes 
the slight amount of deformation caused. Precaution was also taken to 
see that the indenter was clean. As growth proceeded a series of photo- 
graphs was taken at intervals. 


§3. RESULTS 


Several of the crystals were observed over a long period of time and they 
remained as flat plates increasing in size laterally. Buckley (1952) has 
observed this for several hours. The neighbouring liquid was disturbed 
but this had no effect on the crystals, as will be seen later. Only a few 
seconds after the indenter had touched the crystal more rapid growth in 
the [0001] direction was observed. The growth features will be classified 
into three groups : 

(a) growth originating from the indent and along visible cracks which 

may develop during the indentation ; 

(6) growth originating at points far from the position of the indent and 

with no visible markings connecting it with the indent ; 

(c) development of the (1120) faces of an incomplete crystal to the 

first order (1010) close packed faces forming the stable hexagon. 

The growth features of groups (a) and (b) were very similar to the features 
found on the thicker crystals growing in the solution and very often spirals 
were clearly seen. One striking difference is that an intermediate state 
exists, where the growth initiated at certain points has covered part of the 
original crystal with a rather thick layer of new material, but has left other 
parts of it unthickened. This point will be discussed again later. 

The observations belonging to group (c) are very numerous, but no 
information about the mechanism leading to this growth phenomenon has 
been obtained from them. 

The photographs given here show some of the interesting features, and we 
shall describe these in some detail. The first series, figs. 1-8 (Plate LX XX, 
figs. 1-4; Plate LX XXI, figs. 5-8), gives an example of a hexagonal 
plate, which before the indentation, fig. 1, showed no features on the 
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basal plane. The kink in one edge is merely a growth irregularity and 
nothing exceptional happens at this place. Figure 2 is a photograph in 
reflexion of the crystal only a few seconds after indentation. Small 
hexagonal layers are formed at the point of indentation, with the same 
symmetry as the large crystal. The crystal shows uniform tint over 
the whole area except where the hexagonal layers are forming. 

In the next photograph, fig. 3, a spiral is seen originating at the point 
where the indentation was made, and a thick hexagonal plate of the same 
symmetry as the main crystal is seen to have developed on the opposite 
side. 

Figure 4 is a reflexion picture in which the spiral is visible, and more 
complex growth features. The interference fringes show how the thick- 
ness of the crystal is varying, though away from the indent it is still of 
the same initial uniform tint. 

More complex growth patterns with a number of growth hills at the 
crystal edges are seen in fig. 5. Also it may be observed how the kinks 
at the edges appearing in the earlier pictures have moved, tending to 
complete the sides. 

Figure 6 shows in transmission that the plate formed is now very 
thick, but there is no sign of any growth features belonging to the exposed 
basal plane of this plate. 

The reflexion picture, fig. 7, shows how the spiral mechanism is working 
at two of the corners of the hexagon. The part where broad fringes 
still occur must be growing slowly upwards and a filling-in of these areas 
from the edges is now taking place. 

In the last picture of this series, fig. 8, the growth pattern has become 
very complex. The spiral is still visible, and the hexagon plate has 
increased considerably in thickness. 

The next series is an example of the observations on an incomplete 
hexagon. Figure 9 (Plate LX XXII, figs. 9-14) shows the crystal plate 
before indentation. Crystals of this type occur frequently and maintain 
this shape during growth. The long curved edge is unstable on indent- 
ation, however, and growth phenomena have been observed on this edge, 
which have not been seen on any crystal which was left untouched. 

In the following figure (fig. 10) indentation has been made near to the 
edge. Growth has taken place round the indent and there is a thickening 
of the edge of the crystal at this point. 

Figure 11. .A second indentation has now been made in the centre 
of the same crystal. A crack, which appears to be discontinuous, extends 
downwards towards the long edge. The growth by the first indent has 
proceeded, but the thickened edge has begun to dissociate. Further 
along the edge there is a. kink, and crack from which a layer has started 
to form. 

These features have become more marked in fig. 12. The crack from 
the second indentation has formed a layer one edge of which has taken 
up the crystal symmetry. The small facet where the curved edge meets 
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the (1010) face started to increase after the second indentation. It has 
now become a definite face, so tending to form a complete hexagon. 
The original thickened edge has thinned and dissociated. At the kink 
the edge is also tending to form the edge of a complete hexagon, and the 
layer growing here has taken the symmetry of the crystal. 

In fig. 13 the above features have continued to develop, and the layer 
coming from the edge due to the first indentation has met the visible 
end of the crack. At this point a very rapid growth is initiated, and as 
can be seen in the next figure (fig. 14), this seems to have a spiral nature. 

In this picture, under higher magnification (300) than the others 
(x 150), the crystal on the right is the one previously indented. The 
growth at the crack may be noted to have occurred at a far greater rate 
than any other feature on the crystal. The crystal on the left may be 
seen in all the previous pictures to have been unaffected by the distur- 
bances around it. The little facet at the right point has increased 
slightly during the growth. 

This crystal was now indented and the facet may be seen to have 
grown remarkably in size, and new growth has been started by the 
indentation process. The crystal has taken the normal hexagon shape. 
A kink similar to the one formed in the other crystal has occurred on the 
lower side of this one, though in this case it is seen to be joined to the 
main indentation by some crack or flaw. It is worth noticing that 
although this connecting disturbed region in the crystal is visible, there 
is no growth taking place along it. 

The remaining pictures are individual ones taken out of a series. 
In each case the crystal before indentation was a flat plate showing no 
layer growth. 

Figures 15 and 16 (Plate LX XXIII, figs. 15-20). This crystal was 
initially as in fig. 15 and on indentation the long unstable edge changed to 
the shape as in fig. 16, where all the edges are of the stable hexagonal type. 

Figure 17. The photograph shows the thinning effect of the unstable 
edge and building up of the hexagon edges very beautifully. The edges 
completed themselves later and layer growth started at the re-entrant 
point. 

Figure 18 shows the growth at the boundary of a hexagonal crystal 
plate away from the indentation, which can be seen at the top in the 
photograph. Growth around the indentation has hardly started. 

Figure 19 shows the initiation of growth on a crystal in contact with 
the indented crystal. The crystal is of the unstable type, and the 
unstable edge is in contact with the indented crystal. Growth has 
started at this point. Hexagonal layers have formed around the 
indentation. The crystal on the left of the picture has also been indented 
and this shows features typical of those already described. 

Figure 20 shows a crystal where screw dislocations are promoting 
growth at the point of indentation. It is very significant that the 
affected area is limited and its boundary crystallographic. 
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$4. Discussion 


The discussion will be confined to growth of the three types mentioned 
previously. The following points are observed on the hexagonal plate 
crystal (figs. 1-8). A spiral has formed on one side of the plate and a 
hexagonal plate on the opposite side. The point of emergence of the 
spiral is very nearly at the centre of this small plate. If the dislocation 
group responsible for creating the spiral is extending right through the 
mother crystal, it may be that no spiral features are seen on the small 
plate because dissociation into dislocation of lesser strength has taken 
place. The rate of growth in the [0001] direction on the small plate is 
relatively greater than that of the face having the spiral. 

The great number of irregular lines of discontinuity which are observed 
in fig. 8 may be seen to develop higher visibility during growth, and their 
shape is changing only very slightly. A line of discontinuity is most 
likely marking out a surface of misfit, a fault surface, due to presence of 
partical (imperfect) dislocations. If so, it means that fault surfaces 
are created by indenting the crystal. Their visibility is increasing as 
growth proceeds because they form barriers for the advancing layers. 

The existence of fault surfaces extending to the boundary of the crystal 
provides a possible explanation to the observed growth at the edges 
separated from the indent by ‘ unaffected area ’ (fig. 18). 

The spiral growth caused by the indenter does not spread over the 
whole crystal face. This effect is pronounced in fig. 20. Any exposed 
ledge running to the edge of the crystal might be expected to sweep 
over the whole (0001) face, unless obstacles, lines of discontinuity, or an 
equal number of positive and negative dislocations are present. Lines 
of discontinuity are seen in figs. 4-8. In several figures the boundaries 
of the growing region have crystallographic directions, which indicates 
that they are not discontinuity lines. 

In several cases crystals with a curved edge have been indented. The 
main behaviour was always the same. The curved edge at once started 
growing out tending to complete the crystal symmetry. The deformation 
of the crystal by the indenter has caused slip to occur, which is trans- 
mitted to the edge forming steps such that growth in thin sheets parallel 
to the (0001) plane takes place. Figure 17 shows very strikingly how 
the crystal is growing out, the original edge being dissociated into very 
thin layers with almost invisible edges. These layers catch up on each 
other in the directions of close packing and finally complete the crystal 
symmetry. This is also seen in fig. 12. Noe 

It has been pointed out that in fig. 14 a connecting line between the 
indent and the edge phenomenon is visible. In other cases the defor- 
mation caused by the indentation is transmitted to the edge of the 
crystal without any observable change of the crystal material in between. 
Figure 18 illustrates how growth starts at the edge and appears to be 
promoted by screw dislocations, 
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Buckley (1952) says that these crystal plates were so thin that they would 
fold up in the presence of concentration currents, or of other thicker 
crystals. In order to see whether the observed effects could be due to 
disturbances transmitted through the mother liquid, the neighbouring 
crystals were always kept under observation. In the figs. 9-14 one 
neighbouring thin crystal is seen to be unaffected until this very crystal is 
indented. 

Only when neighbour crystals were actually in contact with the 
indented crystal as in fig. 19, growth is initiated so that the mechanical 
deformation seems to be the controlling factor. 

It is admitted that very little can be said about the actual deformation 
process in these experiments, where rather a rough indenter has been 
used, and it is not possible to say whether Frank’s (1951) “ buckle, 
followed by shear mechanism ’’, is in operation. Simple or complicated 
slip processes may well account for the observed facts, and uneven slip 
on a slip plane will provide the screw dislocations for spiral growth. 

Further experiments are in progress with other crystals, some of which 
have not yet been proved to grow by the spiral mechanism. 
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ABSTRACT 


A description is given of experiments which confirm in detail the 
mechanism proposed by Rollin and Simon concerning the behaviour of 
the helium film in pumping tubes of low temperature cryostats. In 
particular the temperature distribution along these tubes which, as a 
result of the film, is of an unusual nature, was determined. Experiments 
were carried out both on tubes having uniform cross sections and tubes 
containing a constriction, and it is shown how the above mechanism 
can be used to measure film transfer rates in each case. The results are 
finally discussed from the point of view of the design of helium cryostats. 


$1. INTRODUCTION 


THE existence of a helium film on surfaces in contact with liquid helium 
below the lambda-point was first suggested by Rollin and Simon (Rollin 
1936, Rollin and Simon 1939) in order to explain the high rate of 
evaporation from vessels kept at constant temperature by pumping off 
the helium vapour, and also to account for the abnormally high warming- 
up rates of such containers when the pump was turned off. They also 
pointed to a possible connection between this behaviour and the puzzling 
‘ distillation ’ phenomenon observed by Kamerlingh Onnes (1922). They 
proposed the following mechanism. There exists a film of liquid helium 
on the walls of the pumping tube which evaporates at its warm end and 
which is being continuously replenished from the bulk liquid by flow in 
the film. This mechanism could account not only for the high rate of 
evaporation but also for the high warming-up rates, as, with the pump 
shut off, the vapour would be forced to condense on the colder bulk 
liquid and warm it up. The well known fundamental experiments of 
Daunt and Mendelssohn (1938) proved directly the transfer of liquid 
helium in a film and fully vindicated the mechanism proposed by Rollin 
and Simon. 

The main purpose of the present experiments was to test in more 
detail this explanation and some of its consequences. In particular, the 
question of the actual point where the helium film evaporates in the tube 
has not hitherto been considered. Ifthe proposed mechanism is the correct 
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one, and if one assumes the pumping tube to be everywhere at the same 
temperature as the liquid helium film on its inside, then since no appreciable 
temperature differences can exist in the film (Kikoin and Lasarew 1938), 
there should be no temperature gradient along the lower part of the 
pumping tube. The film will exist up to the point where the heat influx, 
which is mainly due to conduction along the material of the tube, is just 
sufficient to evaporate the helium creeping up in the film. Taking, for 
instance, a German silver pumping tube of 0-1 mm wall thickness whose 
upper end is at 4-2°xK, one finds, using Berman’s (1951) value for the heat | 
conductivity, that the film should evaporate about 3mm below the warm 
end of the tube. 


Fig. 1 
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Experimental cryostat : pumping tube of uniform cross section. 


Accordingly the experimental arrangement was such that, besides 
measuring the rate of evaporation and the rate of warming up, it was 
possible to obtain information on the temperature distribution along the 
pumping tube. 


§2. APPARATUS AND EXPERIMENTAL METHOD 


Figure 1 shows schematically the general arrangement. The experimental 
cryostat C and the liquid helium container L into which helium was 
liquefied by a small Linde liquefier was situated in a vacuum case 
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surrounded by liquid hydrogen. The liquefier and all its auxiliary 
apparatus, which comprised a small refrigerator compressor for circulating 
the helium, has been described elsewhere (Ambler 1952). 

C could be filled with liquid helium from L through the bellows sealed 
valve V, a method which has several advantages over condensation. It 
is more efficient, it is much quicker and the liquid helium introduced 
into Cis clean. This last point is important for these experiments because, 
as has been shown by Bowers and Mendelssohn (1950), impurities intro- 
duced by condensing helium may increase the creep rate considerably. 
The leakage through valve V with liquid helium at 1 atm pressure in L 
and vacuum in C was about 0-01 ¢.c. of gas (N.T.P.) per minute, and this 
was small compared with the evaporation rate from C of between 0-5 
and 3¢.c. per minute. When L was cooled below the A-point the leak, 
of course, became prohibitively large and, therefore, L was never allowed 
to cool below 2-2°. 

S is a copper shield attached to the thick-walled copper tube K so that 
it was always at the temperature of L. Its function, was to reduce heat 
influx into C through residual gas in the vacuum jacket so that the rate 
of evaporation from C was entirely determined by the creep phenomenon. 
This was proved by measurements of the warming-up rate above the 
A-point and also by the fact that the rate of evaporation was the same 
whether S was at 4:2° or 2-3°. B and T are two iron alum crystals used 
for measuring temperature by the magnetic method; they were so attached 
to the apparatus that B was always at the cryostat temperature and 
T at the temperature of the German silver pumping tube G at a point 
4mm from the top. ; 

In order to obtain a reliable calibration curve for T it was essential 
that during the calibration no temperature gradient existed along G, and 
to ensure this the following procedure was adopted. After C had been 
filled with liquid helium V was closed and the temperature brought to 
the various calibrating points by pumping L instead of pumping C as is 
usual. By throttling the pump the temperature of L was held constant 
at each calibrating point and the temperature of T and B quickly reached 
the same steady value. In every other respect the procedure was the 
same as that described by Hull (1947). 

After the calibration L was allowed to reheat to 4:2°kK and C was 
pumped with a diffusion pump. The final temperatures of T and B and 
the steady evaporation from C were measured. The pump was then shut 
off and temperature determinations were made on T and B at regular 


intervals during the subsequent warming up. 


§3. THE EXPERIMENTS 
(a) Pumping Tube of Uniform Diameter 


When pumping the cryostat through a 45 mm long German silver tube, 
of 3mm diameter and 0-1 mm wall thickness, the temperatures indicated 
by B and T (see fig. 1) were 0-92° kK and 0-93° K respectively. This small 
apparent temperature difference is, however, within the experimental error, 


4T2 


1310 E. Ambler and N. Kurti on Film Flow 


and the result therefore proves that the conclusion mentioned in the 
introduction about the temperature distribution along the pumping tube 
is essentially correct. 

The steady evaporation from the cryostat was 3-5 c.c. gas N.T.P./min. 
which can wholly be accounted for by film flow up the 3mm diameter 
pumping tube. 

The way in which the temperatures of B and T increased with time 
after the pump had been shut off is shown in fig. 2. Up to 1-9° kK the 
temperatures of the two were equal and then T rose more quickly. The 
reason for this is that on approaching the A-point the rate of film transfer 
decreases and consequently the heat of evaporation which must be 


Fig. 2 
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Warming-up curves for cryostat shown in fig. 1. Curve B—temperature 
of cryostat. Curve T—temperature near top of the pumping tube. 


conducted down the pumping tube becomes smaller. Therefore the point 
where the film evaporates recedes down the tube, and when it falls below 
the point of attachment of T the temperature of T rises above that of B. 

From the rate at which the cryostat temperature increases the rate of 
film transfer on the walls of the pumping tube can be deduced. According 
to the mechanism proposed by Rollin and Simon, after the film evaporates 
part of the helium condenses back into the cryostat (say at a rate dm,/dt), 
causing it to heat up, and part fills the dead space (say at a rate dm,/dt) 
so that the pressure there remains equal to the vapour pressure of the 
helium in the cryostat. dm,/dt is calculable from the volumes and 
temperature distribution in the dead space and the rate of increase of 
pressure ; dm,/dt is then obtained from the expression 


dT, 
di 


dm, 


dm dm 
H—<=T, AS (oF ee TE) rie Ding A pet WSO oe 
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where H and 7, are the thermal capacity and temperature of the cryostat 
and its contents respectively, 7,48 is the ‘ fountain’ heat, Z the latent 
heat of vaporization of helium, C, the specific heat of helium gas and 
T; the temperature to which the gas is heated before condensing back. 
It is not possible to determine 7’,, so in using the above formula we have 
adopted the simple and plausible assumption that 7,—T.. 

The transfer rate R is simply given by pat 


and its values are shown in table 1 together with the data used in the 
calculation. The results are plotted against temperature in fig. 3, curve (b), 
which also shows the rate on glass as measured by Daunt and Mendelssohn 
(curve (c)) and the measured rate while pumping (point (a)). 


Table 1 


dm, H aT’. H aT dm, | Transfer rate 
dt at a Rc. 


“K/min. | g/min. | cal/min. _1T ag um g/min. em~! sec"! 
q dt 


(x 104)} (x 10) (x 103) (x104)} (105) 


2-4 4-35 
; 4-2. 
2-9, 

2-4, 
1:78 
0-95 
0-82 
0-76. 


2-60 
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1-65 
1-15 
0-85 
0-70 
0-66 
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Mass of helium in cryostat =0-225 g. 


The fair agreement between (a), (b) and (c) provides quantitative 
confirmation of the interpretation of Rollin and Simon; conversely we 
may say that, using their interpretation, the warming-up process can be 
used to measure transfer rates to a reasonable accuracy. It is worth 
mentioning that, if contrary to our above assumption, 7';>7',, the effect 
is to bring curve (b) into better agreement with (a) and (c). The dis- 
crepancy, anyway, is not significant in view of the higher rates that are 
frequently observed on metal surfaces (Mendelssohn and White 1950). 

Having proved quantitatively the correctness of the Rollin-Simon 
interpretation, we decided to make similar measurements on pumping 
tubes containing a constriction. 

(b) Pumping Tubes with Constriction 

In general the effect of a constriction is to reduce the rate of evaporation 
while pumping the cryostat, so that a lower final temperature can be 
produced (Blaisse, Cooke and Hull 1939), 
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The arrangement is shown in fig. 4, where P is a 0-5 mm diameter German | 
silver capillary soldered into the top of the cryostat C and then burnished. 
The pumping tube G had a diameter of 6 mm and wall thickness 0-1 mm 
and was also made of German silver. In every other respect the arrange- 
ment was the same as in the previous experiment. 

It was found that the rate of evaporation while pumping was only 
reduced to two-thirds of the value found in the previous experiment, 
indicating that the transfer rate per cm perimeter through the capillary 


Fig. 3 
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Transfer rates for cryostat shown in fig. 1: (a) steady evaporation rate, 
(6) calculated from warming-up rates (see table 1), (c) measurements of 
Daunt and Mendelssohn (1938). 


was about four times that for the 3mm diameter tube. Such anomalous 
rates through small apertures have been observed previously by Hudson, 
Hunt and Kurti (Kurti 1948) and by Brown and Mendelssohn (1950). 
With this evaporation we should expect the film to evaporate about 1 em 
from the top of the tube G, i.e. below the level of T, so that the temperature 
of T should be higher than B. This was found to be so, the temperatures 
of T and B being 2-25° K and 0-91° kK respectively. The way in which the 
temperatures of T and B changed when the cryostat was allowed to 
warm up is shown in fig. 5, and it will be seen that there was an initial 
fall in the temperature of T. The transfer rate calculated as before is 
shown in fig. 6 (curve (b)) together with the evaporation rate while 
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pumping (point (a)). The explanation for the lack of correlation between 
those two and for the initial fall in temperature of T is as follows. When 
the cryostat was being pumped the film flow along G, being limited by 
that which could creep through P, was fairly small, so that it evaporated 
below T. When the pump was turned off, however, the helium which 
condensed back did so as film on top of the cryostat and was available 


Fig. 4 


Experimental cryostat : constricted pumping tube. 


for further transport along G, where the rate was subcritical. Hence the 
transfer rate along G got bigger and the film reached a higher point in 
the tube, thus cooling T. The rate reached a maximum value which was 
less than the critical value, as can be seen from fig. 6, where we have 
plotted (curve (c)) the transfer rate for a 6 mm diameter clean glass tube. 
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Fig. 5 
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Warming up curves for cryostat shown in fig. 4. Curve B—temperature of 
cryostat. Curve T—temperature near top of the pumping tube. 


Fig. 6 
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Transfer rates for cryostat shown in fig. 4: (a) steady evaporation, 
(b) calculated from warming-up rates, (c) measurements of Daunt and 
Mendelssohn (1938) 
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The fact that the rate on G was subcritical during the warming-up 
period is of no fundamental importance. It could be brought up to the 
critical value either by decreasing the dead-space so that the amount of 
helium refluxing is bigger or by increasing the transport rate through P. 
In an experiment where the German silver capillary was replaced by one: 
of the same size made of copper, which happened to give a transport rate 
about twice as big, it was indeed found that, while pumping, the rate on G 
was subcritical, but when the pump was turned off the rate quickly built 
up to the critical value. 

The difference between the transfer rates on G with the pump on and 
off arose because the capillary, being an integral part of the cryostat, was 
in good thermal contact with the bulk helium inside. We would expect 
that if this thermal contact were very bad the refluxing helium would 
condense only inside the cryostat (i.e. below the constriction) and the 
transport rate along G would be always the same as that through the 
constriction. This was confirmed by an experiment where the constriction 
was thermally insulated from the top of the cryostat. The constriction 
was a glass capillary 4 mm long and 0-5.mm diameter which was 
attached to the metal parts by means of small Housekeeper seals. 
The final temperature reached while pumping was 0-8°K, with an 
evaporation corresponding to a transfer rate through the glass capillary 
of 7-8 10-5 ¢.c./em sec. The rate calculated from the heating-up process 

_was 8-2 10-5 c.c./em sec at about 1° K, in good agreement with both the 
steady evaporation and the values given by Daunt and Mendelssohn for 
glass. 

§4. Discussion 

Apart from the conclusions already drawn we may use the results of 
the above experiments in the design of low temperature cryostats. If, 
for example, the only connection to the cryostat is the pumping tube, 
then by making the pumping tube longer than a certain minimum length 
we secure no reduction of heat influx into the cryostat. It is only necessary 
to ensure that the film evaporates inside the tube, and with a German 
silver tube of 0-1 mm wall thickness the minimum length is 3 mm if its 
upper end is at 4-2° K and 70 mm if its upper end is at 20° K. 

In the case of cryostats to which there is an external heat influx (e.g. 
along additional tubes or through gas conduction) it has been generally 
assumed that this would produce an evaporation from the bulk liquid 
which is additional to the film evaporation, the two being independent 
(Blaisse, Cooke and Hull 1939). It appears to us, however, from a 
consideration of the temperature and pressure distributions along the 
pumping tube, that this may not be the complete picture, and that the 
effect of an external heat influx is rather to evaporate film than bulk 
liquid, provided of course that it does not exceed the total heat of 
evaporation of the film. 

Even if this were not so (we hope to investigate this question 
experimentally) an arrangement could be devised which in any case 
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would ensure that the external heat influx would evaporate film rather 
than bulk liquid. For example, any additional tubes to the cryostat 
could be brought into thermal contact (e.g. by a copper link) with the 
pumping tube at a short distance above the cryostat, and the heat influx 
along all the tubes will then be used to evaporate the film. Provided the 
heat influx is not too big the film can reach the level of the link so that 
at the point of contact all tubes are at the cryostat temperature and there 
can be no heat flow to the bulk liquid. In the same way heat influxes 
due to radiation and gas conduction can be minimized by surrounding 
the cryostat with a copper shield attached to a suitable point on the 
pumping tube. In brief, the cryostat should be so designed that any 
external heat influx evaporates the film and not the liquid. 


We wish to thank Professor F. E. Simon for his interest in this work 
and the Department of Scientific and Industrial Research for E. A.’s 
maintenance grant. 
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SUMMARY 


Tn this paper the problem of the thermal instability of an incompressible 
fluid sphere heated within and in equilibrium under its own gravitation is 
considered. A general disturbance is analysed into modes in terms of 
spherical harmonics of various orders, J, and the criterion for the onset 
of convection for the first fifteen modes is found both when the bounding 
surface is free and when it is rigid; and it is shown that in both cases 
the mode /=1 is the first to be excited. 


§1. INTRODUCTION 


THE manner of the onset of convection in a horizontal layer of fluid heated 
below has been the subject of experimental and theoretical investigations 
since Benard’s (1900) first discoveries on the subject. The related problem 
of thermal instability of a fluid sphere heated within has not received the 
same attention though it is of interest in a number of geophysical and 
astrophysical connections. However, a beginning in the study of this 
problem was made by Wasiutynski (1946) who derived the basic equations 
of the theory. More recently, Jeffreys and Bland (1951) have re-examined 
the same problem along somewhat different lines; in particular, they 
have attempted to solve the underlying characteristic value problem by 
a variational procedure. In this paper we shall show that Jeffreys and 
Bland have not formulated the variational procedure in its most effective 
form ; and that when this is done it is capable of solving the mathematical 
problems of the theory to a high degree of precision. While this paper 
will be limited to the case of thermal instability in a homogeneous sphere, 
the methods of this paper can be readily extended to the treatment of 
convection in spherical shells—a problem of very specific geophysical 
interest ; it will therefore be developed in that connection in a separate 
paper. 
§2. THE EQquaTIONS OF THE PROBLEM 

While the basic equations of the theory have been derived in the 
papers of Wasiutynski and Jeffreys and Bland, the perturbation in 
temperature which they use as the dependent variable does not appear 
to be the most convenient cne to use; the radial component of the 
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velocity u, (or, rather, ru,) appears to be a more suitable one and has 
the additional advantage that the boundary conditions can be more 
readily expressed in terms of it. For these reasons, it may be useful to 
develop the equations of the problem ab initio. 

Consider then a homogeneous sphere of radius R, in equilibrium under 
its own gravitation and with a uniform distribution of heat sources such 
that in the absence of conduction the temperature at each point will rise 
at a rate «. In the steady state the temperature distribution, 7'(r), 
inside the sphere will be governed by 


PA i a en 


where «x is the coefficient of thermometric conductivity. The solution of 
eqn. (1) appropriate to the problem on hand is 


T=f(R?—r?) where B=e/6x...- . 3, a eee 


In writing the solution in the form (2) we have assumed that 7’=0 at 
r=R; this clearly entails no loss of generality. 
Let the temperature distribution in the perturbed state be given by 
T=B(R*?—17)-+0. -. 2 5 
The equation governing @ is 
00 oT 
ap t Mig, KV, a Pe ene ete Se TS {4\ 
where u; denotes the component of velocity. Further in eqn. (4) (and 
in the sequel) summation over repeated indices is to be understood. 
The equation of motion is 


du; a ap av : 
Pa TP ay Mime aaah aah eNaeele ae 


where V denotes the gravitational potential and the rest of the symbols 
have their usual meanings. Following Rayleigh (1916) we shall allow 
in eqn. (5) for the variation of density only in so far as it modifies the 
external field. Thus, in eqn. (5) we shall replace p which occurs in front 
of dV/dx; by 

p=po(1—«T’), Mere ves ey ke 


where « denotes the coefficient of volume expansion and Po is the density 
at r= (where 7'=0) and regard p occurring elsewhere in eqn. (5) as a 
constant equal to py. On these assumptions eqn. (5) becomes 

du; a 0 [p oV 
ERS rhb iak alte -) pres uma Md rai ch : (7) 


where we have substituted for 7’ according to eqn. (3). With the 
variation of density due to thermal expansion allowed for in this manner, 
we shall from now on treat w,; as a solenoidal vector 


Ca 


a (), 
Ox; : 


(3) 
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When the external field is that due to the action of gravity in a 
homogeneous sphere, 


OV 4 
[vere 37 PGE ;, 5 ° ° . . . . . (9) 


where G denotes the constant of gravitation. With the substitution (9), 
eqn. (7) may be written in the form 


Ou; fa) Ow 9 
Sarr macs ge is 2 3 10) 
where 
4 
Y= 37 PGe, ‘ be Sh! be Ee) 
and 
p 1 r 
= — = V——fby(2R*r—rt), 2. . ww el 612 
oem 4 Py ) (12) 


We shali now suppose that uw; and 6 are small quantities of the first 
order and that we can ignore products and squares of them. On this 
supposition eqns. (4) and (10) reduce to 

= KVL Bue, oO a wey ce ol LS) 
and 
Ou; 
Ot 
Taking the divergence of eqn. (14) and remembering that uw; is solenoidal, 
we get 


——_ Oo yee va. ea Rai” o eae (14) 


P 06 
Coes 
Vo=y (20-42, m=) : were. re. (LD) 
Equations (13) to (15) are the basic equations of this theory. 


§3. THe BounDARY CONDITIONS 
The solution of eqns. (13) to (15) must satisfy certain boundary 
conditions. Some of these conditions will depend on whether the bounding 
surface considered is rigid or free. But in all cases we must require that 


6=0 and w,x,=ru,=0 (ona bounding spherical surface). . (16) 


Additional boundary conditions needed to make the problem determinate 
follow from the equation of continuity ; and these depend on the nature 
of the bounding surface. 
Now the equation of continuity in spherical polar coordinates (r, 9, ¢) 
is 
Oust, Lt, | Ut, cot d 1 Ouy 
ane ree G8 ip r sin & 0p 
where Os ug and u, are the components of velocity along the principal 
elements of arc dr, rd% and rsin 9dgq, respectively. On a bounding 


Order LT) 
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surface which is rigid and on which no slip occurs, U;, us and wu, must 
vanish; and when the bounding surface is spherical, it follows from 


eqn. (17) that 


= =) (on a rigid spherical surface), . . . (18) 
r 
or, equivalently (since u,—0) 
£ (ru,)= aC i2,)=0 (on a rigid spherical surface). . (19) 
r r 


On the other hand on a bounding spherical surface which is free, we 
must require (in addition to u,=0) that the viscous stresses p, and Pg 
also vanish. Now, quite generally, the expression for these stresses are 
(cf. Ramsey 1949, p. 374) 


Guig. 1h Ou 
po=pv (a5 — 24 2), 
Ou, Uy , My 
rds ea (; sindop or =) i 
Since wu, vanishes (identically) on a bounding surface (cf. eqn. (16)), the 
vanishing of p,, and p,, on such a surface requires that 


(20) 


dus peer OU ONG (on a free surface), . . (21) 


From the equation of continuity (17) it now follows that 
(5; — *) (S +2u,] = 0 (on a free surface). . . (22) 


On expanding the left-hand side of this equation and remembering that 
Uu,=0 on a bounding surface, we find ~ 

O7u,, CU noe 0” 
ae TP Gp = Ga MI= Ga 
The boundary conditions (19) and (23) agree with those derived by 
Jeffreys and Bland (1951) by a somewhat less direct line of reasoning. 


Yr 


(u,v ,)=0 (on a free surface). . (23) 


§4. THe Equations GOVERNING MARGINAL STABILITY 


In this paper we shall assume that the principle of the exchange of 
stabilities is valid ; its justification for the problem on hand will be given 
in a later paper. On this principle the equations governing marginal 
stability are those for which 0/dt=0. Hence setting 0/0t=0 in eqns. (13) 
to (15), we have 


kVil=—28u a. 


Ow 
vV7u = on. — 0x ;, . . . . . . e . (25) 


and 


: 00 
Vao=y (s+a50). wi at) +) Sone 
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We can eliminate w from these equations in the following manner : 


, Ow 0 5) 00 
V2 . 2 = oe =>=>—— 2 =v =—— — 
(y0x ,+vV?u,)=V iy Vay ae (s0-+2, =) = ee (27) 
On simplifying the outer parts of this equality and rearranging, we find 
7) 0 0 
4 if ee y 
VU = em - oe on, x,V Je A ser fe 4.) 
On operating this last equation once more by V? we find 
0 0 7] 
6 ts 2.02) 92 ‘ 
Ve oe L (= + —— an fe on, x,V ) Veo. eee C2) 
or making use of eqn. (24), we have 
= 0 0 
8y <= — Ee. Ve 
Vou; 2 (ee -|- ee. a, xV ) (w,x)). ae (30) 


It can be readily verified that the operations 2; and V? are permutable 
when applied to a solenoidal vector. Therefore, multiplying eqns. (28) 
and (30) scalarly by w,;, we ees. : 


0 7) 
4 (a ihe ON PYRG 
Vt (u,e,)= “(nize oa ae aaa Ge, r ve) é eer a(S.) 
and 
2B 0 0 an 
V8 (wan) =— "2! (wear tigen, ae) (um). . (82) 
The differential cee which occurs on the right-hand side of eqns. (31) 
and (32) is 
Q g d 2 2 pd o° 2 g ) 
ee a at var(s Aleut On V 
2 
eee g sin Ree — ere == (Say). 2. (33) 


sin 3 08 oF sin? b dg? 


Hence, eqns. (31) and (32) can be written in the forms 


VE (ujx,)= "1? ae ea Ae (34) 


2 
and V8 (ujn)=— "1 (i) ame Wire are Fe. (OD) 


It is apparent from the forms of eqns. (34) and (35) that the general _ 
solution of these equations must be expressible as a series in spherical 
harmonics. We, therefore, seek a solution of the form 

6=O(r)S (9, p) and uw=ru,=Wi(r)S/(9, ¢), . . - (36) 
where S, (9, g) is a general spherical harmonic of order / and @ (r) and 
W (r) are functions of r only. Since 

L?8(9, p)=UL+1)S(9; &) 
GO le Saas) (37) 
ppd Das ai Carrer Pp? 
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eqns. (34) and (35), for the assumed forms of @ and u,xv,, reduce to 


@ 22d Utley y | 

[atc + | w-u+n Fe oo ae 
and 

@ 2d act 2By 

Ee —— | w= meats A 


In the further discussion a eqns. (38) and (39) it is convenient to 
measure r in units of the radius R of the sphere. With this choice of unit, 
eqns. (38) and (39) become 


d? Dal 
ES rar 


eae 1) 
d d | w= 


26 ee 


and 


det rdr 
where (cf. eqns. (2) and (11)) 


Wi+1)C;W, . . . (41) 


2 4 
oa Et Re SO 
KV 9K" 

A subscript J has been attached to C to emphasize the fact that the 
characteristic values of this number will depend on the order of the 

spherical harmonic considered. 

According to eqns. (16), (19), (23) and (40), the boundary conditions 
with respect to which we must seek a solution of eqn. (41) are that at r=1 
2d ey 

Me ~ Lie rdr | ee 
LW |) RS on 
and either os OPS rae =0 (r=1). 
Further, at r=0 none of the physical quantities must have a singularity ; 
and it can be directly verified from eqn. (41) that this condition requires 


W to behave like r’ asr—-0. . ... 5 Ceeentaas 


A solution of eqn. (41) which satisfies the boundary conditions at 
r=0 and r=1 and is not zero everywhere is possible only when C, takes 
one of a sequence of determinate characteristic values. In the preeeee 


connection we are interested only in the lowest of these characteristic 
values. 


§5. A VARIATIONAL PROCEDURE FOR DrTERMINING THE LOWEST 
CHARACTERISTIC NUMBERS OC, 


Letting 
CMW eae TE a) li defena 
et sae | Yale) e+ | w, 


Pas 2d. Yt end ae ie 
ot apart W=2|5("z) —u+1) ] 4, 
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we can rewrite the differential equation governing W in the form 


l[d/,aF 
=| 5 (" 7) —1(I+ nF | =—iI+1)C;jW. .-. . (46) 
The boundary conditions (43) require that 
F=W=0 and either al or ioe at r=1, Sere, (47) 
dr dr 


Multiply eqn. (46) by r2F and integrate over the range of r. The right- 
hand side of the equation gives 


ig /-. dP oe 
[Pala )er—tety |e red tee ma (18) 


By integrating by parts the first of the two integrals in (48) and 
remembering that F vanishes at r=1 and is continuous at r=0, we obtain : 


= fie (G) tp} a. BE, eter (15) 


Turning next to the left-hand side of eqn. (46), we have to consider 
(cf. eqns. (4 a 


“AWE dr= we Bate dr—l(l-+-1) “WG dr, severest (50) 
dr 


By ee by parts Ae first of the two ae on the left-hand side 
of (50) we get (cf. eqns. (44) and (47)) 


dG fd WN\ 2 
=| ae (: 7) pe Weesiess tise (OL) 
a further integration by parts yields 
dW eride (a dW 
—|G—- = nd Ae Rha 2 
(4 7), +). ¢%(" or az) (52) 


Ifr=1 is a rigid boundary (dW/dr),_,=0 and the integrated part vanishes. 
On the other hand, if r=1 is a free boundary, (dW/dr),_,40 and G 
will have the value (cf. eqns. (45)) 2(dW/dr),_, since W,_,;=0 and 
(d?W/dr?),_,=9. Hence in this latter case the integrated part in (52) 
will have the value —2 (dW/dr)?,_, ; we may clearly take this as the value 
of the integrated part in both cases. Thus, combining (50) and (52) we 
can write (cf. eqns. (45)) 


(7) _ +i ee ee on oe ae (53) 
Hence the result of multiplying eqn. (46) by 7?/ and integrating is 
1: 
{ (2(dF[dr+-U41)F2}dr 


Ul+-1)0 = 2 —______—__ = 7 (say jaes wono4) 
—2(dW/ar)_.+ | 2G? dr 
0 


[i ed 
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Consider now the effect on C, of a variation 5W in W compatible with 
the boundary conditions on W. To the first order, we have 


1 
U(+-1)80 = 7 (S1,—UI+1)0fIn},  . - + + (58) 
2 
where SJ, and 8, denote the corresponding variations in J, and J,: 
sh=2 | {eh oF+UL+ EOE | dr eae 


and 


1 . 
Reet (eee, +2] P63¢ dr. ae 
2 adr dr r=1 0 


The foregoing expressions for 52, and 5I, can be reduced by one or more 
integrations by parts as follows : 


Sti —2[ {5 . G =) —u+yF bor dr 


wees 2| {5 Ui+1) FSF dr es 


(ta) 
yee (x W) +2 ji r(? sr) + n)W ee ar 
a 6 6) od aa 


and 


dw d 
=-4(7 5° 
+2) ws. Ge. Wsrp0)a 
dW d dw i 
ae -4(F, Sp) 4a(Maa), +2] 0 oF ar. ES 


The two terms besides the integral in this last expression for 62, vanish 
if r=1 is a rigid boundary (since in this case dW/dr=0). On the other 
hand, if r=1 is a free boundary then 


(CG Hetsa jw] = (79 W) = a 


and 
~ (6 aeW) +2(G, 9) calms Dra ne 


=1 


Hence in both cases 
1 
51,=2 | SWF dr... 2 8 Sie 
0 
Now combining eqns. (55), (58) and (62), we have 
2 fee Kisdet soe Ao 
UI+-1)8C j= — abe laa G =) F+UI+1)C, wh oF dr. (63) 


Hence to the first order, 53C,=0 for all a arbitrary variations 5F. It is 
evident that the converse of this proposition is also true. Further, it 
follows from (63) that the true solution of the problem gives an extremal 
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value for C;. This last fact enables us to formulate the following variational 
procedure of solving eqn. (41) and satisfying the boundary conditions of 
the problem. 

Assume for F an expression involving one or more parameters A, 
which vanishes at r=1 and behaves like ras r+0. With the chosen foun 
of F determine W as the solution of the differential equation 


a 2d ae 
late AS eee 


which satisfies the conditions at r=0 and r=1. Then evaluate C, 
according to formula (54) and minimize it with respect to the parameters 
A;. In this way, we shall obtain the ‘ best.’ value of OC, for the chosen 
form of F. 

The variational procedure as formulated above differs from that of 
Jeffreys and Bland (1951) in that they do not seek W (rather @ in their 
connection) as a solution of a fourth order differential equation. It is 
necessary to do this since 8C,;=0 only for arbitrary variations 6F i.e., 
only for such variations 6W which satisfy the differential equation 

a@ 2d il+1)7/? 
a ve | SW rae eens.) (85) 


where dF is arbitrary. 


§6. Tor DETERMINATION OF C’; FOR VARIOUS VALUES OF / FOR THE CASES 
(i) WHEN THE BounDING SURFACE IS FREE AND (11) WHEN IT IS 
RiGIp 


The form of the equation to be solved for W suggests that we assume 
for F the trial function 
uf 
i= Woe Aj J pr) 9(%57); . . . . . . (66) 
where J 1/2 denotes the Bessel function of order /-++ 3, the «,’s (j=1, 2, ...) 
are its zeros and the A,’s are the variational parameters. With this choice, 
F vanishes at r=1 and behaves, as required, like r’ as r—0. Also, it may 
be recalled here that the functions J;,1/2 («,;”) satisfy the orthogonality 
relations 


1 
for J pala, (457) J p4112 (a) dr=t 85x [J 1/2 CON yao (67) 


where primes denote differentiation with respect to the argument and 

8,;, is the usual Kronecker symbol. 
With F given by (66) the equation to be solved for W is 

a@ 2d iil+1) 

aay tia 2. | W= GPA Tusale) - + (68) 

Since, 
a 2d Ub+1) T rpolar) 
| aT dr re Vt Be 
4U2 


= Fuse ea, (69) 
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a particular integral of eqn. (68) is 
Liens 
a og Teal als) . ° . . . . ° (70) 
Adding to (70) the complementary function Br’+Cr'” which has no 
singularity at r—0, we obtain the general solution : 

ligeeay 

= pig ane ee . s* Lao a 
3 

(In (71), B and C are constants of integration.) The condition W=0 


at r=1 requires 
B=—C0, °° 2 3 os a 
Leeds : 
and we have V= Reo g Auta) CE eaeaaa 3 dee an 
A 

The constant B is determined by the remaining condition at r=1 namely, 
that here either dW/dr or d?W/dr? must vanish (depending on whether the 
bounding surface at y= 1 is rigid or free). We find that 

B= 392(A 5/037)" r41/ 0(%5)s «tenes ee 
where G=2 for a rigid boundary at r=1 


<a 4 =e . ° (75) 
ieee for a free boundary at r=1. 

Turning next to the evaluation of C, according to formula (54), we find 
that for F and W given by eqns. (66) and (73) (cf. eqns. (48) to (53), (67) 


and (68)) : 
: adF\2 
| {r (=) ++)P dr 
0 r 
: CF 2dF Ul+1)F 
came |e ei, ES NE Da 
e [ert gate dr 7 \ ar 
1 
= J, PEAS peri aleir) 2p pope S r1/2(a gr) ar 
= §2A Por i aiel%s)l | (76) 
and 
1 1 
| rG? dr—2 (Fe). ce | rWE dr 
0 T J r=1 0 


, , 
=32(A P04)" 1.10053) P+ B2;A, iL (nee) eeastare) dr. 

(77)* 

The remaining integral on the right-hand side of (77) can also be evaluated 


if proper use is made of the various recurrence relations satisfied by the 
Bessel functions. We find : 


1 
Ns ee d= F(A Poe IT 14 1/ (05) P+ 2BE,(A j00;2)T 145) 9(05)- » + (78) 
* The transformations used in going from (48) to (49) and similarly from 


(50) to (53) can be used in the reverse fashion since these depend (a 
ee art fi 
definitions) only on the boundary conditions of the problem, : east 
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The resulting expression for C, is, therefore, 
id+-1) C,= 
NA Pa? lI’ r1/9(a,) 27 
Z(A Pos I" 1.1/2(%j) P92 (Aj[o? I 4 1)9(4))>y(A wl ie EF sy 2(% x) 
(79) 


where we have substituted for B according to (74). By minimizing this 
last expression with respect to the A,’s we shall obtain the ‘ best ’ value for 
C’, for the chosen form of F. 

The simplest trial function for F of the chosen form ig 


rt 
Fr= aJr) eeualeat); . . . . ° . . (80) 


Table 1. The Characteristic Numbers C; in Case the Bounding 
Surface if Free 


First Second 
Approximation A Approximation 


SN 


1 
2 
3 
4 
5 
6 
v 
8 
9 


3094-0 

5227-4 

8778-6 

1-3986 x 104 
2-1209 x 104 
3-0853 x 104 
4-3360 x 104 
5-9203 x 104 
7-8887 x 104 
1-0295 x 10° 
1-3196 x 10° 


11-6650 x 10° 


2-0721 x 10° 
25475 x 105 
3-0978 x 10° 


—0-022022 
—0-019119 
—0-016639 
—0-014618 
—0-:012977 
—0-011630 
—0-:010511 
—0-0095692 
—0-0087674 
—0-0080806 
—0-0074825 
—0-0069617 
—0-0065003 
—0-0060922 
—0-0057276 


3091-4 

5224-4 

8775-0 

1-3982 x 104 
21204 x 104 
3-0848 « 104 
4-3354 x 104 
5-9196 x 104 
7-8880 « 104 
1-0294 x 10° 
1-3195 x 10° 
1-6649 x 10° 
2-0720 x 10° 
25474 x 10° 
3-0977 x 10° 


where «, is the first zero of Jj,1;.. For this choice of F there is no 
variational parameter with respect to which we have to minimize and 
eqn. (79) directly gives 
a7” («+) 
U-+1)0 = 1 ees, es + (81) 

os Oy re yja(%1) +9 145/2(%1) 
where it may be recalled that q is given by eqns. (75). Values of C; 
obtained with the aid of formula (81) are listed in tables 1 and 2. 

The values given by (81) can be improved by including a second term in 


F. Thus, with the assumption 


1 
P= FS valet) + Ad reas aleat) }; vo eee (82) 
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where «, denotes the second zero of J ;,4/2 and Aisa variational parameter, 
eqn. (79), after some rearrangement, gives 


Cp (I apalee P+ Aol)’ syrralead?} 
tee U1) 141/ 2\%1 2/41 1+1/ 2\%2 


X [I eas a(%s) {1 Fray 2(%1) 9 J v.5)2(%1) } +29 A (x /%2) Jr, 5/0(01) 
J" t1)2(%2) 
+ A? (x ;/a2)° J’ py 1(%o) {%2 J" 1/242) +9 Tizsi2(%2)}P 2. - + (83) 
The values of C, obtained after minimizing the foregoing expression with 
respect to A are also listed in tables 1 and 2 together with the values of A 
which give the minimum values. 


Table 2. The Characteristic Numbers C, in Case the Bounding 
Surface is Rigid 


First Second 
Approximation A Approximation 


~~ 


8154-0 0-088804 8047-1 

1-0559 x 104 0-098580 1-0403 x 104 
1-5368 x 104 0-10409 1-5132 x 104 
2-2352 x 104 0-10729 2-2006 x 104 
3:1801 x 104 0-10914 3-1315 x 104 
4-4117< 104 0-11015 4-3459 x 104 
5:9759 x 104 0-11062 58892 x 104 
7-9221 x 104 0-11074 78108 x 104 
1-0304 x 10° 0-11061 1-0164 x 10° 
1-3176 x 10° 0-11032 1-3004 x 10° 
1-6600 x 10° 0-10991 1-6390 x 10° 
2-0637 x 105 0-10942 2-0384 x 10° 
2-5357 x 105 0:10888 2-5058 x 10° 
30814 « 10° 0-10830 3:0462 x 10° 
3°7094 x 105 0-10770 3-6685 x 10° 


1 
2 
3 
+ 
5 
6 
7 
8 
9 


From a comparison of the results obtained in the first and the second 
approximations, it would appear that the second approximation gives 
values which are probably correct to one part in 104, As a confirmation of 
the precision reached by formula (83) we may note that according to Jeffreys 
and Bland (1951) the value of C, for the case of a free bounding surface 
obtained by an ‘ exact’ solution of the characteristic value problem is 
3091+-2 (+2 arising from rounding errors); and this value should be 
contrasted with the value 3091-4 given by eqn. (83). 

Finally, we may compare the results obtained by the methods of this 
paper with those derived by Jeffreys and Bland (1951) from their formul- 
ation of the variational principle: The entries 3094-0, 5227-4, 8778-6, 
8154-0 and 10559 in tables 1 and 2 should be compared with their values 
3202, 5640, 9780, 8720 and 12600, respectively. The large discrepancy 
between these values is a reflection of the ineffectiveness of the variational 
procedure as used by J effreys and Bland. 
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§7. ConcLuDING REMARKS 


As we have already stated in §1, the methods of this paper (in particular 
the variational procedure described in §§5 and 6) can be extended, quite 
readily, to determine the criterion for the onset of thermal instability in 
spherical shells. This latter problem has a direct bearing on the question 
of convection in the earth’s mantle and on the interpretation of the 
feature of the earth’s topography to which Vening Meinesz (1951) hag 
recently drawn attention. For, as Vening Meinesz has pointed out (in 
this connection, see also Jeffreys 1952) in a spherical shell, we should 
expect that the mode of disturbance which will first become unstable will 
not be the one with J=1 (as in the case of a complete sphere) but rather 
one with a higher value of 1. The results of calculations which are being 
undertaken to establish these features of convection in spherical shells in 
a quantitative fashion will be published in the near future. It might be 
added here that the writer’s interest in these problems arose from an 
attempt to extend his recent discussion (Chandrasekhar 1952) of the 
magnetic inhibition of convection in a horizontal layer of fluid heated 
below to a spherical geometry with a view to applications to the earth’s 
core (cf. Bullard 1948). But it was felt that before attempting to solve 
the problem of convection in a fluid sphere heated within and in the 
presence of a magnetic field, one must first solve the same problem in the 
absence of a magnetic field. The writer hopes to return to his original 
problem in the near future. A further problem which the writer also hopes 
to consider soon is the extension of the analysis of this paper to a rotating 
sphere. 


In conclusion, I should like to express my thanks to Miss Donna 
Elbert for valuable assistance with the numerical work. 
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By B. H. FLowErs* 
Department of Mathematical Physics, The University, Birminghamy 


[Received October 15, 1952] 


ABSTRACT 


Magnetic moments and quadrupole moments of nuclei are treated in 
a model which considers all particles outside the core of completed levels 
on an equal footing. It is shown that in this way better agreement 
with experimental magnetic moments is obtained than from the simpler 
Schmidt model. The manner in which electric quadrupole moments 
arise in odd-neutron nuclei without the assumption of a deformable core 
is also discussed. 


In considering the magnetic moments of odd-mass nuclei it is customary 
to follow the approach of Schmidt and to assume that the addition of a 
pair of like particles occupying the same orbit has no effect upon the 
properties of the nucleus. It is indeed true that the symmetry properties 
of the low-lying states, particularly the spins, are not changed by such 
an addition; but it is the purpose of this note to point out that 
for properties like the magnetic moment and the quadrupole moment, 
it is essential to consider all particles outside the completely occupied 
levels, whether these particles pair off or not. 

In a recent paper Mizushima and Umezawa (1952) have shown that 
for a number of light nuclei better agreement between measured and 
calculated values of the magnetic moment can be obtained by taking 
into account the interaction between neutrons and protons in jj-coupling. 
They assume that each particle outside the core of closed shells possesses 
a magnetic moment equal to the usual one-particle value as used in the 
Schmidt diagram (and to which we shall refer simply as Schmidt values) ; 
but they combine all such particles, both neutrons and protons, into 
states classified by means of the isotopic spin 7' and the resultant angular 
momentum J. Energy calculations indicate (Flowers and Edmunds 
1952) that for light nuclei, and for forces of reasonably short range, the 
lowest states are usually those having the smallest 7’ consistent with the 
known M,—}(N—P), where N and P are the numbers of neutrons and 
protons outside the core. If the configuration is known, and the resultant 
angular momentum, it is then possible to calculate what the resultant 
magnetic moment should be. In general, as found by Mizushima and 
Umezawa, the values so obtained differ from the simple Schmidt limits. 


ae en mi MTS I< Det ee 
2 Now at the Atomic Energy Research Establishment, Harwell. 
| Communicated by Professor R. E. Peierls, F.R.S. 


On Moments of Odd-mass Nuclei in jj-Coupling 1331 


Making use of the tables of fractional parentage coefficients calculated 
by Edmonds and Flowers (1952), which take into account also the new 
classification of states in jj-coupling according to their behaviour under 
symplectic transformation, the present author has been able to show 
that for three particles in the configuration j?, in the usual ground state 
with T=}, J=j, and seniority 1, which arises from attractive forces of 
reasonably short range, the magnetic moment is given by 


= (67+ 6)-1{(2j—1)u,+(4j-+5)u,} for an odd proton 

and > (1) 
L=(6)-+6)~-1{(4j +5) u_+(27—1)u,} for an odd neutron 

where jin, tp are the Schmidt moments of neutron and proton, respectively, 


appropriate to the given values of / and j. Similar formulae exist for 
more than 3 particles, but for these we shall state only numerical results. 


Table 1. Magnetic Moments of the Light Nuclei of Odd Mass 


Nucleus Configuration ree Schmidt Individual Observed 
ali, (Ps/2)* 4 3-79 3-04 3-26 
iBes (Ps/ 2)? 4 —1-91 —1-16 2 
Be; (Pars) 4 —1-91 —1-16 —1-18 
(Ps/2)7* 3 3°79 3°79 9-69 
11 
2 (Pya)"Puale? F 3-79 3-04 
130, (Py 2) 1 0-64 0-64 0-70 
Ns (Pye)? 3 —0-26 —0-26 —0-28 
140, (ds) 2)2 4 —1-91 = 1-9) —1-89 
cae (d,, 2)*28;/o 4 2°79 2-79 2-63 
ioNe (d52)® 1 1:01 Syl <0 
iiNaie (dsy2)-® 3 (2-87) 2-0 2-22 
i2Mgis (ds) 2) 4 —1-91 yi ee —0-96 
(ds; 2)~* 4 4-79 Fa 
Ale (dsj2)*(8y/202 4-79 3-52 a 
i7Clis (dg/2)8 4 0-13 0-26 0-82 
70l,, (dg/2)-8 8 0-13 0-13 0-68 
$9Sco4 (fo) 2)5 3 5-79 5-10 4-76 
33Mns35 (f/2)p °(Ps/2)n” cs (4:13) 36 3°47 


The results are shown in table 1 for odd mass nuclei of mass number 
A<55. In all cases the experimental moments are taken from the recent 
compilation by Klinkenberg (1952). The following features are worthy 
of notice. The nuclei #0, N, 170 consisting of fully occupied orbits 
plus or minus one particle, have moments in excellent agreement with 
the Schmidt values. On the other hand, 1B, and ?7Al, which are usually 
thought of as closed shell +1 nuclei, are not in such good agreement. 
In the case of 11B one can improve matters by assuming that two neutrons 
are raised from the ps. orbit into the nearby py, orbit, thus implying 
(not very surprisingly) that for such a light nucleus jj-coupling does not 
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strictly apply. Concerning °7Al, it is known that the single-particle 
levels d,/. and 2s,/, lie close together, and it is seen that the configuration 
(d5/2)~3(Sy/2)2 now proposed for this nucleus yields a magnetic moment 
in good accord with experiment. It would clearly be of interest to 
obtain measured values for 0 and ™N in this connection. The nuclei 
23Na and 55Mn are two of the three spin exceptions to the jj-coupling 
model, having J=j—1 rather than J=j in the ground state. For these 
nuclei there properly exist no Schmidt limits, but the values shown in 
brackets in the fourth column are those derived from considering the 
proton configuration only, namely 
j—1 
[pi j Schmidt » 

The corresponding values in the fifth column have not been calculated 
completely because the necessary fractional parentage coefficients are 
not available, and are thus only rough estimates. For *°Mg, Mizushima 
and Umezawa were able to state only that the magnetic moment should 
lie between —1-06 and 2-76. Our unique value of —0-64 is that appro- 
priate to the state of seniority 1. The above authors calculate —0-48 
for the T7=4 state of *Cl and 0-80 for the 7'=3/2 state, and conclude 
from the measured moment that this nucleus has 7=3/2 in the ground 
state. Their values, however, appear to be in error in this case; it is 
clear from formula (1) that since 1, and pz, are both positive for a d3/. 
orbit, the magnetic moment of the 7’'=} state is necessarily positive, 
as indeed experiment shows. 

In all cases, values of the magnetic moment based upon the individual 
particle model are in much better accord with experiment than those 
based upon the naive Schmidt model. The importance of this result 
lies in the fact that, under favourable conditions, it is now possible to 
deduce from the measured moment not only the orbit of the odd particle, 
but the whole configuration giving rise to the moment. Moreover, it 
may not be necessary to appeal to quantized motions of the core, of the 
kind proposed by Aage Bohr (1952), or to large exchange moments as 
considered recently by Russek and Spruch (1952), to explain the very 
large deviations from the Schmidt limits. For the heavier nuclei, it is 
expected that a similar result will hold, and that much of the deviation 
from the Schmidt values can be explained in terms of the coupling between 
neutron and proton configurations. The case of 2°Bi, however, remains 
unsolved. 

Regarding the electric quadrupole moments of odd mass nuclei, 
Rosenfeld (1951) has pointed out three striking features: (i) the quadru- 
pole moment of a nucleus containing one hole (apart from filled levels) is 
approximately equal, but of opposite sign, to that of the corresponding 
nucleus having one particle in this level; (ii) the quadrupole moment of 
an odd-neutron nucleus is roughly the same as that of a nearby odd-proton 
nucleus ; and (iii) heavy nuclei with 100<A<200 have with very few 
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exceptions quadrupole moments about 10 times larger than one would 
expect on the basis of a single particle model. 

The first feature follows immediately from the theory of holes, and we 
shall concern ourselves only with (ii) and (iii). We have ee the 
quadrupole moments of configurations 7° in the states T—1, J=j, and 
aa 1, taking as individual particle values Q,=1 and Won We 

n 


Q(1p-+2n)=(4j+5)/(6j+6) and Q(In+2p)=(2j+7)/(Gj+6). . (2) 
These expressions give Q S 0-7 in both cases for small values of j. Thus 
the quadrupole moment “induced” upon the odd neutron by two 
protons is about equal to that of an odd proton in the presence of two 
neutrons, as observed by Rosenfeld. 

The quadrupole moments of odd-neutron nuclei are therefore found 
to arise quite naturally from the coupling between neutrons and protons. 
It follows that a nucleus consisting of an odd neutron outside a double 
closed shell should have quadrupole moment zero (whatever the spin of 
the nucleus) for in this case there are no protons with which the neutron 
can readily interact. The only experimental evidence concerns 170 which 
has Q=— 0-005 x 10-4 cm?, or about 10-20 times less than those of nearby 
nuclei. It would be interesting to have experimental values for °1Zr, 
and 79°Pb, 

Further, we may enquire what effect the addition of two neutrons to 
an odd-proton nucleus may have upon the quadrupole moment. 
According to (2), the moment should be reduced by about 0-7. In 
table 2 are presented all such ratios known for masses less than 100: 
all are found to be in the right direction and to be surprisingly consistent 
with the theoretical factor bearing in mind that in some cases the neutrons 
and protons occupy different orbits. 


Table 2. Quadrupole Moments of Odd Mass Isotopes 
Differing by Two Neutrons 


Nucleus ee ee Added neutron pair @x 1074 cm? Ratio 
wl —0-07894 

39CU —0:13 

an \ (Ps) 2) * (Ps/2)? or (£5) 2)” —0-12 0-92 
31Ga +0-2318 

ek evar? GasPortea? 4 foaar fo 
Bp +0:26 : 

Br \ (Ps/ 2) * (P12)? OF (9/2)? 40-21 0:8 


Finally, concerning the third feature emphasized by Rosenfeld, it is 
clear that odd-neutron nuclei should have large quadrupole moments 
only when the protons are entering levels of high angular momentum. 
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There is only one odd-neutron nucleus known to have an excessively 
large Q: this is 1Yb with a neutron configuration (f,/.)? according to 
Klinkenberg, and the proton configuration of this nucleus appears to 
contain either 6 or 4 holes in the h,,/. orbit. | We estimate very roughly 
that this could cause an increase in the Q for a single f;). particle by a 
factor of as much as 5, which is at least of the right order of magnitude. 
Similarly, the large quadrupole moments observed for odd-proton nuclei 
in the mass region A & 150 all seem to be connected with the pairwise 
filling of the h,,/. orbits by protons. One would expect similar results 
for the transuranic nuclei where the i,3;. proton levels may be filling in 
pairs. 


This work will appear in greater detail elsewhere, together with 
applications of the methods used here to the calculation of transition 
probabilities for electromagnetic radiation. It is a pleasure to thank 
Professor R. E. Peierls, F.R.S., for the interest he has shown in the 
work. The paper was written while the author was on leave of absence 
from the Atomic Energy Research Establishment, Harwell, and the 
Director’s permission to publish it is gratefully acknowledged. 
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A Counter-Hodoscope Study of Associated Penetrating Particles Underground 


By B. Lrontic and A. W. WourenpaLE 
The Physical Laboratory, University of Manchester* 


[Received September 29, 1952] 


THE production of pairs of associated penetrating particles (APP). has 
been examined using a counter hodoscope in a cave at Stockport under 
26 metres water equivalent of rock. 

The apparatus, shown in fig. 1, consisted of 18 counters arranged in 
two trays AB and C. The master pulse required to trigger the hodoscope 
was produced when at least one counter from each of the groups A, B and 
C, was discharged. Each counter then operated a corresponding neon 
lamp. The counters in the lower tray were separated by 1 cm lead slats 
in order to reduce the effects of knock-on electrons. 

Events in which two counters in group C and one counter in each of 
groups A and B were: discharged, were examined. Events have been 
observed which we interpret as pairs of APP originating in the roof and 
in 2’. Accidental coincidences, due to the passage of two unassociated 
particles through the apparatus, were negligible since the resolving time 
of the coincidence circuit was very short (1 psec). The events are classified 
as follows: adjacent pairs, in which adjacent counters were discharged 
in C and AB; separated pairs in which the counters discharged were not 
adjacent and the trajectories show that the particles diverged from a 
point in the roof or in 2’; parallel pairs, where the trajectories are nearly 
parallel showing that the particles came from the roof. Our results 
are given in the table. The number of adjacent APP’s has been corrected 
for the effects of successive knock-ons in the two trays, by using the 
knock-on frequencies observed for the individual trays. 

The relatively small number of separated pairs compared with that of 
adjacent pairs is to be expected from the known approximate angular 
distribution of the secondary about the direction of the primary, as 
observed in cloud chamber experiments (Braddick et al. 1951) and the 
geometrical bias of the split-tray arrangement. 

The absorption of particles from the roof, measured by putting 20 
and comparing the rates with t=16 cm and t=6 cm, leads io a value for 
the interaction length of the secondary particles 


L=14-+-4 em Pb. 
SS SS ee ee 
* Communicated by Professor P, M. 8, Blackett, F.R.S, 
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This value agrees with that calculated by us from the results given by 


Amaldi et al. (1952) which lead to 
L=10+4 cm Pb. 


The conclusion from this experiment is that the secondary is a nuclear 
interacting particle—possibly a 7-meson. 


Fig. 1 


15 CM Pb 
x 


“GEQnSeoeSS——-¢ 


Appapillai, Mailvaganam and Wolfendale, in work not yet published 
ve te ree ue: evidence in support of this conclusion 
e lack of sufficiently ‘accurate knowledge of istribati 
ge of the angular distribut 
of the secondary about the direction of the primary makes it im oats 
to derive a cross-section for production, : : 
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Number of | Corrected 


Oe t Time adjacent number of | Separated | Parallel 
(cm Pb) | (cm Pb) | (hours) coinci- adjacent 


dences pairs 


pairs pairs 


15 16 448° 221 82 
15 6 253 274 196 
Oar 16 383 219 101 
0 6 330 321 219 


Fig. 2 


°% of single meson flux 


oe 2 4 6 


No, of intervening undischarged counters, 


1338 Correspondence 


Events have also been recorded in which one counter in C and two non- 
adjacent counters in AB were discharged. Similar events have been 
observed by George and Trent (1951) who ascribe them to the production 
of APP’s. An alternative explanation in terms of a simpler phenomenon 
is given here. 

Figure 2 shows the distribution of the counting rate of these 3-fold 
coincidences as a function of the number of intervening undischarged 
counters in the AB tray (simple geometrical corrections have been applied 
to the rates). It will be seen that the rates for t=6cm Pb and 
t—16 cm Pb are virtually the same. The mechanism suggested for these 
events is that a penetrating meson produces near the bottom of ¢ a 
knock-on shower from which low energy photons of energy in the region 
of 3 Mev penetrate the absorber and a counter is discharged by a locally 
produced Compton-recoil electron (c.f. Greisen (1949)). The number of 
photons crossing any counter will decrease with increasing distance from 
the parent meson and so the rate of events will decrease with increasing 
number of undischarged counters. An approximate theoretical analysis 
has been carried out which gives rough numerical agreement with the 
experimental result for both thicknesses of absorber. 

A contribution to these events by the production of APP’s in ¢ is to 
be expected but this is small due to the relatively small number of 
secondaries emitted at large angles. 
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The Existence of a Macroscopic Shear in the Transformation in Cobalt 
By T. R. ANANTHARAMAN and J. W. CHRISTIAN 


The Inorganic Chemistry Laboratory, Oxford* 


[Received September 18, 1952] 
A REGION of f.c.c. crystal may be transformed into a single crystal having 
a h.e.p. structure by the movement of an imperfect dislocation, of type 
- | 
s (112), through every alternate plane of a series of {111} planes. 


The martensitic transformations in cobalt and in iron-manganese alloys 
probably occur in this way, since half-dislocations are present in a stable 
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f.c.c. lattice at the boundaries of extended dislocations (stacking faults) 
in {111} planes. There are six possible (110) slip directions in a given 


set of {111} planes, and the corresponding 5 (110) lattice dislocations 


will dissociate to give three different = (112) dislocations in these 


planes. An identical atomic arrangement is produced by the motion of 
any one of these imperfect dislocations across a glide plane. Transformation 
might occur by the independent extension of stacking fault nuclei, assum- 
ing there are sufficient dislocations to give ~ 1 region of stacking fault 
on each atomic plane, or by the propagation of a half dislocation down a 
series of {111} planes. Bilby (1951) pointed out that the existence or 
non-existence of a macroscopic shear should be a valuable indication of 
the mode of reaction Independent extension of nuclei should produce 
random averaging over the three shear directions, giving zero macroscopic 
shear. If, however, the hexagonal crystal is produced from a single 
stacking fault nucleus, the Burger’s vector of the transformation dislo- 
‘cation would be identical in a series of planes, and the formation of the 
crystal would involve a macroscopic shear of magnitude tan~1 (/2/4). 
This corresponds to the original proposal for transformation (Burgers 
1934), and implies the production of plate-shaped hexagonal crystals 
with a definite habit plane, as in other martensitic transformations. 
Barrett and Hess (1952) have reported that such plates are visible in 
thermally polished cobalt, after cooling to room temperature. 

We have observed the electropolished surface of a polycrystalline 
specimen of pure cobalt under the microscope, while heating in hydrogen 
through the transformation range. At about 415° c thin plates began to 
appear, the specimen surface being tilted about the line of intersection 
of the plates. The plates were similar in appearance to Neumann bands 
in iron, and although very narrow they could occasionally be seen to 
have a lenticular shape; their number increased until transformation 
was complete. Ina single grain there were at most four different directions 
of traces in the surface ; traces were often observed parallel to the sides 
of {111} cubic annealing twins. On cooling to room temperature the 
traces did not disappear completely, presumably because transformation 
in a polycrystalline specimen involves some non-reversible plastic de- 
formation. The appearance of a typical microsection after heating and 
cooling is shown in fig. 1. 

To identify the habit plane, both stereographic analysis and x-ray 
examination were used. In fig. 2 (a), a stereographic solution for the 
grain shown in fig. I is given ; the four trace directions are consistent 
with the habit plane being {111}. The same grain (~ 2 mm in diameter) 
was photographed in an oscillation camera using a narrow X-ray beam. 
By setting each trace vertical in turn, two hexagonal (00-1) poles were 
located exactly and the other two within a range of angles (because of 
‘the specimen thickness which limited the range of oscillation). The 
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results are shown in fig. 2 (b), and the (00-1) poles fall on the {111} poles © 


of fig. 2 (a) to within 1°. It follows that a single cubic grain forms four 
sets of plate-like hexagonal crystals, the habit plane being ails, 

The observations appear to prove the existence of a macroscopic shear, 
and thus suggest that independent spreading of stacking fault nuclei 
cannot account for the transformation. This conclusion is supported by 
unpublished results we have obtained for the variation with grain size 
of the proportion of stacking faults in hexagonal cobalt. The mechanism 


of transformation is still undetermined : Bilby and Cottrell (1951) showed 


that internal climbing of an 5 (112) dislocation is impossible in a f.c.c. 


lattice, and the calculations of Leibfried (1950) and Nabarro (1951) 
show that the reflexion mechanism (Christian 1951) is very improbable. 


Fig. 1 


Surface of electropolished cobalt after heating and cooling 
through the transformation. x 900. 


An alternative possibility, suggested to us by Dr. N. Thompson, is 
that the dislocation movements would be random in a strain-free material 
but have locally preferred directions in our experiments, because of varying 
internal shear stresses. We do not think this very probable, since in 
many regions of our photographs three intersecting sets of plates are 
visible, but it is obviously important to establish whether there is a 
macroscopic shear when a single crystal transforms. The successful 
production of single crystals of hexagonal cobalt (e.g. Myers and Suck- 
smith 1951) implies that the dislocation movements can be confined to 
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one set of {111} planes, but there is no evidence whether or not the 
relative movements of these planes occur randomly in the possible (112) 
directions. 


Fig. 2 


os 


(a) Stereographic projection of area shown in fig. 1. Lines are normal to trace 
directions, open circles are {111} poles of standard projection, and filled 
circles are {111} poles after rotation through 36°. 


iti illati . Two of the poles 
») Position of (00-1) poles from oscillation photographs 
sy enld ates be ee within the angular range represented by the 
thickened lines. 


4X2 
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In some martensitic transformations wide plate-like crystals are 
formed, and the transformation is complete in a small temperature 
interval. It seems probable that in these cases either or both of the 
following conditions must be satisfied: (1) the transformation shear is 
small, so that deformation of the matrix is not unduly large; (2) the 
transformation can proceed by reversing the shear, allowing the strain 
in the matrix to be periodically relaxed. Neither condition is fulfilled 
in cobalt, so that the martensitic plates are thin, and complete trans- 
formation is very difficult. In fine-grained cobalt, it seems virtually 
impossible to transform more than ~ 70%, of the metal, except by cold- 
working. There areso many possible explanations for the faulted structure 
of hexagonal cobalt that we shall not discuss this here. 


We wish to thank Dr. B. A. Bilby for helpful discussion, and Dr. W. 
Hume-Rothery, F.R.S. for his interest. 
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Note on Thermal Radiation at Low Temperatures 


By D. Bru 
Clarendon Laboratory, Oxford* 


[Received October 3, 1952] 


A FUNDAMENTAL result of statistical mechanics and quantum theory is 
the conclusion that, for any system at sufficiently low temperatures, the 
usual thermodynamic concepts (in which energy and temperature are 
regarded as simultaneously determinable quantities) become meaningless, 
Statistical concepts, however, remain significant (Tolman 1938). 

The case of an ideal crystal was discussed many years ago (Planck 
1921, Schaefer 1921) and it was clear that this situation could in general 
be expected only at exceedingly low temperatures at present inaccessible. 
In rather exceptional cases (a cube of diamond with edges of 0-01 cm) 
these temperatures are within the range attainable by adiabatic 
demagnetization (Simon 1939). 
rN i eng 
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The case of thermal radiation enclosed in a cavity can be discussed in 
exactly the same way as an ideal crystal. This does not seem to have 
been done and it is the purpose of this note to fill this gap. 

Let us consider the electromagnetic radiation enclosed in a cavity in 
thermodynamic equilibrium at a temperature 7’. Treating this radiation 
in the usual way as a system of harmonic oscillators we can at once 
write down the partition function of the radiation. We have 


log Z= — log tl exp (Wyk) ieee “(TD 


The summation must be carried out over all the electromagnetic modes 
of vibration of the cavity, the energy is measured relative to the state 
of the system in which every oscillator is in its lowest state. 

With increasing quantum numbers the frequency density increases, 
ultimately becoming practically continuous. At sufficiently high 
temperatures these high states virtually determine the properties of the 
system. In this case it is possible to determine a continuous frequency 
density function (which is independent of the shape of the cavity) and 
the summation in (1) can be replaced by an integration. This leads to 
a manageable expression for log Z and the subsequent straightforward 
calculations of thermodynamic quantities can be found in any book on 
statistical mechanics. 

This procedure is certainly legitimate if for any pair of adjacent 


frequencies v’ and v” 
| A(v'—v")/kT | <1. ee ae ree ee (2) 


Let us consider the lowest possible frequencies in some details. These 
are of the order g/v!*; q is the velocity of propagation of the radiation 
in the enclosure which is of the order of the velocity of light in vacuo, 
and v is the volume of the enclosure. The differences v'—v” are also of 
the same order of magnitude. The actual values of the frequencies 
depend on the shape of the enclosure, but here we need consider only the 
order of magnitude. 
For these frequencies condition (2) becomes 


T > hg/kv, 


where hq/k ~ 1 cm degree. For a cavity with a volume of 1 cm* the usual 
procedure is certainly correct for temperatures much higher than 1° k. 
At temperatures of the order of 1° k or lower condition (2) is no longer 
fulfilled for the lowest frequencies. In this case, however, these frequencies 
virtually determine the properties of the system, as follows from the 
expression for their average number 7%, in a state 7 


j=1] {exp (hv,f{kT)—1}. 
For the lowest state of our example is of the order 1 and for higher 
states “i, decreases rapidly. The system behaves as if the number of 


degrees of freedom is small. Consequently thermodynamic considerations 
which are important for discussing the properties of thermal radiation at 
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high temperatures are no longer relevant, and typical results of the 
thermodynamic treatment such as Kirchhoff’s, Wien’s and Stefan’s law 
are no longer significant. All properties depend upon the shape of the 
enclosure, and, moreover, deviations from Planck’s law must be expected. 

An entirely similar analysis has been given by Planck and Schaefer for 
an ideal crystal. In this case the velocity g is about a factor 10° smaller. 
Therefore the temperature below which thermodynamic concepts break 
down is a factor 10° smaller than in the case of thermal radiation in an 
enclosure of the same size as the crystal. 


I am indebted to the Pressed Steel Company Limited for a Research 
Fellowship during the tenure of which this note was written. 
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(a) Decay event 3. The heavily ionizing particle from the lead plate decays in 
the gas and the charged secondary particle travels away from the camera 
out of the illuminated region. There is marked distortion of the tracks (see 
table 2). 


(b) Decay event 4. The primary particle comes from behind the cloud 
chamber and decays into a slow light meson in the lower half of the 
chamber. Interpretation as the decay of an upward moving V,° is also 
possible (see text), 
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(a) Decay event 6. A heavily ionizing negative particle decays after a flight in 
the chamber of 2-3x10-® sec. The secondary charged particle has a 


specific ionization <2-5 x minimum. 


i 


(b) Decay event 11. A lightly ionizing positive particle decays after a flight in 
the chamber of <4:1x10-!sec. The charged secondary has specific 
jonization <2°-5< minimum. 
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A. KORNDORFFER et al. Phil. Mag. Ser. 7, Vol. 43, Pl. LXXXII. 
Figs. 9-14 
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Figs, 15-20 . 
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